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Equilibrium and Non-equilibrium Quantum Field Theory 

Petr Jizba 
Summary 



This dissertation is concerned with various aspects of equihbriuni and non-equihbrium 
quantum field theory. 

We first focus in Chapter ^ on infrared effects in finite-temperature quantum field theory. 

^ . We propose a simple mathematical method (based on the largest -time equation and the 

Dyson-Schwinger equations) which allows systematic calculations of the change of density 
in energy /particles in heat-bath during a scattering/decay of external particles within the 
heat bath. A careful analysis reveals that the resulting changes in the energy density are 

Q^ ■ finite even in the case of massless heat -bath particles (no infrared catastrophe). 

^ ! As a next point we re-consider in Chapter |] the usual method of pressure calculations. 

O I We use the so called hydrostatic pressure (or pressure at a point) which is defined via the 

Q>^ ■ energy-momentum tensor. We go through all delicate points which one must deal with in 

0^ ■ the context of quantum field theory. Namely the renormalisation of composite operators 

and the vital role of renormalisation for a consistent quantum field theoretical definition of 
■ , pressure are discussed. We finally apply the whole procedure to a toy model system; A^"^ 

Qj \ theory with 0{N) internal symmetry. In the case of the large- A^ limit (also Hartree-Fock 

"-^ ■ approximation) the pressure is an exactly solvable quantity. Using the Mellin transform 

technique we then perform the large-temperature expansion of the pressure to all orders. 

r> I The hydrostatic pressure can be naturally extended to non-equilibrium systems. Us- 

c^ I ing the Jaynes-Gibbs principle of maximal entropy and the (non-equilibrium) Dyson- 

Schwinger equations we derive generalised Kubo-Martin-Schwinger equations and set up 
a calculational scheme for pressure calculations away from thermal equilibrium. As an ex- 
ample we explicitly evaluate in Chapter ^ the pressure for our 0{N) A$^ theory in the 
large- A^ limit in the case of two translationally invariant non-equilibrium systems. 

There follow five appendices which collect together much of the background material 
required in the main body of the thesis. The important part is the detailed analysis in 
Appendix |A.1| of the Dyson-Schwinger equations. The derivation there shows how the 



Dyson-Schwinger equations may be recast into a very useful functional form. In Appen- 
dices ^ and |C| we clarify some finer mathematical manipulations needed in Chapter ^. 
The fundamentals of the information or Shannon entropy are presented in Appendix 0. 
Appendix ^ covers the elements of dimensional regularisation and special functions which 
underlie much of the material presented in the earlier chapters. 
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Notation 



Natural units h = c = ks = 1 will be used throughout this dissertation. The following 
sub and superscripts will be used to label various quantities in the following text. 



Subscripts 



H 
C 

fix 

in 

hb 

r 

V 



Schrodinger picture 

Heisenberg picture 

Closed-time path (also the Keldysch-Schwinger path) 

Gauge-fixing term 

Interaction component 

Heat-bath component 

Renormalised object (operators, vertex functions, etc.) 

Volume 



Superscripts 



T 
L 
D 



Hermitian conjugation 

Complex conjugation (except for the T*-ordering) 
4-dimensional transverse components 
4-dimensional longitudinal components 
Symbol for a dimension 



Either 

/i,z/,A 
a, P 
a, b, c, d 



4-dimensional spacetime index 
3-dimensional space index 

Closed-time path indices :a,P = {+, — } for 
equilibrium, and a,P = {1,2} for equilibrium 
Internal symmetry indices 



non- 



Some of the important fields and functionals which will be used are. 



H= Jd^xHix) 


Hamiltonian 


C{x) 


Lagrange density 


P 


Projection operator onto final particle states 


§ 


§-matrix (§ = I + iT) 


T 


T-matrix 


Uit;t') 


Ket evolution operator 


P 


Density matrix 


Q^u 


General energy-momentum tensor 


er 


Canonical energy-momentum tensor 


$, 0, V 


General fields 



(n) 



Various other notations that will be used are. 

D Full thermal propagator for scalar fields in the closed- 

time formalism (i.e. 2x2 matrix) 

D'^ Connected thermal 2-point Green's function, i.e. connected part of . 

Dp Free thermal propagator for scalar fields in the closed- 

time formalism (i.e. 2x2 matrix) 
Full thermal n-point Green's function 

G^") Non-equilibrium n-point connected Green's function 

G Non-equilibrium propagator in the closed-time path 

formalism (i.e. 2x2 matrix) 

Jb Bose-Einstein distribution 

fp Fermi-Dirac distribution 

^ Information content 

/(. . .) Shannon (information) entropy 

S5(. . .) Amount of information conveyed by a single message 

M Bogoliubov matrix 

XXI Massieu function = — /3x (Elelmholtz) free energy 

0{z) Order z 

p{T) Thermodynamic pressure at the temperature T 

p{x, T) Hydrostatic pressure at the temperature T 

V{T) Hydrostatic pressure for translationally invariant media 

Sp Feynman (causal) propagator for spin-^ particles; (T = 0) 

S'^ Positive energy part of Sp 

S~ Negative energy part of Sp 

Sg von Neumann-Gibbs entropy 
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T Time ordering symbol 

T* The T* (or covariant) ordering 

W Generating functional for G^^^ 

Z (Grand) partition function 

Z The Jaynes-Gibbs partition function 

Z^ Wave function renormalisation 

P Inverse temperature 1/T 

r Generating functional for F*^"^ (i.e. the effective action) 

r*^") n-point vertex function (i.e. IPI n-point Green's function) 

■TTr-\ Left-handed variation 

j£(-^ Right-handed variation 

5^{. . .) Positive (+), negative (— ) frequency parts of Dirac's 

5-f unction 
Ai? Feynman (causal) propagator for scalar fields; (T = 0) 

A+ Positive energy part of Ap 

A~ Negative energy part of A^ 

jj, Chemical potential 

0^ Expectation value of the Heisenberg field operator <l>^ 

in the presence of a c-number source J (i.e. 0° = ($„)) 

E Proper self-energy for the scalar theory; (T = 0) 

S Proper self-energy for the scalar theory at finite tem- 

perature in the closed-time path formalism (i.e. 2x2 
matrix) 

S Self-energy for the scalar theory at finite temperature 

Q Grand (canonical) potential. In the case when the 

canonical ensemble is in question, Vt is (Helmholtz) free 
energy 

(...) Expectation value 

(. . .){pj.} Thermal expectation value. Particles with momenta 

{pk\ are unheated 



The Minkowski metric used throughout is g^^ = g^i, = diag(l, —1, —1, — 1) 
g.ug'"' = 5% = Trig) = 4, g,, = g^^ = -2 



Derivatives with respect to x'^ or x^ are abbreviated as 

d^ = d/dx'^, d^ = d/dXf^ 



m 



Chapter 1 

Introduction and overview 



The development of the theory of quantised fields at finite temperature and density (also 
thermal quantum field theory) over the past forty years or so has led to fundamental 
changes in the understanding of a wide number of physical phenomena. Among those one 
may mention symmetry restoration during high-temperature phase transitions [|l]-0] which 
has found significant application in early universe cosmology, and we may also mention 
applications in neutron-star [Q and supernova astrophysics. However, perhaps one of 
the most important applications of thermal quantum field theory nowadays is to quantum 
chromodynamics (QCD). 

The reason for this interest may be traced to the mid- 70 's when the notion of asymp- 
totic freedom of QCD started to emerge. At zero temperature and chemical potential the 
low-energy and/or momentum transfer behaviour of QCD is characterised by confinement 
(i.e. strong interaction among QCD constituents). The internal scale that determines the 
boundary between small and large energy in QCD is ~ ^qcd ~ 0.2GeV. As the en- 
ergy and/or momentum transfer increases, QCD begins to be characterised by asymptotic 
freedom, i.e., the coupling evolves as 
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{Q^ is the four-momentum transfer, Q^ > 0) and so quarks and gluons behave hke weakly 
interacting, massless particlesQ in high-energy and/or large Q^ processes. This behaviour 
is usually tested experimentally by means of deep inelastic electron-nucleon scattering. 

If one starts to study QCD at finite temperature and/or finite baryon density one 
automatically introduces new (intensive) variables, namely the temperature T and the quark 
chemical potential /i. These bring an additional mass scale with which the coupling Xqcd 
can run. It was Collins and Perry who first showed in 0] that strong interactions become 
weak not only at high energy-momentum transfer, as in the deep inelastic scattering, but 
also at very high baryon density. This reasoning was quickly extended to finite temperature 
[^,^ where it was shown that 

Agc7D(Q' = 0,T,/i)-.0, 

provided that T ^ ^qcd and/or yU ^ ^qcd- Thus at a sufficiently high temperature 
and/or baryon number density QCD systems consist of free quarks and gluons, regardless of 
the energy-momentum transfer. This "deconfined" phase of QCD is called the quark-gluon 
plasma. As was just mentioned, the temperature and/or chemical potential must be greater 
than the QCD fundamental mass scale Aqcd- In practice this means that the temperature 
for creation of the quark-gluon plasma must be at least of the order ~ 0.2GeV ~ lO^^K 

and/or the baryon number density must be of order ~ Aq^j^ ~ 0.8(GeV) ~ lO'^^cm"^. 

^Because quarks become massless, the deconfined phase leads to chiral-symmetry restoration (under 
normal circumstances the chiral flavour group SUl{N) x SUii{N), with N the number of quark flavours, 
is broken to a vector-like subgroup SUv{N)). The chiral phase transition is expected to be particularly 
interesting at high temperatures and/or densities m, uj . 
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Such temperatures prevailed in the very early stage of the universe (~ 10~^s — 10~^s after 
the Big-Bang) and such densities are expected to be present, for instance, in the core of 
neutron stars [Q. 

Apart from the cosmological and astrophysical implications, the quark-gluon plasma 
naturally offers an important laboratory for the study of asymptotic freedom. This fact has 
led to the construction of a Relativistic Heavy Ion Collider (RHIC) at Brookhaven National 
Laboratory, where it is expected that two on-head colliding beams of ^^'^Au will generate a 
sufficient centre-of-mass energy density which, when properly thermalized, will allow the 
formation of a quark-gluon plasma. Similar experiments are planned in the Large Hadron 
CoUider (LHC) at CERN with colliding Pb beams. 

So far we have only outlined possible applications of QFT to systems which are in 
thermodynamical equilibrium. However, during the last several years there has been a 
steadily growing awareness that the usual equilibrium (or thermodynamical) approaches 
in quantum field theory fail to describe correctly such fundamental phenomena as realistic 
phase transitions (both in condensed matter physics and in high-energy physics), quantum 
electrodynamical (QED) plasma evolution (and the related problem of hot thermonuclear 
fusion), heat transfer in stars, etc. Presently, the most important applications of non- 
equilibrium QFT, however, are in the physics of the early universe. The reason why it is 
so important to go beyond conventional equilibrium approaches in describing the universe 
is the extreme conditions present in the early period of its evolution 0. It is now well 
recognised that at the very early times the universe was very hot with energy densities 
approaching the Planck mass density ~ 5.2 x 10^^ g cm~^, and it then cooled rapidly due 
to the expansion. Because of this fact, it is expected that the matter fields went through 
a hierarchy of phase transitions before reaching the present status quo ||TU|. It is obvious 
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that if the matter fields involved in the phase transitions interacted at a rate much smaller 
than the typical cooling (or relaxation) rate of the universe, the evolution would proceed in 
an out-of-equilibrium manner. Such non-equilibrium dynamics may, in turn, lead to many 
crucial physical consequences. One such example is the baryon asymmetry [P, |lT], |12| . If 



one assumes that the relative abundance of baryonic matter over antibaryonic matter is not 
created a priori via the initial-time conditions, then one must find a physical mechanism 
which could generate such an asymmetry. It was A.D. Sakharov [0] who proposed more 
than thirty years ago that, in order to explain a biased production of baryons over anti- 
baryons, one needs to take into account the non-equilibrium evolution of the universe since 
baryons and anti-baryons are produced in equal number in any equilibrium process. This 
is one of the celebrated Sakharov criteria [JTT],0,[T^ for baryogenesis. 

Another extremely interesting example of application of non-equilibrium QFT in cos- 
mology is the production of topological defects during phase transitions. Interest in this 
area stems from the belief that topological defects, such as cosmic strings, might provide 
an explanation of structure formation and the cosmic microwave background radiation 
anisotropies in the universe |jl5|. The domain structure of a ferromagnet is known to be- 



come very different when the sample is cooled adiabatically through the Curie temperature 
compared to when it is cooled rapidly. In analogy one may expect that the production 
of topological defects and their evolution will depend strongly on the cooling rate of the 
universe. The formation of defects and their dynamics can be beautifully mimicked in 
quantum liquids [p!6|,p!7| such as helium ^He or ^He, in liquid nematic crystals |]18[ or in 



superconductors |^ . Recent experiments [jT6[ with ^He have confirmed that the density of 
vortex lines nucleated during the phase transition from a normal ^He liquid to its superfluid 
B-phase depends considerably on the cooling rate. 
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In the framework of non-equilibrium applications of QFT one should not forget the 
currently very important application to relativistic heavy-ion collisions, which, as we have 
mentioned, are expected to produce a sufficient thermal environment for creation of a 
quark-gluon plasma. The major indication that a non-equilibrium description seems to 
be necessary comes from the expectation that the time scale at which the quark-gluon 
plasma should exhibit itself will be too short for its macroscopic equilibration. The point 
is that the energy density in the reaction zone, once the plasma begins to thermalise, will 
create a pressure which will consequently lead to explosive expansion of the plasma. The 
associated expansion time is estimated to be only 10 or 100 times longer than the time 
scale of microscopic equilibration processes ||20|. This introduces a basic uncertainty as 
to whether an equilibrium treatment may be used safely. On the other hand, if a non- 
equilibrium description proves crucial, then one must face the question to what extent a 
small, short-lived, fast expanding system of quarks and gluons can be called a quark- 
gluon plasma. One should honestly admit that there does not presently exist any generally 
accepted model which would satisfactorily describe the quark-gluon plasma dynamics. 

All the mentioned non-equilibrium processes are characterised by the rate at which a 
system changes (e.g. expansion rate, dissipative rate, the rate at which particles are ex- 
changed with an external environment) which is comparable to or greater than the rate 
at which microscopic interactions (i.e. the equilibrating mechanism) are happening. In 
practice this means that the relaxation time scale (i.e. the mean time in which the sys- 
tem relaxes to equilibrium) is much longer than the time scale at which observations are 
performed, and consequently a system does not evolve adiabatically (i.e. it does not go 
through a sequence of states each of which is at thermal equilibrium). It is needless to say 
that many of phenomena in Nature are precisely of this character. 
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Whereas for the description of matter in equihbrium one has at one's disposal a sys- 
tematic and unified approach based upon the formahsm of the Gibbs (micro- or grand-) 
canonical density matrices, no such simple way seems to exist in a non-equilibrium theory 
because of the variety of phenomena and of the complexity of the evolution processes. This 
is to be expected, as there is basically no limitation on the amount of questions that may 
be asked about time-dependent phenomena unless one clearly specifies which degrees of 
freedom (or what degree of reduction) one will adopt for the description of the evolution 
of a system. Within such a reduced description, however, the dynamics has been shown to 
be capable of prediction. In this framework there have steadily crystallised two major (and 
mutually distinct) approaches [plH 23[|. 

The first one, which we mostly refer to in our thesis, is the, so called, 'Maxent school' 
(or maximal entropy school) founded by E.T.Jaynes et al. p^-p7|. The basic philosophy 
behind this approach is rooted in the fact that the prediction of the future macroscopic 
evolution of a system cannot be done with certainty on the basis of the initial macroscopic 
data because of the existing correlations with the parameters or data that are discarded 
in the given reduced description. The corresponding statistical inference about the system 
is not deduced from the underlying microscopic dynamics but instead is rather based on 
information theory. There the algorithm of entropy maximalisation leads to the density 
matrix (or probability distribution) with the least informative content subject to the prior 
knowledge which one has about the system. Generally one may, at different times, adopt 
different macroscopic parameters describing the system (i.e. at different times one may 
choose a different level or reduction in the description). If this is the case, the maximalisa- 
tion of entropy naturally leads to a more complex form of the density matrix. This branch 
of the Maxent school (and its various modifications) is also known as the projection operator 
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technique 

The second school, the so called 'Brussels school', is associated with I.Prigogine, R.Bales- 
cu et al. [^, ^,|3^. The basic emphasis is here put on the microscopic dynamics (i.e. 



Hamilton, Schrodinger or Liouville equations), and all other non-dynamical approaches, 
such as those coming from information theory, are discarded: everything should be derived 
from the dynamics alone. The ultimate aim is a dynamical separation into the parameters 
to be retained on the desired level of description, and those to be discarded. This separation, 
if it exists, should emerge the from dynamics as an asymptotic property valid for the large 
(i.e. observational) time scales. It is well known since the time of Boltzmann that the latter 
may be achieved by introducing a certain hypothesis about the microscopic behaviour of the 
system; e.g. Boltzmann's "Stosszahlansatz" (random collision hypothesis) |31|], Ehrenfest's 
coarse graining hypothesis [0, or Bogoliubov's no initial-time correlations hypothesis [^ 



This line of reasoning is basically inherited and progressed by the Brussels school. It should 
be noted that all the statistical inferences or hypotheses here have a strictly dynamical 
nature (they are directly motivated by the underlying dynamics), but on the other hand 
they are in a certain respect ad hoc, because only certain probabilistic features of the 
microscopic dynamics are emphasised. 

Synopsis 

We first focus in Chapter ^ on infrared effects in finite-temperature QFT. We propose a 
simple mathematical method (based on the largest-time equation and the Dyson-Schwinger 
equations) which allows systematic calculations of the change of energy density (or particle 
density) in a heat-bath during a decay (or scattering) of the external particle (s) within 
the heat bath. The applied method naturally leads to an interpretation of the change of 
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density of energy (particles) in terms of three additive contributions: stimulated emission, 
absorption and fluctuations of the heat-bath particles. This result is completely non- 
perturbative. A careful analysis reveals that the resulting change in the energy density is 
finite even in the case of massless heat-bath particles. This means that there is no infrared 
catastrophe. 

As the next point we re-consider in Chapter |^ the problem of calculating pressure. 
We use the so-called hydrostatic pressure (or pressure at a point) which is defined via 
the energy-momentum tensor. The obvious advantage is a possible extension into a non- 
equilibrium medium. We go through all the delicate points that must be dealt with in the 
context of quantum field theory. Renormalisation of composite operators and in general 
the vital role of renormalisation for a consistent quantum field-theoretical definition of 
pressure is discussed. We finally apply the whole procedure to a toy-model system: A$^ 
theory with 0{N) internal symmetry. In the case of the large- A^ limit (also the Hartree- 
Fock approximation) the pressure is exactly solvable. Using the Mellin transform technique 
we perform the large-temperature expansion of the pressure to all orders in T. 

The hydrostatic pressure can be naturally extended to non-equilibrium systems. Us- 
ing the Jaynes-Gibbs principle of maximal entropy and the (non-equilibrium) Dyson- 
Schwinger equations we derive in Chapter ^ the generalised Kubo-Martin-Schwinger equa- 
tions and set up a calculational scheme for pressure calculations away from thermal equi- 
librium. As an example we explicitly evaluate pressure for the 0{N) A^^ theory in the 
large- A^ limit in two cases of translationally invariant non-equilibrium systems. 

There follow five appendices which comprise much of the background material required 
in the main body of the thesis. The important part is the detailed analysis in Appendix |A.1| 
of the Dyson-Schwinger equations. The derivation there shows how the Dyson-Schwinger 
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equations may be formulated in a very useful functional form. We also outline the connec- 
tion with the more conventional approaches. In Appendices |B| and |C| we clarify some finer 
mathematical manipulations needed in Chapter §. The fundamentals of information or 
Shannon entropy are presented in Appendix 0. Appendix ^ covers the elements of dimen- 
sional regularisation and special functions which underlie much of the material presented 
in the earlier chapters. 

Epilogue 

A reader of this dissertation might be disappointed by the fact that he or she will not find 
a usual pedagogical introduction to the subject. The omission of such an introduction was 
dictated mainly by two considerations. Firstly, we have not felt very competent to provide 
a good account of the fundamentals of both equilibrium and non-equilibrium quantum 
field theory. The subjects themselves are currently immensely vast and the number of 
problems involved quickly approaches the "thermodynamic limit". We therefore take a 
more pragmatic point of view and proceed by sampling a few definite problems, which 
we develop and analyse in great detail. Secondly, and most importantly, we personally 
hold the opinion that a dissertation should reflect student's ability to cope with a subject 
and creatively apply it to practical problems rather than write an essay on mathematical 
or physical foundations of the subject in question. Any such attempts would lead (at 
least in our case) to pure epigonism since we do not feel that we could add anything 
substantially new to existing (and pedagogically excellent) textbooks and review articles; 



see for example Refs. [p^,p3-H2[ for equilibrium quantum field theory and statistical physics 



and Refs. |Q^,^,|^ for non-equilibrium quantum field theory and statistical physics. 
To end, we wish to make one more remark. We are perfectly aware that the presented 



CHAPTER 1. INTRODUCTION AND OVERVIEW 10 

work is incomplete in many respects. Our ignorance or lack of understanding of many 
important topics is of course in part responsible for this weakness. In particular, we refer 
here to discussions and applications of such crucial areas as the imaginary-time (Euclidean 
or Matsubara) formalism, (non-equilibrium) thermofield dynamics, hard thermal loops, 
the theory of temperature-induced phase transitions, linear response theory, stochastic 
approaches to non-equilibrium dynamics, quantum transport equations, theory of transport 
coefficients and many more, which are definitely missing in this dissertation. Although we 
had originally planned to include some of the aforementioned issues (namely those which 
directly concern or resonate with our present research), neither space nor time has allowed 
us to fulfil this wish. 



Chapter 2 

Heat— bath particle number spectrum 



In recent years much theoretical effort has been invested in the understanding relativistic 
heavy ion colhsions as these can create critical energy densities which are large enough to 
produce the quark-gluon plasma (the deconfined phase of quarks and gluons) @,^ . 



A natural tool for testing the quark-gluon plasma properties would be to look for 
the particle number spectrum formed when a particle decays within the plasma itself. 
As the plasma created during heavy ion collisions is, to a very good approximation, in 
thermodynamical equilibrium [Q (somewhat like a "microwave oven" or a heat bath), one 
can use the whole machinery of statistical physics and QFT in order to predict the final 
plasma number spectrum. Such calculations, derived from first principles, were carried out 
by Landshoff and Taylor [^ . 

Our aim is to find a sufficiently easy mathematical formalism allowing us to perform 
mentioned calculations to any order. Because unstable particles treated in |^ can not 
naturally appear in asymptotic states, we demonstrate our approach on a mathematically 
more correct (but from practical point of view less relevant) process; namely on the scat- 
tering of two particles inside of a heat bath. The method presented here however, might 
be applied as well to a decay itself (provided that the corresponding decay rate is much 
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less than any of the characteristic energies of the process). In this chapter we formulate 
the basic diagrammatic rules for the methodical perturbative calculus of plasma particle 
number spectrum d{N {uj)) / du and discuss it in the simple case of a heat bath comprised 
of photons and electrons, which are for simplicity treated as scalar particles. 



In Section ^]T| we review the basic concepts and techniques needed from the theory of 
the largest-time equation (both for T = and T 7^ 0) and the Dyson-Schwinger equations. 
Rules for the cut diagrams at finite-temperature are derived and subsequently extended to 
the case when unheated fields are present. It was already pointed out in [BH] that the thermal 



cut diagrams are virtually the Kobes-Semenoff diagrams |]3^ in the Keldysh formalism [^ 
This observation will allow us to identify type 1 vertices in the real time finite — temperature 
diagrams with the uncircled vertices used in the (thermal) cut diagrams, and similarly type 
2 vertices will be identified with the circled, cut diagram vertices. As we want to restrict our 
attention to only some particular final particle states, further restrictions on the possible 
cut diagrams must be included. We shall study these restrictions in the last part of Section 



2.1 



As we shall show in Section |2.2| , the heat-bath particle number spectrum can be conve- 
niently expressed as a fraction. Whilst it is possible to compute the denominator by means 
of the thermal cut diagrams developed in Section p.l| , the calculation of the numerator 



requires more care. Using the Dyson-Schwinger equations, we shall see in Section 2]3| that 
it can be calculated through modified thermal cut diagrams. The modification consists of 
the substitution in turn of each heat bath particle propagator by an altered one. We also 
show that there must be only one modification per diagram. From this we conclude that 
from each individual cut diagram we get n modified ones (n stands for the total number of 
heat-bath particle propagators in the diagram). Furthermore, in the case when more types 
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of the heat bath particles are present, one might be only interested in the number spectrum 
of some of these. The construction of the modified cut diagrams in such cases follow the 
same procedure as in the previous situation. We find that only the propagators affiliated 
to the desired fields must be altered. 

In Section p.4| the presented approach is applied to a toy model in which a gluon plasma is 
simulated by scalar photons, and we calculate the resulting changes in the number spectrum 
of the 'plasma' particles. Section |2]^ ends with a qualitative discussion of the quark-gluon 
plasma simulated by scalar photons and electrons. 

For reader's convenience, the chapter is accompanied with Appendix |A.1| where we 



derive, directly from the thermal Wick's theorem, the (thermal) Dyson-Schwinger equations 
as well as other useful functional identities valid at finite temperature. 

2.1 Basic tools 

2.1.1 Mean statistical value 

The central idea of thermal QFT is based on the fact that one can not take the expectation 
value of an observable A with respect to some pure state as generally all states have non- 
zero probability to be populated and consequently one must consider instead a mixture of 
states generally described by the density matrix p. The mean statistical value of A is then 

{A)=Tr{pA), (2.1) 

where the trace has to be taken over a complete set of physical states. For a statistical system 
in thermodynamical equilibrium p is given by the Gibbs (grand-) canonical distribution 
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P^TH^Fm^J^^~Z~' ^^'^^ 

here Z is the partition function, H is the Hamiltonian, A^ is the conserved charge (e.g. 
baryon or lepton number), /i is the chemical potential, K = H — fiN, and (3 is the inverse 
temperature: /3 = 1/T. 

2.1.2 Largest— time equation at T=0 

An important property inherited from zero-temperature QFT is the largest-time equation 
(LTE) ||48|-^. Although the following sections will mainly hinge on the thermal LTE, 
it is instructive to start first with the zero-temperature one. The LTE at T = is a 
generally valid identity which holds for any individual diagram constructed with propagators 
satisfying certain simple properties. For instance, for the scalar theory with a coupling 
constant g one can define the following rules: 

X i ] y ~ i Ap (x-y) 

X i ? y ~ i A" (x-y) 

X ?: ] y ~ i A+ (x-y) 

X 1 ? y ~ -i A*p (x-y) 

Here iAp is the Feynman propagator, iA'^ i'iA~) is corresponding positive (negative) energy 
part of iAp, the '*' means complex conjugation and index 1 (2) denotes type-1 {type-2) 
vertex; type-1 vertex has attached a factor —ig whilst type-2 bears a factor ig. Using this 
prescription, we can construct diagrams in configuration space. With each diagram then 
can be associated a function F{xi, . . .x„) having all the 2nd type vertices underlined. For 
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example, for the triangle diagram in FIG.|2.1j we have 



15 




Figure 2.1: A one loop triangle diagram. 



F{xi,X2,X3) 
F{xj_,X2,X3) 
F{xj_,X2,X3) 
F{xj_,X2,X3) 

etc. 



{-igY i^F{xi - X2) iAF{xi - X3) iAF{x2 - X3) 
{-igf{ig) iA+(xi - X2) iA+(xi - X3) iAir(x2 - X3) 
{igf{-ig) (-i)A^(xi - X2) iA+(xi - X3) 2A+(x2 - X3) 
{igf (-z)A^(xi - X2) (-2)A^(xi - X3) (-i)A^(x2 - X3) 



The LTE then states that for a function -F(xi, . . . , x„) corresponding to some diagram with 
n vertices 



I Xi 



+ F{ 



, Xn 



(2.3) 



provided that Xjo is the largest time and all other underlinings in F are the same. The proof 
of Eq.( p.3|) is based on an observation that the propagator i/S.p{x) can be decomposed into 
positive and negative energy parts, i.e. 
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iApix) = e{xo)iA'^{x) + e{-xo)iA~{x), (2.4) 

^A±(x) = j-0-e-^'^e{±ko)6{e-m'). (2.5) 

Incidentally, for Xio being the largest time this directly implies 



iApixj — Xi) = iA {xj — Xi), 

iApiO) = -iA*p{0). (2.6) 

As F{. . . ,Xi, . . .) differs from F{. . . ,Xi, . . .) only in the propagators directly connected to 
Xi - which are equal (see Eq.( |2.6| )) - and in the sign of the Xi vertex, they must mutually 
cancel. 



Summing up Eq. (|2.3|) for all possible under linings (excluding a;,), we get the LTE where 



the special role of the largest time is not manifest any more, namely 



^F(xi,X2,...,x„) = 0. (2.7) 



index 



The sum ^j^^g^. means summing over all possible distributions of indices 1 and 2 (or equiva- 
lently over all possible underlinings) . The zero-temperature LTE can be easily reformulated 
for the T-matrices. Let us remind that the Feynman diagrams for the §-matrix (S = I+zT) 
can be obtained by multiplying the corresponding F{xi, . . . , x„) with the plane waves for 
the incoming and outgoing particles, and subsequently integrate over xi . . .x„. Thus, in 
fixed volume quantisation a typical n-vertex Feynman diagram is given by 



IT ^^»ri . \\ , F(xi,...,Xn). (2. 
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Here the momenta {pj} are attached to incoming particles at the vertices {xm }, while mo- 
menta {qk} are attached to outgoing particles at the vertices {xm,,}- In order to distinguish 
among various functions F{xi, . . . , a:„) with the same variables Xi, . . . , x„, we shall attach 
a subscript /„ to each function F. For instance, the function Fi^{xi, . . . ,0:4) corresponding 
to the diagram 




■4 ^ ^2 



contributes to {qiq2\i"^\'Pi'P2) by 



I I dx, 



Q-i(pi+P2)xi Qi(qi+q2)x4 



■.{iApixi - X2)f{iAF{x2 - X3)f{iAF{x3, - X4)y 



similarly, the function ^24(0:1, . . . ,^4) corresponding to the diagram 




contributes to {qiq2\i'^\PiP2) by 



]^(ia:j 



-i(pi+P2)xi Qi{q\+q2)x4, 



4 = 1 



V^Aup^Up., V^/Auj, 



qi'^q2 



■.iAp{xi - X2)iAF{xi - X3){iAp{x2 - xs))" 



X iAp{x4 — X3)iAp{x4^ — X2), 



etc. 
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This can be summarised as 

Consider now the case \{pj]) = \{<lk}) (let us call it \a)). From the unitarity condition: 

T - Tt = iJ^T, we get 

(a|T|a) - (a|T|a)* = i{a\T^T\a). (2.10) 



On the other hand, by construction F{xi, . . . , x„) = F*{xi, . . . , x„), and thus (see ( |2.7] )) 



F(xi, ...,Xn) + F*{xi,...,Xn) = - ^ F{xi, . . . , x„). (2.11) 



index 



The prime over index in ( |2.11| ) indicates that we sum neither over diagrams with all type 



1 vertices nor diagrams with all type 2 vertices. Using ( |2.9| ), and identifying |{q'fc}) with 
\{Pk}) (= \a)) we get 



(a|T|a) - {a\T\a)* = - ^ (a|T|a), (2.12) 

index 



or see 



(ES)) 



(a|T"fT|a) =i ^ (a|T|a). (2.13) 

index 



Eq. ( p.l2|) is the special case of the LTE for the T-matrices. The finite-temperature exten- 
sion of (p.l3|) will prove crucial in Section |2.3| . 



Owing to the 6'(±fco) in A^(a;) (see Eq.( |2.5| )), energy is forced to flow only towards type 
2 vertices. From both the energy-momentum conservation in each vertex and from the 
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energy flow on the external lines, a sizable class of the diagrams on the RHS's of ( p.l2| )- 
( |2.13| ) will be automatically zero. Particularly regions of either 1st or 2nd type vertices 
which are not connected to any external line violate the energy conservation and thus 
do not contribute (no islands of vertices), see FIG . p .21 . Consequently, the only surviving 
diagrams are those whose any 1st type vertex area is connected to incoming particles and 
any 2nd type vertex area is connected to outgoing ones. From historical reasons the border 
between two regions with different type of vertices is called cut and corresponding diagrams 
are called cut diagrams. 





Figure 2.2: An example of a cut diagram in the ip^ theory which does not contribute to the 
RHS's of ( p.l2D -( 2.13D . Arrows indicate the flow of energy. 



We have just proved a typical feature of T = QFT, namely any cut diagram is divided 
by the cut into two areas only, see FIG. |2.3| . Eq. (|2.12|) , rewritten in terms of the cuttings is 



so called cutting equation (or Cutkosky's cutting rules) p8|-[50[ . 




Figure 2.3: Generic form of the cut diagram at the T = 0. Shadow is on the 2nd type vertex 
area. 



One point should be added. Inserting the completeness relation J2f I/) (/I = 1 i^to the 
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LHS of (|2J3| ), we get 

5^(a|Tt|/)(/|T|a) = z5^(a|T|a). (2.14) 

/ cuts 

It may be shown [^,^ that all the intermediate particles in |/) correspond to cut lines. 



This has a natural extension when (a|T'''T|a) -^ (a|T"''PT|a) with P being a projection 
operator (P = pt = P^) which eliminates some of the states |/). It is easy to see that in 
such case 

{a\T^FT\a) = i^{a\T\a), (2.15) 

cuts 

where tilde over the Ylcuts indicates that one sums over the diagrams which do not have 
the cut lines corresponding to particles removed by P. 

There is no difficulty in applying the previous results to spin-| |^, ^. The LTE 



follows as before: the diagram with only iSp propagators (and —ig per each vertex) plus 
the diagram with only (iSp) propagator Q (and ig per each vertex) equals to minus the 
sum of all diagrams with one up to n — 1 the type 2 vertices (n being the total number of 
vertices) . For gauge fields more care is needed. Using the Ward identities one can show [^ 
that type 1 and type 2 vertices in ( p^.l2D -( p.l3D may be mutually connected only by physical 



particle propagators (i.e. neither through the propagators corresponding to particles with 
non-physical polarisations or Fadeev-Popov ghosts and antighosts). 



^The function iSpix), similarly as (iAi?)*(a;), interchanges the role 5*+ and S . Unlike bosons, for 



fermions iSpix) is not equal to (iSpyix). Despite that, Eq.(2.12) still holds |48 
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2.1.3 Thermal Wick's theorem (the Dyson— Schwinger equation) 

The key observation at finite temperature is that for systems of non-interacting particles 
in thermodynamical equihbrium Wick's theorem is still valid, i.e. one can decompose the 
2r2-point (free) thermal Green function into a product of two-point (free) thermal Green 
functions. This may be defined recursively by 

{T{^{x^) . . . ^(X2„))) = J2 ^pini^ixMxj))) (T( n H^k))), (2.16) 

where ep is the signature of the permutation of fermion operators (= 1 for boson operators) 
and T is the standard time ordering symbol (for generalisation to the contour-time-path 
ordering see Appendix |A.1| ). Note that the choice of "i" in ( p.l6|) is completely arbitrary. 



The proof can be found for example in |^,^|5^. Similarly as at T = 0, Wick's theorem 

i.e. 



can also be written for the (free) thermal Wightman functions ||51|,p3 



{ij{x,) . . .ij{x2n)) =J2ep{ij{x,)ij{x,)) iYlH^k)). (2.17) 






A particularly advantageous form of this is the so called Dyson-Schwinger equation (see 
Appendix |A.1| ) which, at the T 7^ 0, reads 



{G[^Mx)Fm = I dzmx)^{z)) ( G[^]^^\ + J dz{^{z)^{x)) /^^F[^P] 



f2.18) 



where ip{x) is an interaction-picture field and G[. . .] and F[. . .] are functionals of ip. The 
arrowed variations jtt^ are defined as formal operations satisfying two conditions, namely: 
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or 



dlpniz) Olpniz) 0^n{z) 

with 

. , . X =d{x-y)dmn- (2.21) 

(>Vn{y) 

The "p" is for bosons and 1 for fermions; subscripts m, n suggest that several types of 
fields can be generally present. Note, for bosons -^ = -|t- which we shall denote as |^. 



For more details see Appendix |A.l . 



2.1.4 Thermal largest— time equation 



The LTE (|2.13| ) can be extended to the finite-temperature case, too. Summing up 



m 



(|2.13| ) over all the eigenstates oi K {= H — fiN) with the weight factor e ^^' {i labels the 
eigenstates), we get 

(TT"f) =i ^ (T). (2.22) 

index 



Let us consider the RHS of ( |2.22| ) first. The corresponding thermal LTE and diagrammatic 



rules (Kobes-Semenoff rules [|^) can be derived precisely the same way as at T = using 



the previous, largest-time argumentation |53,E5|. It turns out that these rules have basically 



identical form as those in the previous section, with an exception that now (0| . . . |0) -^ {■■■)■ 
Note that labelling vertices by 1 and 2 we have naturally got a doubling of the number of 
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degrees of freedom. This is a typical feature of the real-time formalism in thermal QFT 
(here, in so called Keldysh version p4|). 

We should also emphasise that it may happen some fields are not thermalized. For 
example, external particles entering a heat bath or particles describing non-physical degrees 



of freedom [^. Particularly, if some particles (with momenta {pj}) enter the heat bath. 



the mean statistical value of an observable A is then 



g-/3/fi 



i 

Ptel = lfe})(fe}l. 



which we shall denote as (A) {p.}. From this easily follows the generalisation of ( |2.22| ) 



(TTn{p.}=^EWiP^-}- (2-23) 



index 



Unlike T = 0, we find that the cut diagrams have disconnected vertex areas and no kine- 
matic reasonings used in the previous section can, in general, get rid of them. This is 
because the thermal part of {ip{x)ip{y)) describes^ the absorption of on shell particle from 
the heat bath or the emission of one into it. Thus, at T 7^ 0, there is no definite direction 
of transfer of energy from type 1 vertex to type 2 one as energy flows in both directions. 
Some cut diagrams nevertheless vanish. It is simple to see that only those diagrams survive 
in which the non-thermalized external particles "enter" a diagram via the 1st type vertices 
and "leave" it via the 2nd type ones. We might deduce this from the definition of (T){p }, 
indeed 



2Note that {^{x)ip{y)) = (: ^(x)^(y) :) + (0|^(a;)(^(y)|0) and (: ^{x)ip{y) :) = /^/5(fco)5(fc2 

2y-ik(x-y)^ with /^(fco) = (e'^l'^ol _ l)"!. 
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E Wtei = E E ^(^' fe>iTi^' fe»- (2-24) 

index index * 



Note, we get the same set of thermal cut diagrams interchanging the summation J2index' 
with J2i- It is useful to start then with ^^i^dex' (^! {Pill'^K; {Pj})- This is, as usual, described 
by the (T = 0) cutting rules. In the last section we learned that the general structure of the 
corresponding cut diagrams is depicted in FIG. p.3| , particularly the external particles enter 
the cut diagram via type 1 vertices and leave it via type 2 ones. Multiplying each diagram 
(with the external particles in the state \i; {pj})) with the pre-factor ^ ^ and summing 
subsequently over i, we again retrieve the thermal cut diagrams, though now it becomes 
evident that the particles {pj} enter such diagram only via type 1 vertices and move off 
only through type 2 ones, since the summation of the (T = 0) cut diagrams from which 
it was derived does not touch lines corresponding to unheated particles. Note, the latter 
analysis naturally explains why the unheated particles obey the (T = 0) LTE diagrammatic 
rules even in the thermal diagrams 

Another vanishing comes from kinematic reasons. Namely three-leg vertices with all 
the on shell particles (1-2 lines) can not conserve energy-momentum and consequently 
the whole cut diagram is zero. As an illustration let us consider all the non- vanishing, 
topologically equivalent cut diagrams of given type involved in a three-loop contribution to 
i Yliindex ^)vQ (^^^ FIG.|2.4|)P|. Let us stress one more point. In contrast with T = 0, at finite 
temperature the cut itself neither is unique nor defines topologically equivalent areas, see 
FIG. p.5| , only the number of crossed legs is, by definition, invariant. This ambiguity shows 
that the concept of the cut is not very useful at finite temperature and in the following we 



^Let us emphasise that originally we had the 64 possible cut diagrams. 
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Figure 2.4: An example of non-vanishing cut diagrams at the T ^ 0. The heat -bath consists of 
two different particles. External particles are not thermalized. 

shall refrain from using it. 

In Section ^]3| it will prove useful to have an analogy of (|2.23| ) for {T'^F T) . Here P 







Figure 2.5: The cut diagram from Fig j2.3| c) demonstrates that the cut can be defined in many 
ways but the number of crossed lines is still the same. 



is the projection operator defined as 



P = $^l«;j)(«;j| 



(2.25) 



where "j" denotes the physical states for the heat-bath particles and "a" labels the physical 
states for the outgoing, non-thermalized particles. Let us deal with (T"''PT). Using (|2.15| ), 
we acquire 



PKi 



(TtpT) = z5^^5^(/|T|0 



Interchanging the summations, we finally arrive at 



(2.26) 
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(T¥T)=i^(T), (2.27) 

index 

where tilde over the J2index' Cleans that we are restricted to consider the cut diagrams, with 
only (l-2)-particle lines corresponding to the "a" and "j"particles (i.e. {0\(f{x)(f{y)\0) and 
{il){x)ip{y)) , respectively). The extension of Eq.( |2.27| ) to the case where some external, 
non-thermalized particles {pk} are present is obvious, and reads 

(TtpT)|,,} = z^(T)|,,}. (2.28) 

index 

Finally, let us note that using the LTE, one may extend the previous treatment to various 
Green functions. The LTE for Green's functions is then a useful starting point for dispersion 



relations, see e.g. |34, 54 . 



2.2 Heat— Bath particle number spectrum: 
general framework 



The cutting equation ( |2.28| ) can be fruitfully used for both the partition function Z and 



the heat-bath particle number spectrum d{N {uo)) / duo calculations. To see that, let us 
for simplicity assume that two particles (say $1, $2) scatter inside a heat bath. We are 
interested in the heat-bath number spectrum after two different particles (say 0i, ^2) appear 
in the final state. Except for the condition that the external particles are different from the 
heat bath ones, no additional assumption about their nature is needed at this stage. 

The initial density matrix pi (i.e. the density matrix describing the physical situation 
before we introduce the particles $i(pi), $2(^2) into the oven) can be written as 
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P^ = ^r'5Ze-'^^^|j;pi,p2)(j;Pi,P2|, (2.29) 

j 

where "j" denotes the set of occupation numbers for the heat-bath particles. A long time 

after the scattering the final density matrix pj reads 

Pf = Zj'J2e~^''^FS\j;pi,p2){j;puP2\^^^\ (2.30) 

j 

here P is the projection operator projecting out all the non-heat-bath final states except of 
4>i{qi) , 4'2{(l2) ones. The §-matrix in (|2.3CI|) is defined in a standard way: S =l + iT. The 



Zf in (|2.3CI| ) must be different from Zi as otherwise pf would not be normalised to unity. 
In order that pf satisfy the normalisation condition Tr{pf) = 1, one finds 

Zf = 5^e-^^^(j;pi,p2|S^P§|j;Pi,P2) = (§^PS)p,p. Z, = (TtpT)^,^, Z,. (2.31) 



The key point is that we have used in ( p.31|) the T-matrix because the initial state 



|$i(pi), $2(^2)) is, by definition, different from the final one \(p{qi), 02(^2)) and consequently 
P § can be replaced by iF T. This allows us to calculate Zf using directly the diagrammatic 
technique outlined in the preceding section. 

From ( p.l|) and (|2.3(]| ) one can directly read off that the number spectrum of the heat- 



bath particles is: 



^-^^^ = /|^^"(-^-k^--nE(/l«Kk;-)aKk;c.)p,|/) 



-^^S^u.^ - k^ - 2)(TtPa;(k;^)aKk;.;)T),,, 
(27r)3 ^ '^ (TtPT)p,p2 
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and consequently 

where we have used the completeness relation for the final states |/) and [P; a^a] = 0. 
The subscript "/" denotes which type of heat-bath particles we are interested in. In the 
following the index will be mostly suppressed. 

2.3 Modified cut diagrams 



To proceed further with ( p. 32 ) and (|2.33|) , we expand the T-matrix in terms of time-ordered 



interaction-picture fields, i.e. 



I'M = ^ dxi... dxnanixi . . . Xn)T{^{xi) . . . i){xn)). (2.34) 

Here ip represents a heat-bath field in the interaction picture. Other fields (i.e. 0,0 and 
$) are included^ in the a„. An extension of ( 2.34 ) to the case where different heat-bath 



fields are present is natural. Employing ( p. 34 ) in (T"''PT)p^p2, one can readily see that this 



f actor ises out in each term of the expansion a pure thermal mean value (...). The general 
structure of each such thermal mean value is: (Gm['0]-^n['0]), where F„[. . .] and Gm[- ■ •] are 
the operators with "n" chronological and "?ti" anti-chronological time ordered (heat-bath) 
fields, respectively. Analogous factorisation is true in the expansion of {T'^¥'^a^aFT)p^p^. 
The only difference is that the pure thermal mean value has the form {Gm[i'\a'^a'Fn 



"* When Fermi fields are involved, we have, for the sake of compactness, included in the argument of ?/> 
the space-time coordinate, the Dirac index, and a discrete index which distinguishes tpa from ip^. 
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instead^. In case when various heat-bath fields are present, m = mi + 1712 + . . . + rrin, with 
"m;" denoting the number of the heat-bath fields of /-th type. 

Applying the Dyson-Schwinger equation to (G'm[^]a"'^aF„[-?/']) twice and summing over 
"n" and "m", we get cheaply the following expression (c.f. also (|A.11|) ) 






PlP2 



\ / piP2 



2 ^^ '"'^ '" '"'^"^' ^ ' ' '"'^ '" '"'^"^'^ \ 6My)mx) 



P1P2 



+ {a]ai){rFT)p,p,, (2.35) 

A similar decomposition for {T'^F T) p^^p.^ would not be very useful (cf.( |A.19| )); instead we de- 
fine ((T'I'PT) )p^p2 having the same expansion as (T"l'PT)pjp2 except for the a„(. . .)Pa]„(. . .) 
are replaced by a„(. . .)Fal^{. . .) "'"^''"' . In this formahsm ((T"l'PT)')p^p2 reads 



((T+PT)') 



PlP2 



dxdy{Mx)My)){^^^T^) 
dtpiix) 5%l)i[y) I 



PlP2 



\ / pip2 

+ / ^(T(*(.)^,fe))) (t'P^^^^j^) . (2.36) 

\ / pip2 



^Remember that P = P (g) P = \qi , q2) {qi , q2\ (81 J2j \j){j\- Here P = J^j \j){j\ behaves as an identity 
in the subspace of heat-bath states. 



CHAPTER 2. HEAT-BATH PARTICLE NUMBER SPECTRUM 30 

with the T being the anti-chronological ordering symbol. Comparing ( |2.36D with ( |A.2CI| ), 
we can interpret the RHS of ( p.36| ) as a sum over all possible distributions of one line 



(corresponding to ipi) inside of each given (T 7^ !) cut diagram constructed out of 
{T^FT) p-^p^. As (|2.36| ) has precisely the same diagrammatical structure as 



{T^Fa^aT)p^p^ - (a^a) (T¥T)j 



P1P2 



dA{N(io)) _ d{N{uj))j 



(cf.( |2.35|) ), it shows that in order to compute^ the numerator of ^^ — ^ 
^ ^'''' one can simply modify the usual (T"'"P T)p-^p^ cut diagrams by the following one-line 
replacements (cf.( p.32|) ). 



(i) For neutral scalar bosons: 






+ {(p{x)a{k;uj)){a\k;iu)ip{y))} 



^'^ 6{e-m'){fB{u;){fB{u:) + l) 



(27r)3 
+ 5+(A;o - a;)(l + /bM) - r (A^o + ^)/bH} e'^'^^^"^), (2.37) 



where fsiyj) is the Bose-Einstein distribution: fsi^) = ^FPuj- Term 6'(— A;o)/b(^) de- 
scribes the absorption of a heat-bath particle, so reduces the number spectrum, that is 
why the negative sign appears in front of it. Analogously, 



Here ^^^^ = / ^6+iu;^ - k' - m2)(at(^,k)a(c.,k)), (cf. (|||)) 
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^31, 



+ {a{k]uj)^{x)){a\k]uj)^{y))} 



"^'^ S{e -m')il + fB{uj))fB{uj)e-''^^~y'> 



{2nf 
X (5+(A;o -uj) + 5-{ko + cu)). (2.38) 



Similarly, for A{N) one needs the following replacements (cf. (|2.33| )) 



MxMy)) - I 0^,6{e-m'){fsM{fBM + i) 



+ 9{ko){l + fsM) - 9i-ko)fBM}e-'''^^-y\ 



in^ixMy))) -. I 0^^6{e-m'){l + MuJ,))fBMe-^'^^-y\ 



with the dispersion relation Uk = vl?— m^ 



(ii) For Dirac fermions: 



(2.39) 



The Dirac field is comprised of two different types of excitations (mutually connected via 
charge conjugation), so the corresponding number operator N{uj) = Nb{uj) + Nd{uj) with 
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Thus, the one-hne replacements needed for dA{N{,{uj))/duj are 



/^3u _ 

- {^lj,{x)b4k;uj)){bi{k;u;)^M)} 



6^{P - m^) 6{ko -uj) {/c + m)^ 



X {(1 - /p(c.)) - ^(c.)(l - /p(c.))}e-'^(--^), 



where fpi^) is the Fermi-Dirac distribution: fpioj) = ^(|^|i^) , and 



(2.40) 



/^3u _ 

— ^ 5+(.;^ - k^ -m^) {(6Jk;c.)V^,(x))(6l(k;c.)^^(y)) 

-(6Uk;cu)^p(x))(6„(k;c.)^,(y))} 



"^ ^ 6+{k^-m'')6{ko-iu){^ + m)^ 






(2.41) 



Correspondingly, for A{Nb) we need 



d rC -1,-2 2\ 



(V'p(x)^^(l/)) -^ I -n^^ {k -m){^ + m)p^ 
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X {(1 - Mco)) - /pM(1 - ^M}e-^'=(^-^) 



(2.42) 

For the d-tjpe excitations the prescription is very similar, actually, in order to get — (^ i 
the following substitutions must be performed in (|2.4CI|) - (|2.42|) : 6{kQ) -^ 6'(— fco), fp — *■ 
(1 — fp) and jU — > — /i. 



(iii) For gauge fields in the axial temporal gauge {A^ = 0): 

The temporal gauge is generally incorporated in the gauge fixing sector of the Lagrangian 
and particularly 

Cf,, = -^iAo)^;a~^0. (2.43) 

The principal advantage of the axial gauges arises from the decoupling the F-P ghosts in the 
theory. This statement is of course trivial in QED as any linear gauge (both for covariant 
and non-covariant case) brings this decoupling automatically |^. Particular advantage 
of the temporal gauge comes from an elimination of non-physical scalar photons from the 
very beginning. 

Let us decompose a gauge field Ai,i = 1, 2, 3 into the transverse and longitudinal part, 
i.e. Ai = AT + AT with 



iT _ I s <^A\ A ^„A aL _ ^i^3 



M = ( 5., - y J A, and AT = ^^A,, (2.44) 
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and use the sum over gauge-particle polarisations 

i:e['\k)ef\k)=5.,-^, (2.45) 

A=l 

with e^'^\k) being polarisation vectors, then for dA{N'^ (to)) / dcu we get the following one- 
line replacements 

2 r ^31, 

+ {Af{x)a,{k;uj)){a{{k;u;)Aj{y))^} 



{^- - ^) (Eq-(2.37)) 



{nAj{x)Aj{y))) -. ^ J -^^ 6^u' -k' -m'){{Aj{x)a{{k;uj)){Aj{y)a,{k;uj)) 



+ {Ai{x)a,{k;uj)){Ai{x)al{k;uj)) 



= (<5., - ^) (Eq.(^)). (2.46) 

The replacements needed for A{N'^) can be concisely expressed as 



(•••)- ('5.. -^)(Eq.(2.39)) (2.47) 

As for the longitudinal (non-physical) degrees of freedom, it is obvious that 



{Afix)Afiy)); {T (A^ (x) Af (y))) ^ Q . (2.48) 
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Eqs.( |2.37| )-( |2.47| ) can be most easily derived in the finite-volume limit, e.g. for a scalar 
field we reformulate ip{x) as 



¥'(^) = Yl 



j\r 



.Ert+zkrX 



y/2ErV 



-It^r 



+ 



Al 



AErt — i]ir'X^ 



y/2ErV 



rescaling the annihilation and creation operators by defining a{k) = \/2EkVAk in such a 
way that [A^; A^,] = 5,,> (so that (4^,,) = ^.../^(fco)), while / ^ ^ 1 ^^. 

The replacements (|2.37|) - (|2.47|) are meant in the following sense: firstly one constructs 
all the T 7^ diagrams for (T'''PT)pjp2, using the LTE ( p.28| ) and the rules mentioned 
therein. In order to calculate the numerator of ( p.32| ) or ( |2.33| ) we simply replace (using 
corresponding prescriptions) one heat-bath particle line in each cut diagram and this re- 
placement must sum for all the possible heat-bath particle lines in the diagram. If more 
types of heat-bath particles are present, we replace only those lines which correspond to 



particles whose number spectrum we want to compute (see FIG J2.6| ). 







Figure 2.6: The numerator of ( 2.32| ) and ( p. 33 ) can be calculated using the modified cut diagrams 
for (T"fPT)p^p2. As an example we depict all the possible contributions to the numerator derived 



from the cut diagram on Fig. 2.4 c). The wavy lines and thin lines describe the heat-bath particles. 
The crossed lines denote the substituted propagators, in this case we wish to calculate the thin-line 
particle number spectrum. 



The terms in the replacements ( p.37| )-( ^.47| ) have a direct physical interpretation. The 
f{uJk) and (1 + {—ly f{ijjk)) can be viewed as the absorption and emission of the heat- 
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bath particles respectively |^. The term f{ujk){l + {—iyf{ujk)) describes the fluctuations 
of the heat bath particles. This is because for the non-interacting heat-bath particles 
{{^k — {nk))"^) = f{ujk){l + {—iyf{ujk))- The substituted propagators can be therefore 
schematically depicted as 

X \/ y ~ fluctuations 

1 ^^ 1 

X \/ y ~ fluctuations + emissions + absorptions 

2 ^^ 1 

X \/ y ~ fluctuations 

Collecting all the contributions from emissions, absorptions and fluctuations separately, one 
can schematically write 

"'V^V^Hf ^ dV^{^>)i _^ pemissionf\ _^ pabsorption f\ _^ _p/^«c/ N (2.49) 

duj dio 

where, for instance for neutral scalar bosons 



\ I piP2 

Using ( p.38|) , it is easy to write down the analogous expressions for the p"-^^°^p^^°'^ and 



pfiuc^ To the lowest perturbative order, the form ( p.49|) was obtained by Landshoff and 



Taylor p|. 



2.4 Model process 
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2.4.1 Basic assumptions 

To illustrate the modified cut diagram technique, we shall restrict ourselves to a toy model, 
namely to a scattering of two neutral scalar particles $ (pions) within a photon heat bath, 
with a pair of scalar charged particles 0, (p ('muon' and 'antimuon') left as a final product. 
Both initial and final particles are supposed to be unheated. We further assume that the 
heat-bath photons A are scalars, i.e. the heat-bath Hamiltonian has form 

H'' = lid.Af - ^A\ 

In order to mimic scalar electrodynamic, we have chosen the interacting Hamiltonian en- 
tering in the T-matrix as 



2.4.2 Calculations 



We can now compute an order-e^ contribution to the — ,^ . The evaluation of the 



dA{N.y{ul)) 



dA(N^(oj)) 
du) 



is straightforward. In FIG.^]^ we list all the modified cut diagrams contributing 



to an order-e^. Note that diagrams b) and c) are topologically identical. Similarly, diagrams 
e), f), h), i) and j) should be taken with combinatorial factor 2 (corresponding diagrams 
with a heat-bath particle line on the bottom solid line are not shown). Of course, diagram 
g) vanishes for kinematic reasons. 



For instance, in order to calculate the contribution from diagram a) (see also FIG. p.8|) 
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Figure 2.7: The modified cut diagrams involved in an order-e^ contribution to the photon number 
spectrum. Dashed hues: photons. Sohd hnes: ip, cfy particles. Bold lines: <1> particles. 




Figure 2.8: The diagram a) with a corresponding kinematics. 



we go back to Eq. (|2.9|) and to prescriptions (|2.37| )- (|2.38| ), so we get 






X iA" (yi - X2) iA~ (1/2 - X2) 



d^k 



5{e-m')il + fBico))Mco) 



X {6+{ko-u;) + 6-{ko + uj))e-'''^y'-y'^ 



CHAPTER 2. HEAT-BATH PARTICLE NUMBER SPECTRUM 39 



—2q^k + m^ 2q4,k + m^ 2^3/;; + m^ —2q4,k + m^ 
X (5(A;2-m2)5^(-pi-p2 + gi + g2). (2.50) 

We have dropped the ie prescription in the propagators since adding/ subtracting an on- 
shell momenta gi;2 to/from an on-shell momenta k we can not fulfil the condition (k ± 
'?i;2)^ = "^u- As it is usual, we have assumed that our interaction is enclosed in a 'time' 
and volume box (t and V respectively). Analogously one can calculate contributions from 



other diagrams in FIG. p.T) . Let us emphasise that it is necessary to give sense to graphs 
e), h), i) and j) as these suffer with the pinch singularity; the muon-particle propagator 
{p\.2 — Tn?)~^ has to be evaluated at its pole because of the presence of an on-shell line 
(1-2 line) with the same momenta. Some regularisation is obviously necessary. Using the 
formal identity |^ 



-^5{x) = -\5\x) T ir^{5{x))\ (2.51) 

we discover that the unwanted S^ mutually cancel between e) and h) diagrams (similarly 
for i) and j) diagrams). An alternative (but lengthier) way of dealing with the latter pinch 
singularity; i.e. switching off the interaction with a heat bath in the remote past and future, 
is discussed in [^ . Evaluating all the diagrams (note, we should attach to each digram the 



factor ^ coming from a Taylor expansion of the T-matrix), we are left with (c.f. Eq.( p.49|) ): 
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;^em^ss^o^n / 



otion I 






{TFT)p,p,V8ujp,ujp,{27r) 



d^k 5{k^ - Tni) 5iko - uj) 



X / d%d%6^{ql-ml)6^{ql-ml: 



X {K,{l + fB{oo))6\~Q + qi + q2 + k) 
-K2fB{uj)6\~Q + qi + q2-k)} 



(2.52) 



and 



pfluct^^^ 



tAV/BM(i + /BH) 

{TFT)p,p,V8iJp,iUp,i2nf 



d^k 5{k^ - ml) 5{ko - to) 



X / d%d%6^{qf~ml)6^{ql-ml) 

X {S\-Q + qi + q2 + k)Ki + 5\-Q + gi + ga - k)K2 

-25\-Q + q^ + q2)K^] 



+ 



tAV/^M(l + /^M) /■^4. .,.2 2w.- ^ 

— — — — — — dk5[k - m 5 fco - ^ 

(T P TT)p^p2 V8u;p,t^p2 (27r)5 J ' 

K j d%d%5\-Q + q,+q2) 

2giA;-m2 2qik + m^^J ^^^ ^'' dml ^^^ ^' 



HQi ^92)} 



(2.53) 



with Ki 



Ko 



Q = Pi + P2- The relevant (i.e. order-e°) term (see FIG. p.9|) 



^ + 252fc-m2 ) ' -^3 - (2gife-m2)(2g2fc+m2) ^'^'^ 
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Figure 2.9: The lowest-order cut diagram for (TPT"f)p^p2. 



for (TPTt)p^p2 reads 



(TPTt)p^p, 



\H 



41 



16 Vujp^ujp^{2'KY 



UVujp,ujp,\Q\{'^T^) 



I d%d% S^ql - ml) 5^{ql - m^ 5\-Q + q, + q,) 
Q2 _ Ami. (2.54) 



Eqs.( |2.5^ ) and ( |2.53| ) are analogous to the result obtained in [^] for the decay. In order 
to understand their structure, let us deal with the number spectruni[] for small a;'s. To do 
this, we change the integration variables 



q2^ q2T ^k. 



These changes lead to 



(2.55) 



{2q,k ± m2) 6+{ql - mj) 6+{q2 - mj) 6\-Q + q, + q^ ± k) 

-^ 2q,k 6+{ql -M^tX) 6+{ql -M^tY) 6\-Q + gi + gs), (2.56) 



where M^ 



m. 



1 2 



;, X = qik and Y = q2k. In addition, transformations (|2.55|) have 



unite Jacobian. If one Taylor expands ( p.56| ) in terms of X and Y then one gets successively 
^So we implicitly assume that the photon mass m^ is sufficiently small. 
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higher cj-contributions to ( p.52| )-( p.53D . Expanding ( p.52| ) to the first order in X and Y, 



and keeping only temperature-dependent pieces, we have 



(TPT"f)p^p2 ^^^Pi^P2 ^pemission^^^ ^ pabsorption ^^^^ 



2^2 



with 



fsioo) f d\d% 6{e - m^) 6{ko - u) A, (2.57) 



^ = ^ / d%d% 5^{ql - Ml) 5^{ql - M^ (AKX) 5\Q ~ q, - q,) \ ^^^,,^^,, . 

Here K = ( ^-j + j^ j (we have performed transformation qi *-> q2 in order to express 
( |2.57| ) solely in terms of X). As A is a Lorentz scalar, it must depend on k only via product 
(kQ). One can thus evaluate A in the frame where Q = {Qo,0) and then replace ujQq by 
{kQ) (see also [^). Straightforward calculations show that 

A = ^=^= ^ '^ ^' 2 • (2-58) 

|Q|^¥-M2 (^(fcg)2 + ^2(^M!_MQ2^^ 

Recalling ( |2.54| ), we get 



:^emission/ \ _i_ rpabsorption / \ 



Q'fBico)e' 



7r2M2 VQo - QWQ' - ^M^^Q' - Amf^ 



M 



a;go + |k||Q|)2 + m2^(f -M^^ 




2^^ mUQ^-M^ 



uQ, - |k||Q|)2 + m2g (f - M2) j ^{cuQo - |k||Q|)2 + m^ (f - M^) 
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^ ' ' (2.59) 



A/2 



12 



(u;Qo+|k||Q|)2 + m2(^f 



M2 



with |k| = y/Lu"^ — m^ and |Q| = a/Qo ~ Q^- Eq.( |2.59| ) takes a particularly simple form if 
m^ is negligibly small (i.e. if m^ <^ cj), then 



;^emissiont, .\ _, rpabsorption / 



[uj] + P '■^■^-^[uj] 



Similarly as in the previous case we can evaluate F^^'^'^*. Performing transformation ( |2.55| ), 
and expanding (|2.53| ) to the first order in X and Y , we get 



(TPTt)p^p,^^:^^^^^i^F^^-*(^) 

6 

\ 2 2 

~ 7^/^(^)(l + fB{^)) I d%d% SiP - ml) 5{k, - cu) B, (2.61) 



with 



B = I d%d% 5\Q -q,- q^) (^^) 5+{ql - ml) 5^{ql - m^ 
+ jd%d%5\Q-q,-q,) H^ 



2 (9^\ ^ ^l ] 5^{ql - Ml) 6+{ql - M 
\qik) dM^dM2] ^^ 1^ ^^2 



2) I 

2 / I Mi=M2=M 



Direct calculations lead to 
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B 



271 {kQf 



M2 



2\ ^2 



-M2)(2M2-Q2)^: 



Q\Q- - 4M^)I I (Mj(fcg)2 + ^2 j^f - M^)) (^Mj^^g)^ + m^ (^ _ m^ 



TT 



271 



|g|v/g2-4M2 |Q|(Q2_4M2) 
After some analysis we finally get 



i?/« 



uct I 



UJ 



fBiuj){l + fBJuj)) m^e^ 
47r2 M2 v/Qg - Q2 Vg2 _ 4^2 



M 

M 



arctg 



^Mo + |k||Q| 



m,V^-M2 



— arctg 



^Mo-lkllQI) 



S!_M2 



771^ 



7V 4 



+ 



(2M^ - g^) m^ 
2v/g2 - 4M2 



co-Qo + |k||Q| 



Af2 



:^Qo + |k||Q|)2 + m2(^-M2) 



tuQo — |k||Q| 



M2 



f^(c.go-|k||Q|)2 + m2(^-M2) 



/sH(l + /sH)|k|e2 



7r2(Q2 



4m2) 



(2.62) 



Expression ( p.62| ) considerably simplifies in the limit m^ — i> 0. In the latter case 



pfi 



uct 



fB{uj){l + fB{uj))uje' 



(2.63) 



7r2(g2 _ 4^2) 

so the leading behaviour for F-/^^"^* at small u and m^ <^ cu is dominated by u;~^. Note that 
separate contributions to the 0-th order of a Taylor expansion of F-^'"'^ behave as uj'"^ but 
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they cancel between themselves leaving behind parts proportional at worst to uj~^. The 
minus sign in ( p.63| ) reflects the fact that the fluctuations tend to suppress an increase in 



the particle number spectrum when cu is small. On the other hand, from (|2.60|) we see that 
the emissions and absorptions stimulate an increase in the particle number spectrum for 
small u. 



A result similar to (|2.6CI|) and ( p.63| ) has been derived by Landshoff and Taylor ^5 



for a decay using proper scalar electrodynamics, though in their case a contribution from 
the emission and absorption dominated over fluctuations for small uj. Note that in our 
model both contributions are of comparable size at tu ~ 0. The former feature is inherently 
connected with the fact that our 'photons' are scalar particles. If photons were vector 
particles an additional photon momentum fc^ would go with each three-line photon-muon 
vertex and so one might expect that the contributions ( |2.6CI| ) and ( |2.63| ) would be 'soften' at 
small UJ. We have checked explicitly that for zero-mass photons in the axial temporal gauge 
(i.e. v4° = 0) this is indeed the case, and it was found that F^'"*""^"" + F'^''"°^p**°'^ oc uj-^ 
whilst F-^'"^* oc UJ. 

Until now we have supposed that our heat bath contains only (scalar) photons in thermal 
equilibrium. However, one could similarly treat a heat bath which is comprised of photons 
and charged particles, let say electrons, mutually coexisting in thermal equilibrium. To be 
more speciflc, let us assume that the heat-bath photons A and electrons \E' are both scalars 
so the heat-bath Hamiltonian takes form 

2 

1 Tfl 

H^ = -^[d^Af~^A\ (2.64) 



and the T-matrix interacting Hamiltonian Hm reads 
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i7,„ = ^$2^0t + (e^ + ^A')^^^ + {eA + ^A' 



It is usually argued |]57| , |58| that the interacting pieces in H'^'' can be dropped provided that 
ti — i> — oo and tf -^ oo. Since we assume that 'pions' are prepared in the remote past and 
'muons' are measured in the remote future, we shall accept in the following this omission. 
We can now approach to calculate both the photon and electron number spectrum, i.e. 
— ,^ and — ^ "^^'f respectively. As for — ,^ , an order-e^ contribution is clearly 
done only by diagrams in Fig. |2.7| as there are no relevant graphs with electron vertices 
contributing to this order, so ( |2.59| ) and ( |2.62|) still remain true. On the other hand, there 
is no order-e^ contribution to — ^j^^ ■ The lowest order in e (keeping A^ fixed) is e^. This 
brings richer diagrammatic structure then in the photon case. In FIG. p.lO| we list all the 
generating thermal diagrams contributing to an order-e^. 






Figure 2.10: The generating thermal diagrams involved in an order-e^ contribution to the electron 
number spectrum. Dashed lines: photons. Thin lines: <j), (fr particles. Bold lines: <1> particles. 
Half-bold lines: electrons. 



It is easy to see that out of these 6 generating thermal diagrams we get 43 non-vanishing 
and topologically inequivalent modified cut diagrams; for example from FIG. |2.1(i| c) we have 



only those diagrams which are depicted in FIG. |2.1i . 
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Figure 2.11: The non-vanishing modified cut diagrams from FIG 



Note, the graphs of FIG. |2.1i] must be multiphed by a factor of four as there are two 
equivalent insertions of the modified electron line and two equivalent distributions of the 
photon-muon vertex (so together with ^ from a Taylor expansion of the T-matrix we get 
the symmetry factor 2). Analogously we get 10 inequivalent modified cut diagrams from 
FIG. |2.1(]| a); 7 from b); 8 from d); 6 from e) and 8 from f).) The actual electron number 
spectrum calculations are thus rather involved. Nevertheless, one might evaluate fairly 
quickly p'^''n'^s'^^°^[ijj^"^ -|- jpabsorpuon^^^^ ^^ there are only three diagrams which contribute, 

namely: 






Let us remind that in the finale state we must have, apart from the heat-bath particles, 
only two 'muons', and so the diagram 




can not contribute to 



duj 



Subtracting a temperature independent part, we are left 
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with 



with 



:^emissioni 



UJp 



tion. 



I rpabsorpti' 



UJ. 



I d% d% S-'iq', - ml) 5{ql - ml) 5{ql - u. 



{TVT^),,,,ViUp,iUp,i27T)^ 
X I d% d% 5+{ql - ml) 5^{ql - mj) [K, 5\-Q + q, + q, + q,- g,) 

-K2 5\-Q + q2 + qz + q7 + qe)], (2.65) 



Ki 



1 



{-q2Q + Q^ + ie) {~q-iQ + Q^ - It) {-2qjq, + 2ml - ml] 



K2 = Ki{qe -^ -qe)- 



If we are interested in the quahtative behaviour of ( p.65|) at small tu's, we need to perform 
an integration over pe only. In order to keep our calculations as simple as possible, let us 
assume that mg = m^ = 0. Eq.( |2.65| ) can now be handled in a similar way as in the photon 
heat bath case. We first perform a transformation 



qe -^ qe- 



So ( p. 651) now reads 



(2.65) 



tX^e' fniuJe) 



(TPTt)p^p,\/a;p^^p,(27r)8 
1 

X 



j d%d%5\ql-ml)5+{ql-ml) 



;-g2Q + Q2 + ze) [-q^Q + Q2-^e) 



S, 



(2.66) 



CHAPTER 2. HEAT-BATH PARTICLE NUMBER SPECTRUM 



49 



where 



B 



7 " ye 



{5+(g? + X) - 5+(g2 - X)} 6{ql) 6{q', - u^) 5\-Q + q2 + qz + qi), 



{'^qiqe? 

with X = 2qiqe. As before we might expand B in terms of X. First surviving term reads 



B 



I d% d% (5,.5(g?)) 6\-Q + q2 + q3 + ^r) 



2X 



i^qrqe 



-uol d^. I d% 6{q', - m') S\-Q + q2 + q3 + q?) ^In f ^|— }^ 

Iq?! V?7 + lq7l 



m=0 



(2.67) 



so B cc cijg, and consequently i^^^^^'^^on^j^^^ _|_ ^absorption ^^^^^ o^ "^eT^. Straightforward apphca- 
tion of the previous mathematical operations to F-^'"^*(ti;e) reveals that F'^'"^*(ci;e) oc u!~^ as 
well. Let us mention that the separate contributions present in pf^'^^-''s^on^^^-^^ jpabsorpuon^^^-^ 
and F'^'"^*(ti;e) behave as uj~'^ but they mutually cancel leaving behind terms proportional 
at worst to uj~^. 



Surprisingly enough, we have found that, for small u, our heat bath ( p.64|) changes 
due to scattering $$ -^ (jxp in such a way that the rate of change in the electron number 
spectrum has qualitatively similar behaviour (i.e. uj~^) as the rate of change in the photon 
number spectrum. This is so provided one assumes that both electrons and photons are 
massless particles. Clearly, cu"^ behaviour would be disastrous as it would suggest that 
the energy density udN/du of the heat-bath particles behaves as W^ which would, if 
integrated, produce an infinite contribution to the total energy carried off by the heat-bath 
particles. 
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2.5 Conclusions 

In this chapter we have formulated a systematic method for studying the heat-bath particle 
number spectrum using modified cut diagrams. In particular, for the quark-gluon plasma 
in thermodynamical equilibrium our approach should be useful as an effective alternative 
to the Landshoff and Taylor |4^ approach. The method used in |^ (i.e. to start from 
first principles) suffers from the lack of a systematic computational approach for higher 
orders in coupling constants. One of the corner stones of our formalism is the largest- 
time equation (LTE). We have shown how the zero-temperature LTE can be extended to 
finite temperature. During the course of this analysis, we have emphasised some important 
aspects of the finite-temperature extension which are worth mentioning. Firstly, many of 
kinematic rules valid for zero-temperature diagrams can not be directly used in the finite- 
temperature ones. This is because the emission or absorption of heat-bath particles make 
it impossible to fix some particular direction to a diagrammatic line. It turns out that 
one finds more diagrams then one used to have at zero temperature. The most important 
reductions of the diagrams have been proved. The rather complicated structure of the 
finite-temperature diagrams brings into play another complication: uncutable diagrams. 
It is well known that at zero temperature one can always make only one cut in each cut 
diagram (this can be viewed as a consequence of the unitarity condition). This is not 
true however at finite temperature. We have found it as useful to start fully with the 
LTE analysis which is in terms of type 1 and type 2 vertices. This language allows us to 
construct systematically all the cut diagrams. We have refrained from an explicit use of the 
cuts in finite-temperature diagrams as those are ambiguous and therefore rather obscure 
the analysis. 

The second, rather technical, corner stone are the (functional) thermal Dyson-Schwinger 
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equations. We have developed a formalism of the arrowed variations acting directly on field 
operators. This provides an elegant technique for dealing in a practical fashion with expec- 
tation values (both thermal and vacuum) whenever functions or functionals of fields admit 
the decomposition (|A.1|) . The merit of the Dyson-Schwinger equations is that they allow 



us to rewrite an expectation value of some functional of field in terms of expectation values 
of less complicated functionals. Some illustrations of this and further thermal functional 



identities are derived in Appendix [A.l 



When we have studied the heat-bath particle number spectrum, we applied the Dyson- 
Schwinger equations both to numerator and denominator of corresponding expression. The 
results were almost the same. The simple modification of one propagator rendered both 
equal. We could reflect this on a diagrammatical level very easily as the denominator was 
fully expressible in terms of thermal cut diagrams. Our final rule for the heat-bath particle 
spectrum is 

dAjNju;)) ^ (TtPT)^^^^ 

duj (TtPT)p,p/ ' 

with T being the T-matrix, P being the projection operator onto final states, pi,p2 being 

the momenta of particles in the initial state, f3 being the inverse temperature and M being 

abbreviation for the modified diagrams. Modification of the cut diagrams consist of the 

substitution in turn of each heat-bath particle line by an altered one. This substitution 

must be done in each cut diagram. Replacement must be only one per modified diagram. 

Our approach is demonstrated on a simple model where two scalar particles ('pions') scatter, 

within a photon heat bath, into a pair of charged particles ('muon' and 'antimuon') and 

we explicitly calculate the resulting changes in the number spectra of the photons and. 

It is also discussed how the results will change if the photon heat bath is replaced with 
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photon-electron one. 



Chapter 3 

Pressure at thermal equilibrium 



3.1 Introduction 

A significant quantity of pliysical interest tfiat one may want to calculate in field theory 
at finite temperature, either at equilibrium or out of equilibrium, is pressure. In thermal 
quantum field theory (both in the real- and imaginary-time formalism) where one usually 
deals with systems in thermal equilibrium there is an easy prescription for a pressure cal- 
culation. The latter is based on the observation that for thermally equilibrated systems the 
grand canonical partition function Z is given as 

Z = e-^^ = Tr(e-''(^-'^'^»)), (3.1) 

where Q is the grand canonical potential, H is the Hamiltonian, Ni are conserved charges, 
fii are corresponding chemical potentials, and (3 is the inverse temperature: (3 = 1/T. Using 
identity /3^ = —T-^ together with ( |3.1| ) one gets 

t(^] =n-E + ix,N,, (3.2) 

with E and V being the averaged energy and volume of the system respectively. A com- 

53 
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parison of ( |3.2| ) with a corresponding thermodynamic expression for the grand canonical 
potential [0,0,|^ requires that entropy S = — (|p) y, so that 

dn = -SdT - pdV - Nidfx, ^ p=-(—j . (3.3) 

For large systems one can usually neglect surface effects so E and A^j become extensive 
quantities. Eq. (|3^ ) then immediately implies that Q is extensive quantity as well, so (|3.3| ) 

simplifies to 

n \nZ , ,, 

p= = . 3.4 

The pressure defined by Eq.( p.4|) is the so called thermodynamic pressure. 

Since InZ can be systematically calculated summing up all connected closed diagrams 
(i.e. bubble diagrams) 0, ^, 0, the pressure calculated via ( |3.4| ) enjoys a considerable 
popularity |3^, |36|, ^, ^ . Unfortunately, the latter procedure can not be extended to out 



of equilibrium as there is, in general, no definition of the partition function Z nor grand 
canonical potential Q away from an equilibrium. 

Yet another, alternative definition of the pressure not hinging on thermodynamics can 
be provided; namely the hydrostatic pressure which is formulated through the energy- 
momentum tensor Q'^'^. The formal argument leading to the hydrostatic pressure in D 
space-time dimensions is based on the observation that (0°-'(x)) is the mean (or macro- 
scopic) density of momenta p-' at the point x^. Let P be the mean total (D — l)-momentum 
of an infinitesimal volume V'^^^^^ centred at x, then the rate of change of j-component of 
P reads 



dPHx) f .n_. , d "" ' 



Jl/(D-1) OX" ^ JdY^D-^) 



dt 

1=1 
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In the second equality we have exploited the continuity equation for (9'^-^) and successively 
we have used Gauss's theoremQ The dV^^^^'' corresponds to the surface of V^^^^\ 

Anticipating a system out of equilibrium, we must assume a non-trivial distribution 
of the mean particle four- velocity U^{x) (hydrodynamic velocity). Now, a pressure is by 
definition a scalar quantity. This particularly means that it should not depend on the 
hydrodynamic velocity. We must thus go to the local rest frame and evaluate pressure 
there. However, in the local rest frame, unlike the equilibrium, the notion of a pressure 
acting equally in all directions is lost. In order to retain the scalar character of pressure, one 
customarily defines the pressure at a point (in the following denoted as p{x)) [Q, which is 
simply the 'averaged pressure' Qover all directions at a given point. In the local rest frame 
Eq.( p.5| ) describes j-component of the force exerted by the medium on the infinitesimal 
volume V^^^^\ (By definition, there is no contribution to dP^{x)/dt caused by the particle 
convection through dV^^^^\) Averaging the LHS of ( |3.5| ) over all directions of the normal 
n(x), we getfl 



D-1 

= TTT^^Y. / ds{Q\ix')), (3.6) 

where dil{n) is an element of solid angle about n and 5*^ ~^ is the surface of {D — 2)- 

sphere with unit radius {J dil{n) = S^^^ = 27r^~/r(-^^)) . On the other hand, from the 

-'^The macroscopic conservation law for (6^'') (i.e. the continuity equation) has to be postulated. For 
some systems, however, the later can be directly derived from the corresponding microscopic conservation 
law H. 

^To be precise, we should talk about averaging the normal components of stress p^ . 

•^The angular average is standardly defined for scalars (say, A) as; / A dfl{n)/ J dn{n), and for vectors 
(say, A*) as; J2j I ■^"' ^'' dn{n)/ J dil{n). Similarly we might write down the angular averages for tensors 
of a higher rank. 
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definition of tlie pressure at a point x'^ we miglit write 



D-l 



iSF^T' E / ^^ ^'' ^^(") = -P(^) / ds, (3.7) 

fiere tlie minus sign reflects that the force responsible for a compression (conventionally 
assigned as a positive pressure) has reversed orientation than the surface normals n (point- 
ing outward). In order to keep track with the standard text-book definition of a sign of 
a pressure |l37|,|6^ we have used in (|3.7|) the normal n in a contravariant notation (note, 
n* = —rii). Comparing ( p.6|) with (|3.7|) we can write for a sufficiently small volume V^^~^^ 

p{x) = -j^^^^{eMx)). (3.8) 

We should point out that in equilibrium the thermodynamic pressure is usually identified 
with the hydrostatic one via the virial theorem I^B],^. In the remainder of this chapter we 
shall deal with the hydrostatic pressure at equilibrium. We shall denote the foregoing as 
V{T), where T stands for temperature. We consider the non-equilibrium case in the next 
chapter. 

The plan of this chapter is as follows. In Section \i.'2\ we review the necessary math- 
ematical framework needed for the renormalisation of the energy-momentum tensor. (For 
an extensive review see also refs. [^,|6^,^.) The latter is discussed on the 0{N) $^ the- 



ory. As a byproduct we renormalise $^, ^a^b and O^*^ operators. The corresponding QFT 
extension of (|3.^ ) is obtained. 

Resumed form for the pressure in the large- A^ limit, together with the discussion of both 



coupling constant and mass renormalisation is worked out in Section p.3| . The discussion is 
substantially simplified by means of the Dyson-Schwinger equations. 
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In Section p.4| we end up with the high-temperature expansion of the pressure. Cal- 
culations are performed ioi D = 4 (both for massive and massless fields) and the result is 
expressed in terms of renormalised masses m^,(0) and mr{T). The former is done by means 
of the Mellin transform technique. 

This chapter is furnished with two appendices. In Appendix we clarify some math- 
ematical manipulations needed in Section |3.3| . For the completeness' sake we compute in 
Appendix |C| the high-temperature expansion of the thermal-mass shift 5m? {T) which will 
prove useful in Section Ol 



3.2 Renormalisation 



If we proceed with ( |3.8| ) to QFT this leads to the notorious difficulties connected with the 
fact that Q^'^ is a (local) composite operator. If only a free theory would be in question then 
the normal ordering prescription would be sufficient to render (G^'^) finite. In the general 
case, when the interacting theory is of interest, one must work with the Zimmerman 'normal' 
ordering prescription instead. Let us demonstrate the latter on the 0{N) $^ theory. (In 
this section we keep A^ arbitrary.) Such a theory is defined by the bare Lagrange function 



a=l \a=l / 



(3.9) 



we assume that m-Q > 0. The corresponding canonical energy-momentum tensor is given 
by 



ef = ^ a'^$,9"<i>, - ^^'^z: . (3.10) 

a 

The Feynman rules for Green's functions with the energy-momentum insertion can be easily 
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explained in momentum space. In the reasonings to follow we shall need the (thermal) 
composite Green's function^ 

D'^'^ix^ly) = {T* {$,(xi) . . . $.(a;„)er (y)}) . (3.11) 

Here the subscript r denotes the renormalised fields in the Heisenberg picture (the internal 
indices are suppressed) and T* is so called T* product (or covariant T product) |]53| , |67H6^1 . 
The T* product is defined in such a way that it is simply the T product with all differential 
operators "D^- pulled out of the T-ordering symbol, i.e. 

T*{V;mxr) . . . V;:^xn)} = V{zd{,y)r{<l>{x,) . . . ^x^)} , (3.12) 

where V{id^^}) is just a useful short-hand notation for V^^^V^^ . . .V^^- ^^ the case of 
thermal Green's functions, the T might be as well a contour-ordering symbol. It is the 
mean value of the T* ordered fields rather than the T ones, which corresponds at T = and 



at equilibrium to the Feynman path integral representation of Green's functions p^,\i 
A typical contribution to Q'^'^{y) can be written as 



so the typical term in (|3.11|) is 



2)(z%i) {T* {<^r{x,) . . . <^rMHyi) . . . Hyk)}) \y^=y • 

Performing the Fourier transform in (|3.11|) we get 

''By $ we shall mean the field in the Heisenberg picture. The subscript H will be introduced in cases 
when a possible ambiguity could occur. 
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D^'^{p-\p)= E /(n 



(27r)^ 



(2vr)^ d^ip-Y.q,) l^[r)(gM) D{p-\q% (3.14) 



where T^'^^^i- . .) is a Fourier transformed differential operator corresponding to the quadratic 
(k=2) and quartic (k=4) terms in Q'^" . Denoting the new vertex corresponding to ^^^^(• • •) 
as ®, we can graphically represent ( ^.111 ) through ( ^.141 ) as 




Bi'''{p^\p) 



P, 



Pn Pi 



Figure 3.1: The graphical representation of D^'^[p'" 



For the case at hand one can easily read off from ( ^.l(j| ) an explicit form of the bare composite 
vertices, the foregoing are 
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(For the internal indices we do not adopt Einstein's summation convention.) We have 



tacitly assumed in FIG. |3.1| that the vacuum bubble diagrams present in the shaded blobs 
are divided out. We have also implicitly assumed that summation over internal indices 
is understood. Note that in the case of thermal composite Green's function, the new 
vertices are clearly of type-1 as the fields from which they are deduced have all a real-time 
argument^ (type-1 fields). 

Renormalisation o/$a(x)$(,(x) 



Now, if there would be no G^*^ insertion in ( p.ll| ), the latter would be finite, and so it is 



natural to define the renormalised energy-momentum tensor [O^'^] (or Zimmermann normal 
ordering) in such a way that 

Dr(x"|y) = (T* {^r{xi) . . . ^.(x^) [er]}), 

is finite for any n > 0. To see what is involved, we illustrate the mechanism of the 
composite operator renormalisation on $a(x)$f,(x). We shall use the mass-independent 
renormalisation (or minimal subtraction scheme - (MS)) which is particularly suitable for 
this purpose. In MS we can expand the bare parameters into the Laurent series which has 
a simple form []53| , [66| ,|70|, namely 



A„ = /-A,(l + |:^^j (3.15) 

Here Oq and h^ are analytic in D = 4. The parameter /x is the scale introduced by the 
^For a brief introduction to the real-time formalism in thermal QFT see for example pSj]. 
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renorinalisation in order to keep A^ dimensionless. An important point is that both a^s 
and 6fc's are mass, temperature and momentum independent. 

It was Zimmermann who first reahzed that the forest formula known from the ordinary 
Green's function renormahsation |5^, ^ can be also utilised for the composite Green's 



functions rendering them finite |6^,|71[ . That is, we start with Feynman diagrams expressed 
in terms of physical (i.e. finite) coupling constants and masses. As we calculate diagrams 
to a given order, we meet UV divergences which might be cancelled by adding counterterm 
diagrams. The forest formula then prescribes how to systematically cancel all the UV loop 
divergences by counterterms to all orders. However, in contrast to the coupling constant 
renormahsation, the composite vertex need not to be renormalised multiplicatively. We 
shall illustrate this fact in the sequel. Let us also observe that in the lowest order (no loop) 
the renormalised composite vertex equals to the bare one, and so to that order A = [A], 
for any composite operator A. 

Now, from ( ^.15| ) and ( |3.16D follows that for any function F = F{mr, Xr) we have 



OF dml OF ml OF 



dm^ dm^ dml mj dml 
So particularly for 

F = D{xi, ...,Xn) = (T*{<l>,(a;i) . . . $.(x„)}), 

one reads 

d d 

""' d^^^""^' . . . , x„) = m^ Q^D{x,, . . . , x„) 

N 

^ 1 A/ / d^x 

a=l 



A/" / d^x ^ l^4> 0r(xi) . ..(prixn) ml(j)l{x) exp{iS[(j),T]) 

^ a=l ^ 

( ~2 ) / ^^^ 5Z ^"(^1' ■■■,Xn\x; ml) . (3.17) 

^ ^ -^ a=l 
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Here A^"^ is the standard denominator of the path integral representation of Green's func- 
tion. We should apply the derivative also on Af but this would produce disconnected graphs 
with bubble diagrams. The former precisely cancel the very same disconnected graphs in 
the first term, so we are finally left with no bubble diagrams in ( p. 171) . In the Fourier space 
([3T7| ) reads 






^(pi,---,Pn)= (-^)i^ Daipi,...,Pn\0;ml). (3.18) 



As the LHS is finite, there cannot be any pole terms on the RHS either, and so ^^ ^^o'^a 
is by itself a renormalised composite operator. We see that mg precisely compensates the 
singularity of ^^^^ $2. 

Now, it is well known that any second-rank tensor (say Mab) can be generally decom- 
posed into three irreducible tensors; an antisymmetric tensor, a symmetric traceless tensor 
and an invariant tensor. Let us set Mat = ^a^t, so the symmetric traceless tensor Kab 
reads 



K,b{x) = $„(a;)$b(x) - 6ab/N J^ *c(^) , (3-19) 

c=l 

whilst the invariant tensor lab is 

N 

2/ 



Iab{x)=Sab/NY,'^ 



c=l 



Because the renormalised composite operators have to preserve a tensorial structure of the 
bare ones, we immediately have that 

Kab = A,[Kab] and 4^ = ^[4^] , (3.20) 
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where both Ai and A2 must have structure (1 + ^ (poles)). The foregoing guarantees that 
to the lowest order Kab = [Kab] and lab = [lab]- As we saw in ( |3.18| ), mllab is renormalised, 
and so from (|3.2(]| ) follows that m^Iab = C [lab]- Here C has dimension [m^] and is analytic 
in D = 4. We can uniquely set C = m^. because only this choice fulfils the lowest order 
condition lab = [lat] (cf. Eq.( |3.16|) ). Collecting our results together we might write 



i:*? 



;S$2 



E^ 



^S$2 J][$ 



cj 5 



(3.21) 



with Zj^^2 = A2 = ^^. In the second equality we have used an obvious linearity p5[ of 

rriQ 

[...]. From ( p39D and (lOTI ) follows that 



Sa 

N 



N 



M^)Mx) = A,[<l>a{x)<l>b{^)] - "-^{A, - Z^^.) Y}^l{x)] 
So particularly for $^ one reads 



(3.22) 



c=l 



$ 






(3.23) 



From the discussion above it does not seem to be possible to obtain more information about 
Ai without doing an explicit perturbative calculations, however, it is easy to demonstrate 
that Ai 7^ 2'2#2. To show this, let us consider the simplest non-trivial case; i.e. N=2, 
and calculate Ai to order A^. For that we need to discuss the renormalisation of the n- 
point composite Green's function with, say, (^\ insertion. To do that, it suffices to discuss 
the renormalisation of the corresponding IPI n-point Green's function. The perturbative 
expansion for the composite vertex to order A^ can be easily generated via the Dyson- 
Schwinger equation |72| and it reads 
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(3.24) 



where 
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Here cross-hatched blobs refer to (renormahsed) IPI (n + 2)-point Green's function, circled 
indices mark a type of the field propagated on the indicated line, and uncircled numbers 
refer to thermal indices (we explicitly indicate only relevant thermal indices). The coun- 
terterms, symbolised by a heavy dot, are extracted from the boxed diagrams (elementary 
Zimmermann forests). In MS scheme one gets the following results: 



(D 



I (hw {Dii(g)Dii(-g) - Di2(g)Di2(-g)} |ms pol 



e term 



J {2pi 

49^2 ( ^(^ A, /x^-^ mf-2 ) Ims = -A,/i4-^/2 {D - 4) (Any 



■'mi 



D 
(47r)T 



3^«>^ 



-A./i^-^/e {D - 4) (An) 



Here Dn and D12 are the usual thermal propagators in the real-time formalism p^[36| ,p8 
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(see also Section p.3| ). From (4.24) we can directly read off that 



[*?] = ( 1 - o.n^^'^]w.„^. + ^(A.^)l ^l + f-TTTT^^^ + ^(A?) ) ^l 



2(D-4)(47r)2 ^^-'^^y ^ \^ 6(D-4)(47r)2 

As the coefficient before $2 is not zero, we conclude that Ai 7^ Zj^^2. It is not a great 
challenge to repeat the previous calculations for the $i$2 insertion. The latter gives 



Eg. ( p. 231 ) exhibits the so called operator mixing [Q; the renormalisation of $2 cannot be 
considered independently of the renormalisation of $^ (c 7^ a) . The latter is a general fea- 
ture of composite operator renormalisation. Note, however, that $a*^fe (ct 7^ b) do not mix 
by renormalisation, i.e. they renormalise multiplicatively. It can be shown that compos- 
ite operators mix under renormalisation only with those composite operators which have 
dimension less or equal |53|,|65|,[7T[|. 

Unfortunately, if we apply the previous arguments to ra = 0, the result is not finite; 
another additional renormalisation must be performed. The fact that the expectation 
values of [. . .] are generally UV divergent, in spite of being finite for the composite Green's 
function^, can be nicely illustrated with the composite operator [$2] in the A^ = 1 theory. 
Taking the diagrams for D{0\0) and applying successively the (unrenormalised) Dyson- 



Schwinger equation |72] we get 



^Also called the matrix elements of 
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+ 




(3.25) 



Eq.(3.25) might be rewritten as 



D(0|0) 



^(0|0)|ao 



1 

'36 



6 



d''q^ 



n(i^^''(E^^)^""'(^'|o)iA?^(^')^ 



i=l 



(3.26) 



where D'^'^P{q"^\0)\xk is the m-point amputated composite Green's function to order A^, 
and D{q'^) is the full ?7i-point Green's function. The crucial point is that we can write 
il'(OlO) as a sum of terms, which, apart from the first (free field) diagram, are factorised to 
the product of the composite Green's function with n > and the full Green's function. 
(The factorisation is represented in (2.17) by the dashed lines. ) 

Now, utilising the counterterm renormalisation to the last two diagrams in (3.25) we 



get situation depicted in FIG. |3.2| . Terms inside of the parentheses are finite, this is be- 
cause both the composite Green's functions (n > 2 !) and the full Green's functions are 
finite after renormalisation. The counterterm diagrams, which appear on the RHS of the 
parentheses, precisely cancel the UV divergences coming from the loop integrations over 
momenta qi . . .qi which must be finally performed. The heavy dots schematically indicates 
the corresponding counterterms. In the spirit of the counterterm renormalisation we should 
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finally subtract the counterterin associated with the overall superficial divergence^ related 
to the diagrams in question. But as we saw this is not necessary; individual counterterm 
diagrams (Zimmermann forests) mutually cancel their divergences leaving behind a finite 
result. 
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Figure 3.2: Counterterm renormalisation of the last two diagrams in Eq.(3.25). (Cut legs indicate 
amputations.) 



So the only UV divergence in Eq.(3.25) which cannot be cured by existing counterterms 
is that coming from the first (i.e. free field or ring) diagram. The foregoing divergence 
is evidently temperature independent (to see that, simply use an explicit form of the free 

thermal propagator ©n). Hence, if we define 

^A simple power counting in the $^ theory reveals | |53| that for a composite operator A with dimension 
LOA the superficial degree of divergence lo corresponding to an n-point diagram is cj = lda ~ "■• 
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($')renorm = ([$1) " (0|[$1|0), (3.27) 



or, alternatively 



($')renorm = ([$']) " ([$'])|free fields , (3-28) 

we get finite quantities, as desired. On the other hand, we should emphasise that 

($2) - (0|$'|0) = Z^2 {([$2]) - (0|[$2]|0)} ^ finite in D=A. (3.29) 

An extension of the previous reasonings to any A^ > 1 is straightforward, only difference is 
that we must deal with operator mixing which makes ( |3.27|) and ( p.28| ) less trivial. 



The important lesson which we have learnt here is that the naive "double dotted" normal 
product (i.e. subtraction of the vacuum expectation value from a given operator) does not 
generally give a finite result. The former is perfectly suited for the free theory {Zy,^2 = 1) 
but in the interacting case we must resort to the prescription ( p.27|) or ( p.28|) instead. 



Renormalisation of the energy-momentum tensor 

In order to calculate the hydrostatic pressure, we need to find such {Q'^^) Irenorm which apart 
from being finite is also consistent with our derivation of the hydrostatic pressure performed 
in the introductory section. In view of the previous treatment, we cannot, however, expect 
that Q'^^ will be renormalised multiplicatively. Instead new terms with a different structure 
than 0^^ itself will be generated during renormalisation. The latter must add up to G^'^ in 

order to render D^'^{x^\y) finite^. 

^In fact it can be sliown ||36|,|65[| tliat tlie Noether currents corresponding to a given internal symmetry 
are renormalised, i.e J° = [J°], however, this is not the case for the Noether currents corresponding to 
external symmetries (like G^'' is). 
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Now, the key ingredient exploited in Eq.( p.5| ) is the conservation law (continuity equa- 
tion). It is well known that one can 'modify' G^f*^ in such a way that the new tensor B^*^ 
preserves the convergence properties of Q^'^ . Such a modification (the Pauli transformation) 
reads 

For scalar fields, ( |3.3CI| ) is the only transformation which neither changes the divergence 



properties of Q^ nor the generators of the Poincare group constructed out of Q^ |3^,5i 
|55| . Because the renormalised (or improved) energy momentum tensor must be conserved 
(otherwise theory would be anomalous), it has to mix with Q^ under renormalisation only 
via the Pauli transformation, i.e. 

[Of] = ef + dxX^^^ . (3.31) 

In order to determine X^^^ ^ we should realize that its role is to cancel divergences present in 
Q^ . Such a cancellation can be, however, performed only by means of composite operators 
which are even in the number of fields (note that Q^ is even in fields and Green's functions 
with the odd number of fields vanish). Recalling the condition that renormalisation can 
mix only operators with dimension less or equal, we see that the dimension of X^^^ must 
be D — 1, and that X^'^" must be quadratic in fields. The only possible form which is 
compatible with tensorial structure ( p.30| ) is then 

N 

X"'" = J2 ('-biXr-, D) {d^g"" - d^g'^") $„<l>6 . (3.32) 

a,b=l 
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Prom the fact that both B^'^ and [Q^"] are 0{N) invariant (see Eq.( plOD ), ^xX^^"' must 
be also 0{N) invariant, so Cab = SabC. Thus, finally we can write 



N 



[0 



flU] 



ef + c(A,; D) J2 {^'"d'' - g^'d^) <I>1 , 



(3.33) 



a=l 



with c = Co + ^(poles), here cq is analytic in D. Structure of c{\r] D) could be further 
determined, similarly as in the A^ = 1 theory, employing a renormalisation group equation 



| 66| . We do not intend to do that as the detailed structure of c will show totally irrelevant 
for the following discussion, however, it turns out to be important in non-equilibrium case. 
Now, similarly as before, [O^'^] gives the finite composite Green's functions if n > but 
the expectation value ([Oc'^]) is divergent (discussion for the A^ = 1 theory can be found in 
Brown p^). The unrenormalised Dyson-Schwinger equation for D^'^{0\0) reads 



p = 




p = 



( 



( 



36 




) 



(3.34) 



The structure of the composite vertices in (3.34) is that described at the beginning of this 
section. Note that the amputated composite Green's functions in individual parentheses 
are of the same order in A^. Performing the counterterm renormalisation as in the case of 
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([$^]), we factorise the graphs inside of parentheses into the product of the renormahsed 
2- (and 6-) point composite Green's function and the renormahsed fuU 2- (and 6-) point 
Green's function. The latter are finite. The UV divergences arisen during the integrations 
over momenta connecting both composite and full Green's functions are precisely cancelled 
by the remaining counterterm diagrams. The only divergence comes from the free-field 
contribution, more precisely from the T = ring diagram. Defining 

(er)lrenorm = ([6^) " (0|[er]|0), (3.35) 

or 

(er)lrenorm = ([6^) " ([er])lfree field, (3-36) 

we get the finite expressions. Note that the conservation law is manifest in both cases. 
In equilibrium (and in T = 0) we can, due to space-time translational invariance of (...), 
write 

([©r]) = (©D + d,{X'n = (©D • (3.37) 



Using ( [3.351) or (|3.36 ) we get either the thermal interaction pressure or the interaction 



pressure, respectively. This can be explicitly written as 



D-l 



nh.int.(T) = ViT) - no) = - p^ E {(®-) - (0|©-|0)} ' (3-38) 

or 

1 ^"' 
VintXT) = V{T) - Pfree fieldl^) = "^3^1 ^ ^^®'^^^ " (©c.)lfree field} • (3-39) 



In order to keep connection with calculations done by Drummond et al. in [^ we shall in 



the sequel deal with the thermal interaction pressure only. If instead of an equilibrium, a 
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non-equilibrium medium would be in question, translational invariance of (. . .) might be 
lost, in that case either prescription ( p.35| ) or ( p.36| ) is obligatory, and consequently c{\r; D) 



in (|3.33|) must be further specified. 



3.3 Hydrostatic pressure 

In the previous Section we have prepared ground for a hydrostatic pressure calculations in 
the 0{N) $^ theory. In this section we aim to apply the previous results to the massive 
0{N) <l>^ theory in the large-iV limit. Anticipating an out of equilibrium application, 
we shall use the real-time formalism even if the imaginary-time one is more natural in 
the equilibrium context. As we aim to evaluate the hydrostatic pressure in 4 dimensions, 
we use here, similarly as in the previous section, the usual dimensional regularisation to 
regulate the theory (i.e. here and throughout we keep D slightly away from the physical 
value D = A). 

Let us start first with some essentials of our model system at finite temperature. 

3.3.1 Mass renormalisation 

In the Dyson multiplicative renormalisation the fact that the complete propagator has a 



pole at the physical mass leads to the usual mass renormalisation prescription ||53[| : 



ml = ml + S(m^) , (3.40) 

where rrir is renormalised mass and S(m^) is the proper self-energy evaluated at the mass 
shell; p'^ = m^. In fact, Eq. (|3.4CI|) is nothing but the statement that 2-point vertex function 
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(2) 

Tr evaluated at the mass-shell must vanish. The Dyson-Schwinger equations correspond- 
ing to the proper self-energies read |72|-[75|: 



E" 



2 Z^b=l 





E 



N 
b=l 




(3.41) 



S'"^|a^c = ; 



K/N 



where hatched blobs represent 2-point connected Green's functions whilst cross-hatched 
blobs represent proper vertices Tr (i.e. IPI 4-point Green's functions). As S"° are the 
same for all a, we shall simplify notation and write S instead. In the sequel the following 
convention is accepted: 



bd 



^(n) abed.. 



Note that the second term in (3.41) does not contribute in the large- A^ limit. It is easy to see 
that the third term does not contribute either. This is because each hatched blob behaves 

at most as N^ whilst F*^^^ goes maximally^ as A^^^ . Consequently, various contributions 

^ In the <&'' theory there is a simple relation between the number of loops (i), vertices {V) and external 
lines (E); AV = 21 + E. Together with the Euler relation for connected graphs; L = I — V + 1 (here / is 
the number of internal lines), we have L — V = ^^-- As each loop carries maximally a factor of A'' (this 
is saturated only for 'tadpole' loops) and each vertex carries a factor of N^^, the overall blob contribution 



behaves at most as N 



L-V 



N^ 
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from the first graph in (3.41) contribute at most N^, whereas in the second graph the 
contributions contribute up to order A^~^. So the first diagram dominates, provided we 
retain only such 2-point connected Green's functions which are proportional to N^ (as 
mentioned in the footnote, these are comprised only of 'tadpole' loops.). After neglecting 
the 'setting sun' graph. Eg. (3.41) generates upon iterating the so called superdaisy diagrams 
61],^,|7|. 



Let us now define E(m^) = Aq M.{ml). Because the 'tadpole' diagram in (3.2) can be 
easily resumed we observe that 

■M[m,) --J j^ q2_ml- S(m2) + le ~2j (2^ q^ - m^ + ie ' ^ ' ^ 

hence we see that S is external-momentum independent. So if we had started with the 
renormalisation prescription: iVr {p^ = 0) = — m^, we would arrived at ( p.40|) as well (this 
is not the case for A^ = 1!). 

At finite temperature the strategy is analogous. Due to a doubling of degrees of freedom, 
the full propagator is a 2 x 2 matrix. The latter satisfies, similarly as at T = 0, Dyson's 
equation 

D = Di7 + Dj7 (-iS) D . (3.43) 

An important point is that there exists a real, non-singular matrix M (Bogoliubov matrix) 



3J,^,^ having a property that 



^f = m['%- _f^. )m. (3.44) 

Here Ap is the standard Feynman propagator and '*' denotes the complex conjugation. 
Consequently, the full matrix propagator may be written as 
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p2— jtiq — Sy+ie 





(3.45) 



p2— mg — Ey— ie 



Similarly as in many body systems, the position of the (real) pole of D in p^ fixes the 
temperature-dependent effective mass mr{T) [Q, ^. The latter is determined by the 
equation 



m: 



\{T)=ml + Re{T.T{ml{T))) . 



(3.46) 



From the explicit form of M it is possible to show |^,^ that ReSn = ReSr. As before, 
the structure of the proper self-energy can be deduced from the corresponding Dyson- 
Schwinger equation. Following the usual real-time formalism convention (type-1 vertex 
~ — iAo, type-2 vertex ~ zAq ), the former reads: 



-m 



11 



where 




-m 



22 




(3.47) 



1 1/-N1 1 1 2/-^- 2 1 
*. f ^ •. ^ 



and similarly for D22. In (4.49) we have omitted diagrams which are of order 0{1/N) or less. 
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Note that the fact that no 'setting sun' diagrams are present imphes that the off-diagonal 
elements of S are zero. Inspection of Eq.(3.47) reveals that 



^11 — -^r 



I (Cp ^"(^' ^) ""^ ^'^ = - Y / (0^ ^^^(^' ^^ ■ ^^-^ 



It directly follows from Eq. 



that both Sii and S22 are external-momentum inde- 



pendent and reaQ If we define T,T{rnf.(T)) = Aq Airi^-iT))^ then Eq.( |3.46| ) through 
Eq.( ^:iHD imphes that 



ml{T) =ml + Xo MTiml{T)) . 



(3.49) 



A resumed version of Dn is easily obtainable from (|3.45|) |3^ , ^ and consequently ( 3.48| ) 
yields 



-MtK2(T)) 



1 

2 






q^ — mliT) + ie 



+ (47r),5+(g2-m2(T)) 



eQoP - 1 



Im 



(27r)^ eiof^ - 1 g2 _ ^2(7-) + ^g 



(3.50) 



Let us remark that ( p.50| ) is manifestly independent of any particular real-time formalism 
version. This is because the various real-time formalisms [Q, ^ differ only in the off- 
diagonal elements of D. 

In passing it may be mentioned that because Sii(m^) is momentum independent, the 

wave function renormalisation Z$ = 1. (The Kallen-Lehmann representation requires the 

^"Reality of Sn can be perhaps most easily seen from the largest-time equation. The LTE states that 
Sii + S22 + S12 + S21 = 0. Because no 'setting sun' graphs are present, S12 + S21 = 0, on the other hand 



Sii + S22 = 2iImSii (see (p4|) ). 
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renorinalised propagator to have a pole of residue i at p^ = m'^. The former in turn imphes 
that Z$ = (1 — S'^j^(p^)|p2=^2)^^ = 1.) Trivial consequence of the foregoing fact is that 

rF = r(2) and r^') = r^. 

3.3.2 Coupling constant renormalisation 

Let us choose the coupling constant to be defined at T = 0. This will have the advan- 



tage that the high temperature expansion of the pressure (see Section p.4|) will become 
more transparent. In addition, such a choice allows us to stay on a safe ground as the 
renormalisation of the coupling constant at finite temperature is rather delicate [^ . 

By assumption the fields $a have non-vanishing masses, so we can safely choose the 
renormalisation prescription for A^ at s = (s is the standard Mandelstam variable). For 
example, one may require that for the scattering aa — > bb 

r{4) aabb(^^ _ Q) _ -\r/N, {b ^ a) . (3.51) 

The formula ( p.51| ) clearly agrees with the tree level value T\J^^'^ (s = 0) = —\q/N. Let us 
also mention that Ward's identities corresponding to the internal 0{N) symmetry enforce 

p(4) aaaa ^^ ^^^^ ^J^g COUStraiutQ 

+ rW''«"''(pi;p2;P3;P4), (3.52) 



"Actually, Ward's identities read H^ J d^x ^gy M^) = J d^'x ^^ M^) (here (f>a = ^■, W 
is the generating functional of connected Green's functions). Performing successive variations with respect 
to (j)a{v),(t)a{z),(f>a{y) and (f>b{'w), taking the Fourier transform, and setting the physical condition <j)c = 0, 
we get directly (3.52|). 
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for any b ^ a. The structure of T^^^ is encoded in the Dyson-Schwinger equation which 
reads ra,|71: 



a b 






fl._ b 




oT-"' 








,1.6 



N 






^ 



2 ^-'c:i \ 




M ~ "^JN 



(3.53) 



The sum X]i=i schematically represents a summation over various scattering channels. 
Similarly as before, we can argue that the last three graphs contribute at most A^~^, whilst 
the second ('fish') graph may contribute up to order A^^^. So in the large- A^ limit the last 
three diagrams may be neglected, provided we keep in the 4-point vertex function only 
graphs proportional to A^~^. However, the former can be only fulfilled if we retain such a 
'fish' graph where summation over internal index on the loop is allowed. Remaining 'fish' 
graphs (describing t and u scattering channel interactions) are suppressed by the factor 
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A^ ^ as the internal index on the loop is fixed. In this way we are left with the relation 



- _^ _ !^ 

~ N 2N 



J2 I S^T^'^'^'^s) 



(2vr; 



Aq AoA,(iV - 1) [^ 

' N 2A^2 



(g2 — ml + it) {{q — QY — m^ + ie) 



dx 



{2tt)^ (g^ - m^ + a;(l - a;)s + ie) 



• \2 



s=0 



(3.54) 



with Q = Pi + P2 and s = Q^, Pi,P2 are the external momenta. To leading order in 1/A^ 
we may equivalently write 



Ar = Ao + AoA, M'{ml) , (3.55) 

the prime means differentiation with respect to m'^; A^(??i^) is defined by (|3.42|) . Evaluating 
explicitly A^'(m^), we get from (p.55|) 



An 



Af- 



l-A,r(2-f) (m,)^-V2 (47r)- 



(3.56) 



Assuming that both Aq > and A,, > 0, we can infer from ( p.56| ) that 



> A, > 



2(47r)T(m,. 

r(2-f) 



,4-D 



(3.57) 



so for D = 4 we inevitably get that A^ = 0. The latter indicates that the theory is 
trivial |^,^0, or, in other words, the 0{N) $'^ theory is a renormalised free theory in 
the large-iV limit. This conclusion is also consistent with the observation that the theory 
does not posses any non-trivial UV fixed point in the large-iV limit [0,|T^ . 
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On the other hand, if we were assuming that Aq < 0, we would get a non-trivial 
renormalised field theory in D = 4 (actually, from ( p.56| ) we see that Aq — > 0_ , provided that 
Xr is fixed and positive and D -^ 4_). However, as it was pointed out in refs. [^ ,|6l|,fr^,[78|, 
such a theory is intrinsically unstable as the ground-state energy is unbounded from below. 
This is refiected, for instance, in the existence of tachyons in the theory [^,^,|7^, therefore 
the case with negative Aq is clearly inconsistent. 



The straightforward remedy for this situation was suggested by Bardeen and Moshe [ffZ 



They showed that the only meaningful (stable) 0{N) $^ theory in the large- A^ limit is that 
with Xr, Xq > 0. This is provided that we view it as an effective field theory at momenta 
scale small compared to a fixed UV cut-off A. The cut-off itself is further determined by 
( |3.55| ) because in that case (assuming rrir ^ A) 



Ar 



1 - i^H 



A2 

m'i 



(3.58) 



2 ^ 



which implies that for A,., Aq > we have A < m^ exp(^^|^). The case A = m^ exp{^Y~) 
corresponds to the Landau ghost [|7S| (tachyon pole [^|7^). For reasonably small A,-, A is 
truly hugeQ and so it does not represent any significant restriction. The following discussion 
will be confined to such an effective theory. 

3.3.3 The pressure 

The partition function Z has a well known path-integral representation at finite tempera- 
ture, namely 



^^For example, if A^ = 1 and to,, sa lOOMeV, we get A < lO^'^^MeV or equivalently A < lO^^^K (this is 
far beyond the Planck temperature - lO'^^K). 
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Z[T] = exp(m[T]) = V(j) exp{iS[(j);T]) 



S[^;T] 



d^x C(x) 



(3.59) 



c 



Here VH = —^Q is the so called Massieu function []2S|,^ and Jf^d^x = J^dxo Jyd^~^x 
with the subscript C suggesting that the time runs along some contour in the complex 
plane. In the real-time formalism, which we adopt throughout, the most natural version is 



the so called Keldysh-Schwinger one [p4|,|36[|, which is represented by contour in FIG. p.3| 



Ret 



Imt 



C, 



t ^ + CO 



Figure 3.3: The Keldysh-Schwinger time path. 



y 



Let us mention that the fields within the path-integral ( |3.59|) are further restricted by the 
periodic boundary condition (KMS condition) [^|36|,^ which in our case reads: 



4>a(ti -i/5,x) = 4>a(ti,^) ■ 



As explained in Section |3.2| , we can use for a pressure calculation the canonical energy- 
momentum tensor Q'^'^. Employing for G^'^(x) its explicit form ( |3.10| ) together with ( |3.14| ), 
one may write 
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(0r) = f /^(2gV-^^^(g'-m^))m)ii(9;r) +^9''' (\y,<Pim \, (3.60) 



where Dn is the Dyson-resumed thermal propagator [34,3^, i.e. 



^-(^' ^) = g^-m^iT)+^e + ^'"^ ^^'^ ' "^^^^^^ ^M^ " ^^'''^ 



Note that we have exploited in (|3.60|) the fact that the expectation value of 6(f''(x) is x 



independent. On the other hand, in (|3.61|) we have used the fact that m"^ is q independent. 
In order to calculate the expectation value of the quartic term in Eq. (|3.6CI|) , let us observe 



(c.f. (|3.59| )) that the derivative of Ht with respect to the bare coupling Aq (taken at fixed 
mo) gives 



.9m [T] 



d\o 8N 

which implies that 



y|^rf^x/('f^$^(0)'j V (3.62) 



TV ^ 2 



E<^^^])-w'-W- 



.a=l 



The key point now is that we can calculate ITI[T] in a non-perturbative form. (The latter 
is based on the fact that we know the Dyson-resumed propagator Dii(g;T) (see ( p.61|) .) 
Indeed, taking derivative of tYl with respect to ttiq (keeping Aq fixed) we obtain 






dml 2 Jc \^ "- " 2 J [271 
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(3.64) 



m[T;Xo;m'o]=f3VN / dm^ A^T(m'(T)) + m[T; Aq; oo] 



(3.65) 



Let us note that tXl[T] Aq; cxd] is actually zero[f] because tXl[T; Aq; rnl] has the standard loop 



expansion l3y,|72[ depicted in FIG.p^ 



16 



16 




Figure 3.4: First few bubble diagrams in the ITl expansion. 



It is worth mentioning that in the previous expansion one must always have at least one 
type-1 vertex [^. The RHS of FIG. p.4| clearly tends to zero for rrio ^ oo as all the (free) 
thermal propagators from which the individual diagrams are constructed tend to zero in this 



limit. The former result can be also deduced from the CJT effective action formalism 76 



or from a heuristic argumentation based on a thermodynamic pressure [^. Note that in 
the large-iV limit the fourth and fifth diagrams in FIG. p.4| must be omitted. 
The expectation value (|3.63| ) can be now explicitly written as 



^■^To be precise, we should also include in FIG j3.4| an (infinite) circle diagram corresponding to the free 
pressure |56|,|7^. However the later is Aq independent (although toq dependent) and so it is irrelevant for 



the successive discussion (c.f. Eq.( ^.63 )). 
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y^lm )=8N' / dml /T^^r^Im -^^^^ -, ] . (3.66) 



In fact, the differentiation of the proper self-energy in (|3.66|) can be carried out easily. 
Using ( |3.49| ), we get 



dXo Xo ' dXo dXo Xoil-XoM'^)' 



Prom Eq.( |3.49| ) it directly follows that 

dmJiT) 1 



dml (1 - Ao M't) ' 

which, together with the definition of A4t, gives 






,a=l 



^2(r) J (27r)^ 650/5 -1 (^2_^2+,,)2 

/■oo 

-8A^^ / dmlMTiml) 






dml 
AN^ Ml{ml{T)) , (3.67) 



where we have exploited in the last line the fact that A4'^{m,l -^ oo) = 0. Let us mention 
that the crucial point in the previous manipulations was that m,r is both real and momentum 
independent. Collecting our results together, we can write for the hydrostatic pressure per 
particle (cf. Eq.(p38D) 
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p(T)-P(o) = - j^^j^ {{eu) - {oHM) 






2 7 (27r)^-i VP-l)/ e5o/3_i 

6^{q -m^(0)) 



2 i (27r)^-i V(^-l) 

^ (S|K(T))-E^(m^(0))). (3.68) 



2A 



Applying the Green theorem to the last two integrals and eliminating the surface terms 
(for details see Appendix P) we find 



nT)-vio) ^ ^/(^^%^--^m) 



--^(S^(m,2(T))-S2(m^(0))) 

ZAq 



MT{ml{T))-M{mlm + -^ {^l{ml{T)) - T?{r^,m) , (3.69) 



where we have introduced new functions N'xi'nT'liT) and /^{m. 






^T{ml{T)) = i [ v:^ -1^ Oiq^ - mliT)) 



2 J (27r)^-i e*'^ - 1 
Ar(m2) = hm^Ar^(m^(T)). (3.70) 



Eg. ( p. 691) can be rephrased into a form which exhibits an explicit independence of bar 



quantities. Using the trivial identity: 
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^ (s|K(T))-s^K(o))) 

= ^ (StK^(T)) - S(m^(0))) (SHm^?(T)) + SK^lO))) 

ZAo 

= -^ {MAmliT)) +M{mlm) . (3.71) 

we get 



V{T)-V{Q)=MT{ml{T))-M{mlm + ^^ (A<^(^2^T)) + ^K'(O))) , (3.72) 



where 5m iT) = ml{T) — m^(0). The resuh (|3.72|) has been previously obtained by authors 



pT| in the purely thermodynamic pressure framework. 



3.4 Hydrostatic pressure in D = A (high— temperature 
expansion) 

In order to obtain the high-temperature expansion of the pressure in D = 4, it is presumably 
the easiest to go back to equation ( p.68| ) and employ identity ( p.71| ). Let us split this 



task into two parts. We firstly evaluate the integrals with potentially UV divergent parts 
using the dimensional regularisation. The remaining integrals, with the Bose-Einstein 
distribution insertion, are safe of UV singularities and can be computed by means of the 
Mellin transform technique. 

Inspecting (|3.68|) and ( p.71|) , we observe that the only UV divergent contributions come 



from the integrals: 

1 f d'^q 



+ (D^/(2^^' r(g^-m^(T))-5+(g^-m,^(0))) 
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^ 4 J '{2^ \q^ - ml{T) +it^ q'- m2(0) + ie J ' ^^'^^' 



which, if integrated over, give 



(3.73)= + ^ ' ^ \\ ^^ , (K(T))t-K(0))f 

(D-l)r(^)(47r)f ^^ '^^ ^^ 

+ ^""'^r/!'"^^ U<iT))^-' + K(0))--^ ) . (3.74) 

4(47r)^ V 



Taking the hniit D = A — 2e ^ A and using expansions 



r(-n + £) = L_^^ / ^ 



n\ 



\ k=l 



a"+" = a" 



(l + £lna + 0(£^)) , (3.75) 

(7 is the Euler-Mascheroni constant) we are finally left with 

(3-73)1.^. = - '"-'°]';r' In (^) + S'-HT) (mUT) + »J(0)) ^ . (3.76) 



The fact that we get finite result should not be surprising as entire analysis of Section 172 
was made to show that V{T) — V{0) defined via Q'^'^ is finite in D = 4. 

We may now concentrate on the remaining terms in ( p.68|) , the latter read (we might, 
and we shall, from now on work in D = 4) 
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3 7 (2^^ ^^;^^^3T^(^ -"^^(^)) + ^^y (2^d^;^^^^^(^ --.m)- (3.77) 

Our following strategy is based on the observation that the previous integrals have generic 
form: 



dx (x^ - 1)^ , (3.78) 



2 7r2 J^ '' '" ' e^y -I 

with i/ = 0, 1 and y = 771^(3. Unfortunately, the integral ( 3.78| ) can not be evaluated exactly, 



however, its small y (i.e. high-temperature) behaviour can be successfully analysed by 
means of the Mellin transform technique [^6|,|73[. Before going further, let us briefly outline 
the basic steps needed for such a small y expansion. 

The Mellin transform f{s) is done by the prescription [BB, 



5U-g3 



f{s)= / dxx'-^ f{x), (3.79) 

with s being a complex number. One can easily check that the inverse Mellin transform 
reads 



-| ptoo+a 

fix) = -77rT ds X-' f{s) , (3.80) 

K^'^) J-ioo+a 



where the real constant 'a' is chosen in such a way that f(s) is convergent in the neigh- 
bourhood of a straight line (— ioo + a, ioo + a). So particularly if /(x) = ^j^ry °^^ "^^^ ^^^ 
formula 1.3.19) that 
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fis) = Tis)as)y~ 



(Res > 1) 



(3.81) 



where ( is the Riemann zeta function {({s) = X]^i^~^)- Now we insert the Melhn trans- 
form of f{x) = ^xy_i to ( p.78| ) and interchange integrals (this is legitimate only if the 
integrals are convergent before the interchange). As a result we have 



OO -| 

dx q(x) 

^ ^ ^ e^?^ - 1 



-° ds 

ha ^(27r; 



Tis)as)y-'g{l - s) 



(3.82) 



with g{x) = 9{x — 1) (x^ — 1)2. Using the tabulated result ( ||83|; formula 6.2.32) we find 



g{l-s) = ^B{-u-l + ls;l + iy) (Res > 2 + 2z/) , (3.83) 

with B{ ; ) being the beta function. Because the integrand on the RHS of ( |3.82| ) is analytic 
for Res > 2 + 2z/ and the LHS is finite, we must choose such a that the integration is 
defined. The foregoing is achieved choosing a > 2 + 2z/. Another useful expressions for 
^(1 - s) are ( H; formula 1.2.34 or 1.2.37) 



^(1-s) = 5(| + z/;-2-2z/ + s) 2Fi[-\-u;-2-2u + s;-l-u + s;-l] 



25+- 5(1 + z/; -2 - 2z/ + s) ai^ihl - z/; | + i^; -| - z/ + s 



■11 

' 2i ' 



where 2F1 is the (Gauss) hypergeometric function [B3]. Using identity 



22x-l 

T{2x) = ^^T{x)T{x + l), 



we can write 
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(3.82) 



4v^ 



joo+a 



ds 
i(27r) 



r(|^)C(^)Q2/)"^r(-^-i + H 



(3.84) 



The integrand of ( |3.84| ) has simple poles in s = — 2n (n = 1, 2, . . .), s = 1, s = —2n + 
2z/ + 2 (n = 0, 1, . . . , z/) and double pole in s = 0. An important point in the former pole 
analysis was the fact that C{s) has simple zeros in —2m {m > 0) and only one simple pole 
in s = 1. The former together with identity 



shows that no double pole except for s = is present in 



Now, we can close the 



contour to the left as the value of the contour integral around the large arc is zero in the 
limit of infinite radius (c.f. [^ and ||8^; formula 8.328.1). Using successively the Cauchy 
theorem we obtain 



4v^ (3.82) 

r(| + ^) 



n=0 

oo 



2n-2u-2 



TX 



-2n+2.+2/_^^^), (_1)"|5_ 



2n+2j^+2| 



n! (-2n + 2z/ + 2)! 2^''-^^-^ 
7r-2" {2n)\ C(l + 2n) (-l)"+'^+i 



n=l 



n\in + l + u)\ 24"-i 



+ y 



_i 7r(- 1)^^+1 (z/ + l)! 22-^+3 2(-l) 



(2Z/ + 2)! 



[y + 



9^f(^)^^-^|i}.(3-) 



where -Bq's are the Bernoulli numbers. Let us mention that for C,{2n + 1) only numerical 
values are available. 
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Inserting (j^Wj ) back to (^), we get for V{T) - V{0) 



1 Sm^(T) 

V{T) - V{0) = (3.76) + - h{mr{T)) + -^ h{mr{T)) 



— - m^[T) 



90 24 V 2 / 4 7r V 3 4 



^ m,^(T) m,^(0) / / m,(0) \ ^ _ l^ <(0) 



32 7r2 V \T An J ' 1) 128 vr^ 

m2"+2(T) 7r-2"-2 (2n)! C(l + 2n) (-l)i+" 



^ (m,^(T) - (^ 5m^(r) 



n=l 



T2"n! (n + 2)! 24"+4 

(3.86) 



Note that (|3.76|) cancelled against the same term in | l2{mr{T)) + "^^ ' Io{mr{T)). One 
can see that (|3.86|) rapidly converges for large T, so that only first four terms dominate at 
sufficiently high temperature. The aforementioned terms come from the poles nearby the 
straight line {—ioo + a, ioo + a) (the more dominant contribution the closer pole). It is a 
typical feature of the Mellin transform technique that integrals of type 

dxg{x) — — -, 
e^y — 1 

can be expressed as an expansion which rapidly converges for small y (high-temperature 

expansion) or large y (low-temperature expansion)^. 

For a sufficiently large T we can use the high-temperature expansion of 5m? {T) found 

in Appendix |C[ Inserting ( |U.6| ) to ( p. 861 ) we obtain 



^"^By the same token we get the low-temperature expansion if the integral contour must be closed to the 
right. 
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^ ' ^ ' 90 24 127r 



8 V 144 247r IGyr^ / \ ^ v ; i j.^^ j j v ; 



Up to a sign, the result ( p.87| ) coincides with that found by Amehno-Cameha and Pi [35| for 
the effective potentiaQ Actually, they used instead of the N -^ oo limit the Hartree-Fock 
approximation which is supposed to give the same Veff as the leading 1/A^ approximation 



Note that the discussion of the mass renormalisation in Section 3.1 can be directly 
extended to the case when 771^(0) = (this does not apply to our discussion of A^!)- Latter 
can be also seen from the fact that (|3.86| ) is continuous in mr{0) = (however not analytic). 
The foregoing implies that the original massless scalar particles acquire the thermal mass 
■m'^{T) = 6m'^{T) . From ( p.86| ) one then may immediately deduce the pressure for massless 



fields $a in terms of 6m{T). The latter reads 



V(T)-V(0) = ^ L^^ + ^ L^^ 

^ ^ ^ ^ 90 48 48 TT 



^ (<5m(T))2"+4 vr-2"-2 (2n)! C(l + 2n) (-l)"+i 



n=l 



T^'^ (n-1)! (n + 2)! 2^^ 



^■^Let us remind \75,7t, BSl that from the definition of Veff the thermodynamic pressure is —Veff- In 



order to obtain (3.87) from Veff in pq], one must subtract the zero temperature value of Veff and restrict 



oneself to vanishing field expectation value and positive bare mass squared. 
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This result is identical to that found by Drummond et al. in ||61|| . 



A noteworthy observation is that when the energy of a thermal motion is much higher 
then the mass of particles in the rest, then the massive theory approaches the massless one. 
This is justified in the first (high-temperature dominant) term of ( |3.86D and (|3.88|) . This 
term is nothing but a half of the black body radiation pressure for photons |^^ (photons 
have two degrees of freedom connected with two transverse polarisations). One could 
also obtain the temperature dominant contributions directly from the Stefan-Boltzmann 
law |^,|3^,|5^ for the density energy (i.e. (B°°)). The formal argument leading to this 
statement is based on the noticing that at high energy (temperature) the theory at hand 
is (classically) approximately conformally invariant, which in turn implies that the energy- 
momentum tensor is traceless [^]. Taking into account the definition of the hydrostatic 
pressure ( p.8| ), we can with a little effort recover the leading high-temperature contributions 
for the massive case. 

3.5 Conclusions 

In the present chapter we have clarified the status of the hydrostatic pressure in (equilib- 
rium) thermal QFT. The former is explained in terms of the thermal expectation value 
of the 'weighted' space-like trace of the energy-momentum tensor Q^^ . In classical field 
theory there is a clear microscopic picture of the hydrostatic pressure which is further 
enhanced by a mathematical connection (through the virial theorem) with the thermody- 
namic pressure. In addition, it is the hydrostatic pressure which can be naturally extended 
to a non-equilibrium medium. Quantum theoretic treatment of the hydrostatic pressure 
is however pretty delicate. In order to get a sensible, finite answer we must give up the 
idea of total hydrostatic pressure. Instead, thermal interaction pressure or/and interaction 
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pressure must be used (see ( |3.38| ) and ( p.39| )). We have established this result for a special 
case when the theory in question is the scalar $^ theory with 0{N) internal symmetry; but 
it can be easily extended to more complex situations. Moreover, due to a lucky interplay 
between the conservation of B^'^ and the space-time translational invariance of an equilib- 
rium (and T = 0) expectation value we can use the simple canonical (i.e. unrenormalised) 
energy-momentum tensor. In the course of our treatment in Section |3.2| we heavily relied on 
the counterterm renormalisation, which seems to be the most natural when one discusses 
renormalisation of composite Green's functions. To be specific, we have resorted to the 
minimal subtraction scheme which has proved useful in several technical points. 

We have applied the prescriptions obtained for the QFT hydrostatic pressure to ^^ 
theory in the-large N limit. The former has the undeniable advantage of being exactly 
soluble. This is because of the fact that the large-A^ limit eliminates 'nasty' classes of 
diagrams in the Dyson-Schwinger expansion. The surviving class of diagrams (superdaisy 
diagrams) can be exactly resumed, because the (thermal) proper self-energy S, as well as 
the renormalised coupling constant A,, are momentum independent. We have also stressed 
that the 0{N) $'^ theory in the large- A^ limit is consistent only if we view it as an effective 
field theory. Fortunately, the upper bound on the UV cut-off is truly huge, and it does 
not represent any significant restriction. For the model at hand the resumed form of the 
pressure with mr{0) = was firstly derived (in the purely thermodynamic pressure context) 
by Drummond et al. in [^ . We have checked, using the prescription ( |3.38|) for the thermal 
interaction pressure, that their results are in agreement with ours. The former is a nice 
vindication of the validity of the virial theorem even in the QFT context. 

The expression for the pressure obtained was in a suitable form which allowed us to take 
advantage of the Mellin transform technique. We were then able to write down the high- 
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temperature expansion for the pressure in D = 4 (both for massive and massless fields) in 
terms of renormahsed masses mr{T) and m,.(0). We have exphcitly checked that all UV 
divergences present in the individual thermal diagrams 'miraculously' cancel in accordance 
with our analysis of the composite operators in Section |3.2| . 



Chapter 4 

Pressure in out— of— equilibrium media 



In recent years significant progress has been made in understanding the behaviour of QFT 
systems away from thermal equihbrium. Motivation for the study of sucli systems comes 
both from the early universe as well as from the quark-gluon plasma (deconfined phase 
quarks and gluons). Non-equilibrium effects are expected to be relevant in the relativistic 
heavy-ion collisions planed at RHIC and LHC in the near future |]86H89|. 

One of the significant physical variables, in the context of non-equilibrium QFT, is 
pressure. Pressure, as an easily measurable parameter^, is expected to play an important 
role in a detection of phase transitions. This is usually ascribed to the fact that the pressure 
should exhibit a discontinuity in its derivative(s) when the local phase transition occurs. 
The aforementioned has found its vindication in solid state physics and fluid mechanics, 
and may play a crucial role, for instance, in various baryogenesis scenarios. 

It is well known that for systems in thermal equilibrium, the pressure may be calcu- 
lated via the partition function p^,p6|,|62|. However, this procedure cannot be extended to 



off-equilibrium systems where is no such thing as the grand-canonical potential. In this 



"'^In this connection we may mention the piezo resistive sihcon pressure sensors used, for instance, in 



superfluid hehum [QQ or neutron (X-ray) diffraction technique used in sohd state physics ]9^, 93 

96 
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chapter we consider an alternative definition of pressure, based on tlie energy-momentum 
tensor. Tliis, so called, hydrostatic pressure is defined as the space-like trace of the energy- 
momentum tensor (see Chapter |^), and in equilibrium, it is identified with the thermody- 
namical pressure via the virial theorem |36|. There are several problems with the validity 
of this identification on the quantum level, indeed gauge theories suffer from a conformal 



(trace) anomaly and require special care ||36|. However, we will avoid such difficulties by 



focusing on a scalar theory which is free of the mentioned complications []36| , p4| . The major 



advantage of defining pressure through the energy-momentum tensor stems from the fact 
that one may effortlessly extend the hydrostatic pressure to non-equilibrium systems (for 
discussion see Chapter H) . 

The aim of this chapter is to provide a systematic prescription for the calculation of the 
hydrostatic pressure in non-equilibrium media. This requires three concepts; the Jaynes- 
Gibbs principle of maximal entropy, the Dyson-Schwinger equations and the hydrostatic 
pressure. In order to keep our discussion as simple as possible we illustrate the whole 
procedure on a toy model system, namely the 0{N) $^ theory. The latter has advantage 
of being exactly solvable in the large- A^ limit, provided that we deal with a translationally 
invariant medium. As a result the hydrostatic pressure may be expressed in a closed form. 

In order to provide meaningful results also for readers not entirely familiar with the 
off-equilibrium Dyson-Schwinger equations and the Jaynes-Gibbs principle of maximal en- 
tropy, we briefly summarise in Section ^]l| the basic essentials. (Detailed discussion of the 
equilibrium case may be found in Appendices [A.1| - |A.3| or in Refs. ^,|9^,|96|.) As a byprod- 
uct we construct the generalised Kubo-Martin-Schwinger (KMS) conditions. Section [4.2| is 
devoted to the study of the (canonical-) energy momentum tensor in the 0{N) $^ theory. 
If both the density matrix and the full Hamiltonian are invariant under 0{N) symmetry one 
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obtains Ward's identities in a similar manner as in equilibrium. We show how these drasti- 
cally simplify the expression for the pressure. In Section [4.3| we concentrate our analysis on 
the large- A^ limit. In this setting we derive a very simple expression for the pressure - pres- 
sure depends only on two-point Green functions. Section \i.4\ then forms the vital part of 
this paper. Owing to the fact that the infinite hierarchy of the Dyson-Schwinger equations 
is closed (basically by chance) we obtain simple equations for two-point Green functions - 
Kadanoff-Baym equations. These are solved exactly for three illustrative density matrices 
p. We choose deliberately translationally invariant p's. The reason is twofold. Firstly, for a 
non-translationally invariant medium one must use the improved energy momentum tensor 
instead of the canonical one |^4[. This is rather involved and it will be subject of our future 



work. Secondly, the Kadanoff-Baym equations turn out to be hyperbolic equations whose 
fundamental solution is well known. As a result we may evaluate, for the density matrices 
at hand, the hydrostatic pressure explicitly. The chapter ends with a discussion. 

4.1 Basic formalism 

The key object of our interest is the energy-momentum tensor Q^u{x). A typical contribu- 
tion to Q^u{x) can be written as 

Here $ is a field in the Heisenberg picture and V^. stands for a corresponding differential 
operator. Since P^^$(x) and V^^^{x) generally do not commute for i ^ j, one must 
prescribe the ordering in Q^u- Our strategy is based on the observation that one can 
conveniently define such ordering via the non-local operator 
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\\m^V{td{^y)Tmxi) ■ ■ ■ $(x„)} , (4.1) 

where V{id{^}) is just a useful short-hand notation for V^-^^T)^^ . . .V^^, and T* is defined 
in Section p.2| . It is clear that both T* and T coincide if all the arguments Xi are different, 
so (T* — T)(. . .) is an operator with a support at the contact points. The T* ordering is in 
general a very useful tool whenever one deals with composite operators. In the sequel we 
shall repeatedly use this fact. 

4.1.1 Off— equilibrium Dyson— Schwinger equations 

Let us now briefly outline the derivation of the Dyson-Schwinger equations for systems 
away from equilibrium. For simplicity we illustrate this on a single scalar field $. 

We start with the action S[^] where $ is linearly coupled to an external source J{x). 
Working with the fields in the Heisenberg picture, the operator equation of motion can be 

written as 

^^ .[$ = <|.-^] + J(x) =0, (4.2) 



5$(x) 
where the index J emphasises that $ is implicitly J-dependent. This dependence will be 

made explicit via a unitary transformation connecting ^"^ (governed hj H — J$) with $ 

(governed by H). If J{x) is switched on at time t = U we have 

^J{x) = B-\t;U)^{x)B(t-U) 

= B{U-tf)B{tf-t)^{x)B{t-U) 

= Tc U{x) exp {i j dS J{ymy))\ . (4.3) 
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Here we have used the fact that 



B{t,;t2)=T[exp{t I d'^y J{yMy)) ] t, > h 



(4.4) 



The close-time path C runs along the real axis from ti to tf {tf is arbitrary, tf > ti) and 



then back to tj. With this setting we can rewrite ( [4. 21) as 






.f\5S[^-- 



J±\ expitjd'y (J+(l/)$+(l/) - J_(2/)<^-(2/)))) = 0: 



(4.5) 



where, as it is usual |^,0], we have labelled the field (source) with the time argument on 
the upper branch of C as $+ ( J+) and that with the time argument on the lower branch of C 
as $_ {J~)- Introducing the metric (cra)^/? (o"3 is the usual Pauh matrix and a, /3 = {+; — }) 
we can write J+$+ — J_$_ = Jq (0-3)"^ $^ = J" {0-3)^13 ^^. For the raised and lowered 
indices we simply read: $+ = $+ and $_ = — $^ (similarly for Jo). Taking Tr(p . . .) with 
p = p[$, 9$, . . .] being the density matrix, we obtain 



63 



Z[J] 6^{x) 



<l>"(x) 



6 



—I 



SJJx) 



Z[J] 



-rix) 



(4.6) 



with Z[J] = Tr {p T^exp [i J^d^y J(?/)$(|/))} being the generating functional of Green's 
functions (which in the non-equilibrium context we shall denote az G). Employing the 



commutation relation: —ijj-Z = Z 
form, namely 



— t 



5 J a 



), we may recast ( [4.6|) into more convenient 



r(x) 



63 
6^{x) 



r{x)+t d^yG'l{x,y){asf^ 



5<P^{y) 



(4.7) 
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where the symbol I indicates the unit. As usual, the mean field, (pa, is defined as the 
expectation value of the field operator: 0q,(x) = (<l'a(x)). Defining Z[J] = exp{iW[J]), 
two-point Green functions are then defined as Gai3{x,y) = —5'^W/5J"{x)5J^{y). Setting 
J in ([4.6|) - ([4.7|) to (i.e. physical condition) we obtain a first out of infinite hierarchy 
of equations for Green's functions. Successive J variations of (|4.6| )- (|4.7| ) generate higher 
order equations in the hierarchy. The system of these equations is usually referred to as 
the Dyson-Schwinger equations. 

For the future reference is convenient to have the expression for the effective action 
r[0c]- This is connected with W[J] via the Legendre transform: 



T[<P,] = W[J]- / d'yJ{y)<P,{y). (4.8) 

Jc 

Following the previous reasonings, one can easily persuade oneself that the expectation 
value of Q^^ reads 

{Q^^{x)) = {B'^^mx)]) = Q^^[cP+{x)+ijd'yQ+p{x,y) {a.f^ ^^]^- (4-9) 

We have automatically used the sub-index '+' as the fields involved in Q^^ have, by defi- 
nition, the time argument on the upper branch of C. 

4.1.2 The Jaynes— Gibbs principle of maximal entropy 

In this section we would like to review the Jaynes-Gibbs principle of maximal entropy (also 



maximum calibre principle) pj,|2^, 0, ^, which we shall employ in the following. The 
formalism is a generalisation from an ordinary Gibbs's principle of maximal entropy to 
the systems out of equilibrium. The objective of the principle is to construct the 'most 
probable' density matrix which fulfils the constraints imposed by experimental/theoretical 
data. 
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The standard rules of statistical physics allows us to define the expectation value via 
the density matrix p as 

(...)=Tr(p...), (4.10) 

with the trace running over a complete set of physical states describing the ensemble in 

question at some initial time ti. The specification (or reasonable approximation) of the 

density matrix is thus crucial for determining the macroscopic properties of a given system. 

The usual approaches [p^lOO|| trying to determine p start with the Schrodinger picture. 



The merit of this procedure is that one transforms the whole time dependence of the expec- 
tation value into the density matrix itself. Thus the time evolution of such an expectation 
value is equivalent to the time evolution of p, and we need not to solve separate equations 
for each observable^. 

Suppose that the system is prepared in such a way that at time t the probability of the 
system being in the state \4'n]'t) is Pn(t). The density matrix has then the standard form 

Ps(t) = 5^P„(t) |Vn;i)(V'n;t|. (4.11) 

n 

(by ^^ we mean; sum over discrete spectrum and integrate over continuous one). We should 
stress that the ensemble {\ipn',t)} in ( [4.1 1| ) not necessary consist of mutually orthogonal 
states, although the density matrix can always be formally diagonalised by selecting its 
eigenbasis (polar basis). Applying the Schrodinger equation to ( [4. Ill ), the evolution of ps 
reads 

*^ = .|M*),//l + M^. (4.12) 

where 



^This advantage docs not seem to be relevant for Green's functions with different time arguments. 
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— ^^ = 2^^^ l^n;t)(^n;t|. (4.13) 



For a closed Hamiltonian system p„(t) cannot be changed, and so ( [4.12| ) reduces to the 
celebrated von Neumann-Liouville equation 

^-^ = i[Ps{t\H]. (4.14) 

Let us mention that the vanishing of the time derivative of Pn{t) implies that the von 
Neumann-Gibbs entropy Sq = — Tr(pslnps) is time independent. Indeed, in the latter case 
the evolution of ps can be formally written (c.f. ( [4.12| )) as 

Ps{t) = U{t-t')ps{t')U{t-t')-\ (4.15) 

where the evolution operator is determined from the Schrodinger equation 

i^U(t-t') = H{t)U(t-t'), U{t',t') = l, 
and so using the property of a trace operation we have, 

Scit') = -Tr(p,(t')lnp.(t')) = -Tr(p,(t)lnp,(t)) = Scit) . (4.16) 

The time dependence of Pn(t) reflects the fact that the system in question interacts with 
an exterior - a heat bath (whose structure and dynamics are usually unknown). In order 
to determine Pn(t) one would need to resort to some physical model [|100|| describing the 
dynamics of the environmental system - task by no means easy. In order to avoid these 
difficulties it is customary to apply a stochastic description of a system-bath interaction 



101, 102], or equivalently use an irreversible "master equation" |102| for p. 
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To keep our ideas simple, we shall from the very outset assume that p„ is constant. This 
means that either no external environment interacts with our system or that the interaction 
is reflected only in the Hamiltonian H (e.g. effective external fields or time dependences) 
but the probability of the population of any state stays unchanged - adiabatic interaction. 

Solving the von Neumann-Liouville equation would be a formidable task. One usually 
resorts either to variational method [[10CI| , |103| | with several trial p's, or one may recast ( [4.14| ) 



into an infinite hierarchy of integro-differential equations for two-particle, etc. distribution 
functions [^ (i.e. Bogoliubov-Born-Green-Kirkwood-Yvon or BBGKY hierarchy). In 
the latter approach one hopes that one can device an effective truncation of the hierarchy 
allowing to close the system of equations, and then solve it (usually perturbatively). How- 
ever, in both aforementioned approaches we must specify the initial value data in order 
to solve uniquely the evolution equation(s). This is quite delicate task, and one usually 
uses the most simple (and not always physically well motivated) choices of p (e.g. Gibbs 
(grand-) canonical distribution ||100| |, Gaussian distribution [ P5| , p5| ,pU|). 

Rather than following the previous path, we shall use the Heisenberg picture instead. 
This seems to be more suitable for our purpose. In our particular case the polar form of 
Ph reads 

PH = ^Pn\i^n){lpn\- (4.17) 

n 

In the following we shall denote pn simply as p. In order to determine p explicitly we shall 



resort to the Jaynes-Gibbs principle of maximal entropy [El, 12^, PT], P3 . The latter allows 



one to construct p in a way which naturally extends the original Gibbs prescription [pi,P^ 
which refers only to equilibrium. The basic idea is 'borrowed' from the information theory. 
Let us assume that we have criterion of how to characterise the informative content of p. 
The most "probable" p is then selected out of those p which are consistent with 'whatever' 
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we know about the system and which have the last informative content (Laplace's principle 
of insufficient reasons). Consistency with anything known about p must be kept; while to 
chose more informative p is to presume an extra information which we do not control, and 
to make unjustified implicit assumption concerning the information we do not know about. 
It remains to characterise the information content (measure) ■^[p] of p. This was done 
by C.E.Shannon ||104|| , L.Brillouin ||105|| and L.Szilard ||106|| with the result 



^[p] =Tr(plog2p) 

The density matrix is then chosen to minimise '6[p]- It is quite surprising that up to 
a (negative) multiplicative constant ■^[p] coincides with 5'g'[p], and thus the principle of 
insufficient reasons basically turns out to be the maximal entropy principle^]. Note that 
no assumption about the nature of p was made; namely there was no assumption whether 
p describes equilibrium or non-equilibrium situation. (We were interested only in the 
information content of p, and not in the underlying dynamics of a system). To put more 
flesh on the bones, let us rephrase the former into more physical language ||2^,|25|. What we 
actually need to do is to maximise Sq subject to the constraints imposed by our knowledge 
of expectation values of certain operators Pi[$, 9$], . . . , Pn[^, d^]- This yields a density 
matrix p which incorporates the fact that all the quantum states which are permitted 



by the constraints have equal probabilities [^. In contrast to equilibrium, all Pfc[---]'s 



may be operators which are not the constants of the motion (both the position and time 

dependences are allowed) . So namely if one knows that 

•^ It Appendix O we show that the Shannon entropy (which is proportional to Sg) equals (in base 2 of 
logarithm) to the expected number of binary (yes/no) questions whose answer take us from our current 
state of knowledge to the one certainty. So, the bigger Sq, the more questions one has and consequently 
one possess a higher ignorance about a system. 
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{Pk[^,d^])=gk{xi,X2,...), (4.18) 

the entropy inaxiinalisation leads to 

with the 'partition function' 

Z[Ai, . . . , A„] = Tr I exp(- J2 [ H ^"^^^ ^^^^i' • • •)^^['^' ^^^H ' (^'^O) 



It is possible to show that the stationary solution of 5*^ is unique 1^^ and maximal. The 
latter goes on the account of the fact that Sg is a concave functional (see Appendix 0). 
Both in ( [4.19| ) and ( [4.20| ) the time integration is not either present at all (so fk is specified 



only in the initial time tj), or is taken over the gathering interval (— r, ti) (i.e. as J ' dt . . .). 
If, instead, one has certain partial knowledge about the expectation value of Pk[^,d^] at 
some discrete times prior ti, the corresponding integration must be replaced by a discrete 
summation. Note that if we have no prior knowledge about the system, then p = 1/W, 
where W is the number of accessible quantum states. Equivalently we may say that the 
no prior restrictions mean our total ignorance about the system and as a consequence we 
must affiliate with each quantum state an equal probability. 

The Lagrange multipliers A^ might be eliminated if one solves n simultaneous equations 



S InZ 
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Using the definition of tlie von Neumann-Gibbs entropy togetlier witli ( 4.19| ) we get 



\max 

n 



SG[9u---,9n]\max = -Tr(plnp)|, 

= lnZ[Ai,...,A„] + ^ /"jJrfV, A,(xi,...)(7.(a:i,...)- (4.22) 

i=i -^ j 

So one may view Sq max as tfie Legendre transformation of ln2^. In the equihbrium case the 

former is the standard relation between entropy and the grand-canonical potential. Having 

( |4.22| ), the explicit solution of ( [4 .211 ) may be formally written as 



6gk{xi,...) 
Now, in order to reflect the density matrix ( [4.19|) in the Dyson-Schwinger equations, we 
need to construct the corresponding boundary conditions. This may be done quite straight- 
forwardly. Using the cyclicity of Tr(. . .) together with the relation 



e^^ 



Be — y —:Cn, Co — B,Cn — [A, Cn-l] , 

^-^ n! 

n=0 

we can write the generalised KMS conditions for the n-point Green function as 

mXi) . . . ^Xn)) = ($(X2) . . . $(X„)<I>(X1)) + J2 ^y ^•^(^2) . . . HXn)Ck{Xi)) , (4.24) 

fc=l 

where A = ln(p), B = $(xi), xio = tj. So for the two-point Green functions we haveQ 

G+_(,T, y) = G_+(a;, y) + V -Tr(p Hx)Cn{y)) . 



71=1 



''in special cases when p = |0)(0| oi p — e ^^ /Z{(3) the boundary conditions are the well known 
Feynman and KMS boundary conditions, respectively. 
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In this chapter we aim to demonstrate that conditions ( |4.24| ) together with the causahty 
condition are sufficient to determine Green's functions uniquely. 



4.2 The 0{N) <|)^ theory 

The 0{N) $^ theory is described by the bare Lagrangian 

N ^ / N 



2 



a=l \a=l 



2 



The corresponding canonical energy-momentum tensor is O^*^ = ^^d^^"'d'^(^"' — g^'^ C, and 
from Eqs. (|4.1| ) and ( [4 .91 ) its expectation value is 



er($+(x)) = (er(x)) 

c c 

+ w^ E { ((<^+(^))' + ^^+%(^. ^)) ((</'+ (^))' + «++(^' ^)) 

+ 2 (</.;(x)0^(x) + 2Gf+(x, x)) zG;^+(x, x) } 

+ terms with T^ and T^ . (4.25) 

Before proceeding further with our development, we want to show how one can significantly 
simplify Eq. ( |4.25| ) provided that both the density matrix and the Hamiltonian are invariant 
under rotations in the iV-dimensional vector space of fields. This situation would occur if 
the system was initially prepared in such a way that no field $" was favoured over another. 
The fact that p is invariant under 0{N) transformations means that 

f/(e)p[<l>, 9$, . . .]U-\e) = p[<^, 9$, . . .] , (4.26) 

where the fields $" transform under A^-dimensional rotations: U{e)^°'U~^{e) = i?~^(e)"^$^ = 
[exp(ejTj)]'*''$^, where R{e) is the rotation matrix in the A^-dimensional vector space, and 
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the generators Tj are real and antisymmetric N x N matrices. It is obvious that the previ- 
ous relation for ^'^ can be satisfied for all times only if the full Hamiltonian, which governs 
the evolution of $", is also invariant against the A^-dimensional rotations. 

Let us now consider the generating functional Z[J] corresponding to the 0{N) symmet- 
ric density matrix. Employing the cyclic property of Tr(. . .) together with the infinitesimal 
transformation, 6R{e) = 1 + eiTi, we obtain that the variation of Z must vanish. The latter 
implies the following (unrenormalised) Ward's identities: 



.4,. „,.jmj] f ,, ,,,..-.mj] 



/''''^«^jiW)^Jc''''^"^sm- ^'-''^ 



Taking successive variations with respect to source J, we obtain the following results (see 
also |]96|): n-point Green functions with n odd vanish, while for n even (n = 2k, k = 1,2, . . .) 
one has 



Gr.^«;S(^l'---'^2.)= J2 U^a.a,G^"t^..^Jx„...,X,,), (4.28) 

all pairings i<j 



where G*^^'^) is a universal 2A;-point Green function. Similar results can be obtained for 

■paia2...a2fc ( \ 

Finally note that these results enable the expression for the expectation value of the 
energy momentum tensor to be simplified significantly to 



N 



t-g^'^{d^^dl G++(x,y)U^-m^ G++(x,x)) 
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+ terms with T^^^ . (4.29) 

In the rest of this chapter we shall confine ourselves only to situations where both p and H 
are 0{N) invariant. 

4.3 The large-A^ limit 



Let us now examine behaviour of (|4.29| ) to the order 0{1/N). For this purpose it is im- 



portant to know how either G*^"^ or F*^"^ behave in the A^ -^ oo limit. At T = or 



m 



equilibrium the Feynman diagrams are available |6^,|7J] and the corresponding combina- 
torics can be worked out quite simply. On the other hand, the situation in off-equilibrium 
cases is more difficult as we do not have at our disposal Wick's theorem. One may devise 



various diagrammatic approaches, e.g. kernel method |^, cumulant expansion |5^, cor- 



relation diagrams 0, etc. Instead of relying on any graphical representation (as we done 
in Section |3.3| ) we show that for both equilibrium and off-equilibrium systems, the situa- 
tion may be approached from far more general standpoint using only Ward's identities and 
properties of F and W. 

In order to find the leading behaviour at large A^ it is presumably the easiest to consider 



the Legendre transform (|4.8| ). The explicit A^ dependence may be obtained by setting 
(j)'^ = 0, which implies J'^ = J for all the group indices. Eq. (|4.8| ) then indicates that both 
F[0] and W[J] are of order A^. If we expand F[0"] in terms of (j)'^ taking into account Ward's 
identities we get 

F[0] = F[0] + In [ d^x d^y F(2)(a;, y) (j){x)(j){y) 



^ Jc 

+ ^N^ f d\ dS d^z d^q F(^) (x, y, z, q) 
4! Jc 

X <P{x)<P{y)<P{z)<P{q) + ■ ■ ■ . (4.30) 
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Since the LHS of ( [4.30| ) is of order A^, F^^^ must be of order A^'', T^'^^ of order A^ ^, and, in 
general, F^^"-* must be of order N^~^. This suggests that in the expression for the energy- 



momentum tensor (|4.29|) , terms containing F^'*^ can be ignored. The above argument can 



be repeated in a similar way for W. 

Hence, collecting our results together, the expectation value of the energy-momentum 
tensor to leading order in A^ may be written as 

.N 



(er(a;)) = zNd>id;G+^{x,y)l^^--Xog^''{G++{x,x)f 



- ^^9'\d^^dl G++(a;,2/)L^^ - ml G++(a;,a;)) . (4.31) 

This result is surprisingly simple: the expectation value of the energy-momentum ten- 
sor, and thus the hydrostatic pressure, is expressed purely in terms of two-point Green's 
functions. The latter can be calculated through the Dyson-Schwinger equations ( [4.7|) . Fur- 
thermore, these equations have a very simple form provided that both the large-A^ limit 
and Ward's identities are applied. If we perform a variation of ( |4.7|) with respect to Jfiiy) 
we obtain, to order 0{1/N), the following Dyson-Schwinger equations for two-point Green 
functions: 

( D + mo + — ^ G„„(x, x) J Qap{x, y) 

= -6{y-x){a^)^p. (4.32) 

These dynamical equations for two-point Green functions are better known as the Kadanoff- 
Baym equations |]40[| . 



Let us mention one more point. The generalised KMS conditions for G±=p are signifi- 
cantly simple in the large- A^ limit. This is because in sum ( |4.24|) only terms of order 0{N^) 
contribute. This implies that only quadratic operators Pj[$,(9$] in the density matrix are 
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relevant. As a result, the Jaynes-Gibbs principle naturally provides a vindication of the 



popular Gaussian Ansatz p8|-p.0CI||. 



4.4 Out— of— equilibrium pressure 

The objective of the present section is to show how the outlined mathematical machinery 
works in the case of the hydrostatic pressure. In order to gain some insight we start 
with rather pedagogical, but by no means trivial examples; translationally invariant, non- 
equilibrium density matrices. We consider the more difficult case of translationally non- 
invariant density matrices in our future work. 

4.4.1 Equilibrium 

As an important special case we can choose the constraints 

(Pfc[$, d^])\u =gk = constant , (4.33) 



where tj is arbitrary. Eq.( |4.33D then implies that P^'s are integrals of motion. Since in the 
finite volume systems the spatial translational invariance is destroyed, the only integral of 
motion (apart from conserved charges) is the Hamiltonian. Thus the system is in thermal 



equilibrium and the laws of thermodynamics [^,^ prescribe that g = Jq dT'CyiT') (C- 



'V 



is the heat capacity at constant volume V and T is temperature). Eq. ( [4.23|) determines 



the Lagrange multiplier; A = 1/T = (3. The density matrix maximising the Sq is then 
the density matrix of the canonical ensemble: p = — ^,~f . Due to the translational 
invariance of p, the Kadanoff-Baym equations read 

(n, + ml{T)) G^p{x -y) = -5{x - y){a:,)^p , (4.34) 
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where the temperature-dependent renormahsed mass is (see, for example |^,|73|); m'^iT) 



ml + if G„„ 



(0). The corresponding KMS boundary condition is 



G+_ (x, U; y , 0) = G_+(x, U - i(3] y, 0) 



(4.35) 



Because mr{T) is a spatial constant, a Fourier transform solves equations ( |4.34| ). The 
solutions of ( |4.34| ) subject to condition ( [4.35| ) are then the resumed propagators in the 
Keldysh-Schwinger formalism 






k"^ — ml ± ie 



+ 2'KfB{%\)6{k'-ml) 



2Tx{e{Tk,) + fB{\k,\)}5{k''-ml) 



(4.36) 



with /_b(x) = (exp(/5x) — 1) ^ being the Bose-Einstein distribution. 



Now, the total hydrostatic pressure in D dimensions is classically defined as [p6| , |74| , pU7 

1 



p{x,T) 



(D-l 



■(e^(x)). 



Because Q'^^ is a composite operator, a special renormalisation is required p6|, |65|, |66|, pi 



As we have shown in [^, for translationally invariant systems the renormahsed pressure 



coincides with the, so called, thermal interaction pressure V (see Eq.( 3.38 )). The latter 
reads 



V{T) =p{x,T)-p{x,0) 



1 



(D 



YA(^'c^) - meUo)} . 



(4.37) 



Let us remind that the energy-momentum tensor O^'^ need not to be the canonical one 
(however, the canonical one is usually the simplest one), c.f. Section p.2| . As a second point 
we should mention that prescription ( ^.37| ) retains its validity for non-equilibrium media 
as well. This is because the short-distance behaviour of G++{x, y), which is responsible for 
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the singular behaviour of Q^'^ , comes from the particular solution of the Kadanoff-Baym 
equation ( [4.34| ). The latter can be chosen to be completely independent of the boundary 
conditions (actually it is useful to chose the Feynman causal solution). On the other hand, 
the homogeneous solution of (|4.34| ), which is regular at |x — y| —>■ 0, reflects all the boundary 



conditions. One may see then that the UV singularities which trouble Q^'^ may be treated 
in the same way as at the T = 0. Incidentally, the former is an extension of the well known 
fact that in order to renormalise a finite temperature QFT, it suffices to renormalise it at 
T = 0. 

Inserting the solution ( [4.36| ) into the expression for the energy-momentum tensor ( [4.31| ) 



we arrive at the thermal interaction pressure per particle (compare Section |3.4] , see also 



85,9 



rm = ^'^' T^m^jT) ^ Tml{T) ^ X^ f T' T' mr{T) , T^ ml{T) 



90 24 127r 8 V 144 247r IQtt 



r2 



4/^A T„ I'^riT) 



+ O [mt{T) In ^-j^jj, (4.38) 

where the renormalised coupling constant A^ comes from the T = renormalisation prescrip- 
tion: r|jj_^^^(s = 0) = —Xr/N (s is the usual Mandelstam variable). A direct calculation of 
r*^^) in the large- A^ limit were performed in Section p.2| with the result: 



A. = Ao - ^AoA.^^^r [2 - D/2] [mlf"-' . (4.39) 

To regularise the theory we have used the usual MS scheme with the dimensional regular- 
isation (D 7^ 4). The high temperature expansion of the pressure ( |4.38|) to all orders can 



be found in Section ^^ where the calculations, however, were approached from a different 
standpoint. 
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4.4.2 Off— equilibrium I 

The next question to be addressed is how the above calculations are modified in the non- 
equilibrium case. To see that let us choose the following constraint 

g{k) = {n{k))k = {Him,^ . (4.40) 

Here t, is arbitrary, and function g{]<.) is specified by theory /experiment. The 7Y(k) is a 
quadratic operator fulfilling the condition {H)\n^oo = ('^)- As g(k) is finite, both H and 
7i must be renormalised (i.e. we must subtract the zero temperature part). Obviously 

n{k)=ujka^{k)a{k), 
with ujk = Vk'^ + M"^ and 



M' = ml+'-X,l0^^GMk'). (4.41) 

In the large-iV limit the corresponding density matrix ([4.19 ) reads 



with /5(k)/2(27r)^co'fc being the Lagrange multiplier to be determined. Using Eq.( [4.21| ), we 
find that 

sW-p^JfBl^. (4.43) 



where V denotes the volume of the system. Clearly, expression ( |4.43| ) can be interpreted as 



the density of energy per mode. The fact that /?(k) is not constant indicates that different 
modes have different 'temperatures'. 
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The Kadanoff-Bayin equations coincide in this case with those in ( [4.34| ) provided 



rririT) -^ A4. The boundary condition can be worked out simply using the prescription 
(|42^ ). This gives 

G+_(A;) = e-'^(^)^°G_+(fc) . (4.44) 

The fundamental solution of the Kadanoff-Baym equation is 



G^,ik) = ^,}2\^^l^^^ - 2^^Sik' - M')Uik) 



eai3 = e(cr3)a/3 ■ (4.45) 

Let us mention one more point. The boundary condition ( |4.44| ) is not by itself sufficient to 
determine /^/j's. (This fact is often overlooked even in the equilibrium QFT.) It is actually 
necessary to substitute this condition with an additional condition, namely the condition 
of causality. The causality condition, i.e. vanishing of the commutator [$(x), $(?/)] for 
(x — yY < 0, importantly restricts the class of possible /a/3's. To see this, let us look 
at the Pauli- Jordan function ( [$(x), $(?/)]). The latter is the homogeneous solution of 
the Kadanoff-Baym equation with the initial-time value data: {[^{x),^{y)])\xo=yo = 
(i.e. causality condition) and do{ [$(2;), $(?/)]) |xo=j/o — "(^^(x — y). Thus the Pauli- Jordan 
function is uniquely determined and its Fourier transform reads 

F.T.(([<l>(x), <^{y)])){k) = -i27c6{k^ - M^) e{ko) . (4.46) 

Relation ( 4.46|) immediately implies that 



U_(k)=9{-ko) + f{k) 
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/_+(A;) = e{ko) + J{k) 
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(4.47) 



with / being, so far arbitrary, and for both /+_ and / ^ identical, function. The / is then 

fixed via the generahsed KMS condition ( [4.44|) . Similarly, the causality condition specifies 
that G++(A;) - G__(A;) = PP{ l/{k'^ - M^)} (the symbol 'PP' denotes the principal part). 
By inspection of the definition of Gq^ one may easily realise that 

G++ + G__ = G+_ + G_+ 
G+_(fc) = -(G+_(-fc))* 
G++(A;) = -(G__(-A;))* . (4.48) 

From these relations follows that /++ = / = / and f{k) = {f{—k))*. The / is the same 

as in (14.471 ). 

Since the divergences in Gq,/3 come from the first term in ( |4.45| ) (i.e. from the particular 
solution), we can shift the corresponding (zero temperature) poles at Z) = 4 to the bare 
mass. In this case we can write 

M^ = ml + 6m^ , (4.49) 

where rrir is the renormalised mass in the vacuum and 6m is the mass shift due to an 
interaction with the non-equilibrium medium. Inserting the '++' components of ( [4.45| ) 
into (14.41|) , we obtain 



M^ = ml + -Ao 



M 



D-A 



(4vr) 



-F [1 - D/2] + NiM"' 



(4.50) 



where 



N{M^ 



d'^k 



(2vr) 



D 



27rS{k'-M')U4k)- 



(4.51) 
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In an equilibrium system /++ would be the Bose-Einstein distribution fs- Note that 
because ( [4.5 1| ) corresponds to a homogeneous solution of the Kadanoff-Baym equation at 
\x — y\ = it is automatically finite. Thus for a translationally invariant medium (both 
equilibrium and non-equilibrium) /++ must act as a regulator in the UV region. 

Let us consider the expression for the expectation value of O^*^ given in ([4.31|) . It is a 
matter of a few lines to show that 



(4.52) 



with G being the T = causal Green function. Using ( [4.39| ), ( |4.41| ) and ( |4.45| ) one may 
directly check that Eq.( |4.52| ) does not contain UV singularities when the limit D ^ 4 is 
taken. This verifies our introductory remark, that (|4.37|) is finite even in non-equilibrium 



context. From the generalised KMS condition ( [4.44| ) and from ( [4.45| ) we obtain that /++ is 



/++(*) -^Sioi^- ("-53) 

So far the results obtained were completely general and valid for any translationally invari- 
ant system. Let us now consider a system in which 5f(k) = -^^exp {—ujk/a). As we shall 
see, this condition corresponds to systems where the lowest frequency modes depart from 
strict equilibrium, whilst the high energy ones obey standard Bose-Einstein statistics. This 



behaviour is typical of many off-equilibrium systems, eg. ionised atmosphere [108], laser 



stimulated plasma [|109|| , hot fusion ||110|| , etc. The a is usually referred to as the ionisation 
half-width of the energetical spectrum. From ( [4.43| ) we can determine /3(k) as a function 
of the physical parameter a: 
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/5(k) = - + — Y ^^ e^""^/'^ . (4.54) 

a ujk ^ n 

n=l 

To proceed, some remarks on the interpretation of /3(k) are necessary. Firstly, Eq.( |4.54j ) 
implies that for a sufficiently large uok {oJk ^ cr) the function /3(k) is approximately constant, 
and equals 1/a. Thus at high energies the distribution /++ approaches the Bose-Einstein 
distribution with the global temperature T ^ a. In other words, only the soft modes were 
sensitive to processes which created the non-equilibrium situation. The interpretation of 
a as an equilibrium temperature, however, fails whenever a ~ 0;^. Instead of a one may 
alternatively work with the expectation value of /3(k), i.e. 



^ /ci^k/3(k)e-^/- 



/#ke 
1 E-.^A-.((n+l)A</.) 



where Kn is the Bessel function of imaginary argument of order n. Because the system is for 
the significant part of the energetical spectrum in thermal equilibrium, l/(/3) approximates 
the corresponding equilibrium temperature to high precision. An interesting feature of 
( 4.55| ) is that it is quite insensitive to the actual value of Ai. Dependence of (/?) on A4 for 



a sample value a = lOOMeV is depicted in FIG. [4.1| . An important observation is that the 

(-1)"+^ _ 



asymptotic behaviour of {P) at cr ^ 00 goes like (/5) ~ a/ a, where a = 1 + |X]^i, 



In2 + 7r724?« 1.1. 
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Figure 4.1: A plot of (/3) vs. M at a = 100 MeV. 



Using (|4.37| ) and ( |4.52|) we get for the pressure per particle 



V{a) 



27r2 



M'a'K2{M/a) 
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(4.57) 



Note that Vq conies from the UV divergent parts of (|4.52| ). Whilst the separate contributions 
are UV divergent, they cancel between themselves leaving behind the finite Vq. As we 
already emphasised, the divergences come from the particular solutions of the Kadanoff- 
Baym equations. Because the former do not directly reflect the boundary conditions, the 
form of Vq must be identical for any translationally invariant medium. The non-trivial 
information about the non-equilibrium pressure is then encoded in terms in the brace [. . .] 
in (^36|) . 

Let us now perform the "high-temperature" expansion of the pressure ( [4.56|) . As a 
"temperature" parameter we may chose a. In this case we have 
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(4.58) 



where 
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M 



2k+2 



The '?/'(...) is Euler's psi function. For a sufficiently large a the leading behaviour of 5m? 
may be easily evaluated. To do this, let us ffist assemble ( [4. 391) and ( [4.49|) -( ^.51j ) together. 
This gives us the (renormalised) gap equation 



with 
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Setting X = 5m? /m^. and s = a /m.^. we obtain the following transcendental equation for x 
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(4.60) 
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The graphical representation of Eg . (^4. 601) is depicted in FIG.4.2. 



b) 




Figure 4.2: A plot of the Eq.( ^.60| ): a) the general shape, b) a small x behaviour. 



From this one may read off that for large x {x > 50) there exists a critical value of A,, above 
which the gap equation has no solution. (The plateau is actually bent downward with a 
very gentle slope.) FIG. |4.2| b clearly shows that if A^ ^ 1/s < 1 then x <^ 1. Using the 
asymptotic behaviour of Ki[z] for z ^ {Ki ~ {z)^^) we arrive at more precise estimate 
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of Xr for which x <^ 1; namely A^ ~ 1/s^ = m'^/a'^. This estimate is very helpful for the 
asymptotic expansion of 5m'^. For a sufficiently high a we may write 
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Inserting this back into (|4.56| ) we get 
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(4.62) 
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It is interesting to note that l/vr^ ?^ 7r^/97 which is almost the Stefan-Boltzmann constant. 
This once more vindicates our interpretation of a as a "temperature" . A plot of the pressure 



as a function of o is depicted in FIG. 4.3 
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Figure 4.3: A plot of pressure as a function of T, a for tti^ = lOOMeV. The gray line corresponds 
to equilibrium pressure, the black line corresponds to pressure ( [4.62 ). 
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Figure 4.4: A plot showing the difference of equilibrium and non-equilibrium pressures for rrir 
lOOMeV. 



It is due to the low frequence modes contribution to the pressure that P{(t) < P{T). This 
is a direct result of our choice of 5'(k), namely that a cannot be interpreted as temperature 
for the low frequence modes (i.e. tUfe < a). The smaller a is the less important contribution 
from non-equilibrium soft modes and so the smaller difference between both pressures. 



The result ( 4.62 ) can be alternatively rewritten in terms of (/5) . Using, for instance, 



the Pade approximation | 111 | for {(3), we arrive at 



V{{(3)) = 0.0681122 (/5)-^ - 0.0415368 (/3)"2A^^ 

+ A^ (-0.000647 (/?)-^ + 0.0000164 (/?)-2A^2^ 

+ 0{MHn{M{f3)y,Xl) . (4.63) 

The coefficient 0.0681122 ^ 7r^/145 is 1.6 times smaller than the required value for the 
Stefan-Boltzmann constant, so the parameter (/?) is a slightly worse approximation of the 
equilibrium temperature than a. In practice, however, it is a matter of taste and/or a 
particular context which of the above two descriptions is more useful. 
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4.4.3 Off-equilibrium II 

As was noted above, it is the specific form of the constraint g{\i.) which prescribes the 
behaviour of /3(k). Let us now turn our attention to systems which depart 'shghtly' from 



equihbrium, i.e. when g(k) in ( [4.40|) deviates a httle from the equihbrium density of energy 



per mode. In this case the constraint (|4.4CI| ) reads 



^(k)=(^^ft;i^ + ^^(k); Sg{k)«g{k), (4.64) 

with Pq = 1/Tq being an inverse of the equihbrium temerature. As a special example of 
(|4.64| ) we choose 



«"-'-(2^ e^-JHk)-l ^ a(k,A)«l. (4.65) 

The inverse mode "temperature" /?(k) is then /5o — ln(a(k))/u;fc. So ln(a) measures (in units 
of Uk) the deviation of the mode temperature from the equilibrium one. The particular 
value of a{k, Pq) depends on the actual way in which the system is prepared. To avoid 
unnecessary technical complications, we select a(k, Pq) to be a momentum-space constant 
(generalisation is, however, straightforward). This choice represents the change in the 
mode temperature which is now inversely proportional to Uk', the deviation is bigger for 
soft modes and is rapidly suppressed for higher modes. Obviously, Tq becomes the global 
temperature if cu^ ^ ln(a;). The generalised KMS conditions ( [4.44[ ) together with solutions 
(|4.45|) determine /++ as 



/-W = P5^- ("-e") 



Eq.( [4.66|) is a reminiscent of the, so called, Jiittner distribution^ with fugacity a [p.l2| , |113 



^It should be mentioned that this similarity is rather superficial. The Jiitner distribution describes 
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Now, using ( ^.37| ) and ( [4.52| ) we get for the pressure per particle 



V{To)+Vo 



27r2 



^' n=l n=l 



where Sm^ satisfies the gap equation 






(4.67) 



(4.68) 



If we set, as before, x = 5m? /m^ and s = l/Pof^r we get the transcendental equation for x 



A;' = ^^ j(l + x)ln(l + x)-x + 8(l + x) f 



dz \Ji 
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oz Vx+T/s 



a 



(4.69) 



The corresponding numerical analysis of ( [4.69p reveals that for x ^ 1 , 1/s ^1. So at 
X -C 1, Ar ~ 5m? /Tq. This estimate is important for the asymptotic expansion of 5m?. 
However, in order to carry out the asymptotic expansion of ([4.68|) (and consequently (f4.67|) ) 
we need to cope first with the sum Yl^=i Oi"'Kk{nJV{.[3Q)/n^ {k = 1, 2) (also called Braden's 
function). Expansion of Kk{. . .) yields a double series which is very slowly convergent, 
and so it does not allow one to easily grasp the leading behaviour in Tq. In this case it is 
useful to resume ( [4.67|) -( p:.68|) by virtue of the Mellin summation technique |^ . (It is well 
known ^,|7^,^ that at equilibrium this resummation leads to a rapid convergence for 
high temperatures.) As a result, for a sufficiently large Tq we get 



5m 



XrT? XrToM 



24 



Att 



.2U/2 



systems which are in thermal but not chemical equilibrium. (As we do not have a chemical potential, 
chemical equilibrium is ill defined concept.) The fugacity then parametrises the deviation from chemical 
equilibrium. 
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— r j 1 — In j — ;- J J — (1 — r^)^" arcsin(r) 

+ O(lnTo), (4.70) 

where we have set r = \n{a)TQ/A4. The corresponding expansion of the pressure ( [4.67| ) 
reads 

2 / . i \ 

— - (l — r'^) ( — r'^ + 3 — 3 (l — r^) ^ ) arcsin(r) 

+ ^^ + 0(lnTo). (4.71) 



where ({3) ~ 1.202. Note that for a = 1 we regain the equihbrium expansion ( 4.38| ). The 
corresponding plot of "P as a function of Tq and a is depicted in FIG. ^4.5| . 
In passing it may be mentioned that the expansion ( [4.71 ) is mathematically justifiable only 
for a ^ ev^^^^ ^ 1. 

4.5 Conclusions 

In order to get a workable recipe for non-equilibrium pressure Jaynes-Gibbs principle of 
maximal entropy, the Dyson-Schwinger equations, and the hydrostatic pressure. The basic 
steps are as follows. 

To find quantitative results for pressure one needs to know the explicit form of the 
Green's functions involved. These may be find if we solve the Dyson-Schwinger equations. 
The corresponding solutions are unique provided we specify the density matrix p (the 
construction of p is one of the cornerstones of our approach, and we tackled this problem 
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To[GeV] 
Figure 4.5: Behaviour of the pressure ( [4.67| ) as a function of a and Tq at rrir = lOOMeV. 

using the Jaynes-Gibbs principle of maximal entropy). With p at our disposal we showed 
how to formulate the generalised Kubo-Martin-Schwinger (KMS) boundary conditions. 

To show how the outlined method works we have illustrated the whole procedure on 
an exactly solvable model, namely 0{N) $^ theory in the large-iV limit. This model is 
sufficiently simple yet complex enough to serve as an illustration of basic characteristics 
of the presented method in contrast to other ones in use. In order to find the constraint 
conditions we have considered two gedanken experiments in which the system in question 
was prepared in such a way that only low frequence modes departed from the strict equilib- 
rium behaviour. Such processes can be found, for example, in ionised atmosphere, in laser 
stimulated plasma or in hot fusion. In both alluded gedanken experiments we were able 
to work out the hydrostatic pressure exactly. Furthermore, after identifying a "tempera- 
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ture" parameter (virtually temperature of high modes) we carried out the corresponding 
high-temperature expansions. 

As it was discussed, one of the main advantages of the Jaynes-Gibbs construction is 
that one starts with the (physical) constraints (i.e. parameters which are really controlled 
and measured in experiments). These constraints directly determine the density matrix and 
hence the generalised KMS conditions for the Dyson-Schwinger equations. This contrasts 
with the usual approaches where the density matrix is treated as the primary object. In 
these cases it is necessary to solve the von Neumann-Liouville equation. This is usually 
circumvent using either variational methods ||10CI| , |103|| with several trial p's or reformulating 



the problem in terms of the quantum transport equations for Wigner's functions [^. It 
is, however, well known that the inclusion of constraints into transport equations is a 
very delicate and rather complicated task (the same is basically true about the variational 
methods) Pi[TTl. 



Chapter 5 

Summary and outlooks 



In this dissertation we have explored various aspects of equihbrium and non-equihbrium 
quantum field theory. 

In Chapter we dealt with certain aspects of infrared effects in finite-temperature QFT. 
It is well known p^,|36| that these effects are far more subtle at finite temperature than 



at T = 0. The latter is ascribed to the fact that there is no finite-temperature Bloch- 
Nordsieck mechanism [p3i|115|| at our disposal. Intuitively, however, one might expect that 
if thermal QFT is a well posed theory, no infrared divergences should be present. This point 
was extensively studied on particles that decay and scatter inside a plasma (or heat bath). 
As a byproduct we proposed an easy mathematical formalism allowing one to calculate 
the plasma particle number spectrum formed when a particle (s) decays (scatter) within 
the plasma. This formalism which is based on the largest-time equation and the Dyson- 
Schwinger equations, is embodied in a one-line modification to the corresponding thermal 
cut diagrams. Such diagrams arise naturally when the imaginary part of a Green function 
is computed, and the whole calculations are immensely simplified by the observation that 
the thermal cut diagrams are virtually the same as the Kobes-Semenoff diagrams in the 
Keldysh formalism. 
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The vital part of Chapter |^ was dedicated to the study of the modifications which must 
be made in various quantum field contexts (scalar fields, spinor fields and gauge fields in 
temporal gauge). It was found that the modified propagators have an easy interpretation 
in terms of emission, absorption and fluctuation of the plasma particles. To demonstrate 
how the method works, we used various heat baths (photon, and electron-photon) in which 
two uncharged scalar particles scatter into a pair of charged particles. We calculated the 
leading u behaviour for the resulting changes in the plasma particles number spectra. We 
also found that the energy density udN/du is finite as u tends to zero. 

In Chapter |^ we have turned our attention to pressure. Pressure is undoubtedly one of 
the important parameters characterising quantum media at finite temperature. An exten- 
sion of the pressure calculations to off-equilibrium systems should enhance our predicative 
ability in such areas as (realistic) phase transitions, early universe cosmology or hot fu- 
sion dynamics. Unfortunately, the standard approaches based on the partition function 
do not allow such an extension because there is no such thing as the grand-canonical po- 
tential away from equilibrium. We considered an alternative approach based on the, so 
called, hydrostatic pressure. As a warming up exercise we started with thermal equilib- 
rium. The whole procedure was then illustrated on a toy model (A$'^ theory with 0{N) 
internal symmetry in the large- A^ limit). The result obtained matches that found in [^ 
for the thermodynamic pressure. (Whilst these two pressures agree in equilibrium , there 
is no thermodynamic pressure out of equilibrium.) 

We took advantage to probe various mathematical techniques (composite operators 
renormalisation, Dyson-Schwinger diagrammatic equations, finite temperature coupling 
renormalisation, Mellin transform) which are indispensable in calculations of the hydro- 
static pressure (most of them also being applicable to non-equilibrium as well). One of the 
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key results obtained was the prescription for the renormahsed pressure (Eqs.( p736D , ( [4.4q ) 
and ( [4.37| )). This was achieved by means of the Zimmerman forest formula. 

The model in question has the undeniable merit of being exactly solvable. In particular, 
one could find the pressure by means of a non-perturbative treatment. This is because 
of the fact that the large-A^ limit eliminates 'nasty' classes of diagrams in the Dyson- 
Schwinger expansions. The surviving class of diagrams (superdaisy diagrams) could be 
exactly resumed because the (thermal) proper self-energy S as well as the renormahsed 
coupling constant A,, were momentum independent. The resulting expression for the pres- 
sure obtained was then in a suitable form which allowed us to take advantage of the Mellin 
transform technique, and we were able to evaluate the pressure in D = 4 (both for massive 
and massless fields) to all orders in the high-temperature expansion. 

As we have already mentioned, there are various motivations to be interested in the non- 
equilibrium pressure. One of them is based on the belief that the pressure (as an easily 
measurable parameter) should exhibit a discontinuity in its derivative(s) when a local phase 
transition occurs. This has been observed in solid state physics and fluid dynamics and it 
may play a central role in, for example, detection of the quark-gluon plasma or in various 
baryogenesis scenarios. 

In Chapter ^ we approached the problem of pressure in a non-equilibrium quantum 
plasma by extending the notion of the hydrostatic pressure to systems out of equilibrium. 
On the equilibrium level the hydrostatic pressure; i.e. p = — 1(0* j) (here Q'^'^ is the energy- 
momentum tensor of a system in question), has a well defined microscopic picture, which can 
be carried over to an off-equilibrium as we have showed in Chapter ^. Unlike in equilibrium, 
the non-equilibrium expectation value is sensitive to the particular choice of the energy- 
momentum tensor (see Chapter 13), though fortunately in translationally invariant media 
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this sensitivity is not relevant. Our strategy for the calculation of pressure was simple. To 
find quantitative results for pressure one needs to know the explicit form of the Green's 
functions involved. These may be found if we solve the corresponding Dyson-Schwinger 
equations. The solution is unique provided we specify the density matrix p, and hence the 
generalised KMS conditions. In order to construct the density matrix we invoke the Jaynes- 
Gibbs principle of maximal entropy. The Jaynes-Gibbs principle is basically the Bayesian 
inference about the most probable density matrix p. This is based on the maximalisation of 
the Shannon (or information) entropy, subject to the prior knowledge which one has about 
the system (usually specified at some initial time ti). 

To show how the outlined method works we have illustrated the whole method on our 
favoured 0{N) ^'^ theory in the large- A^ limit. In order to specify the initial-time con- 
straints we have considered two gedanken experiments in which the system in question was 
prepared in such a way that only low frequency modes departed from the strict equilib- 
rium behaviour. Such processes can be found, for example, in ionised atmosphere, in laser 
stimulated plasma or in hot fusion heated up by ultra-sound waves. Furthermore, after 
identifying a "temperature parameter" (virtually the temperature of the high modes) we 
carried out the corresponding high-temperature expansions. The leading high-temperature 
coefficients coincided to a good approximation with the Stefan-Boltzmann constant which is 
precisely what one would expect: at high "temperature" particles do no feel non-equilibrium 
"background" and consequently their distribution approaches the equilibrium one. 

The next logical step in our investigation would be to calculate the pressure in a non- 
equilibrium, non-translationally invariant medium. In contrast to the previous situation, 
one is confronted now with two conceptual difficulties. As we have shown in Section [4.4.2 



m 



order to get the renormalised pressure one has to consistently use the renormalised energy- 
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momentum tensor (0'^'" is a composite operator). Unfortunately there are infinitely many 
renormalised energy-momentum tensors (generated from each other via Pauli's transfor- 
mation), and these give rise to different pressures. The choice of the "correct" Q^'^ is then 
crucial for the next investigation. (As we discussed in Sections ^^ and [4.4| , in translation- 
ally invariant media pressure does not depend on a particular form of B^'^.) Our future 
strategy is based on the observation that because the pressure is an observable, one should 
obviously work with an energy-momentum tensor which is an observable as well (or the 
corresponding Hermitian operator). We propose to use the Belinfante energy-momentum 
tensor. The fact that the physical significance of the Belinfante tensor (both in classical and 
quantum systems) is greater than that of other energy-momentum tensors was pointed out 
successively by Belinfante ||116|| , Rosenfeld [ p.l71| and Jackiw ||6^. In scalar QFT the Belin- 



fante energy-momentum tensor coincides with the, so called, improved energy-momentum 
tensor [Q. The generic form of the improved energy-momentum tensor for the 0{N) $^ 
theory was derived in Section |^ Eg. ( p. 331 ). Because the corresponding (3 function is ex- 
actly solvable in the-large A^ limit |jl3| we expect that constant c{Xr', D) is exactly solvable 
as well, and consequently the improved energy-momentum tensor can be written in closed 
form. 

The second problem to be solved is connected with the fact that the corresponding 
Dyson-Schwinger equations are far more involved than in the translationally invariant case. 
Even, for the 0{N) $^ theory in the large- A^ limit the Kadanoff-Baym are not any more 
hyperbolic equations. This fact, among others, complicates the renormalisation program 
(one cannot use the momentum space renormalisation). For this purpose we intend to resort 
to differential renormalisation |p.l8|] . Solving the renormalised Kadanoff-Baym equations 



[at least numerically) is the next step in our study. 



Appendix A 

Finite— temperature Dyson— Schwinger 
equations 

A.l Functional formalism 



Eq.( p2.18D gives us an alternative definition of Wick's tlieoreni in terms of the "functional 



derivation" -gjjrpr. We refer to Eq.( ^.18| ) as the Dyson-Schwinger equation because the 
classical T = Dyson-Schwinger equations are implied by it (see Appendix |A.3|) . Let us 
first show that ( p.l8| ) is consistent with Wick's theorem (|2.16| )- (|2.17|) . To be specific, let 
us consider an ensemble of non-interacting particles in thermodynamical equilibrium. In 
order to keep the work transparent, we shall suppress all the internal indices. There is no 
difficulty whatsoever in reintroducing the necessary details. Let us first realize that for any 



(well behaved) functional the following Taylor's expansion holds [^ 



^/<'— / 



^[i^] = ^ dxi... dxna"{xi . . . Xn)ip{xi) . . . 1p{Xn) , (A.l) 

The same is true if ip is an operator instead. In the latter case the «""(...) are not gen- 
erally symmetric in the x's [[ When Fermi fields are involved, we might, for the sake of 
^If X = X[tp, dip], the q;" may also contain derivations working on the various fields. 
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compactness, include in the argument of ip the space-time coordinate, the Dirac index, and 
a discrete index which distinguishes ipa from ip^. In the latter case J dx ^ ^ J dx, where 
summation runs over the discrete indices. With this convention, the expansion (|A.1|) holds 



even for the Fermi fields. An extension of (|A.1D to the case where different fields are present 



is natural. Particularly important is the case when ■?/; is a field in the interaction picture. 



using Wick's theorem and decomposition ( |A. 1| ) one can then write 



{G[^]^P{x)F[^]) = 

(p \ n n In 

/ dx ) «"(xi ...Xn) 5^(±1)"-' {i'{xi)i^{x)) I X{i^{xk)F 

/ dy r{yi ■■■ym) $^(±1)'-' mx)i^{xi)) ( G[i,] n i'iy, 



I \k^l 

m I m 

k') ) ■ (A.2) 



with [J dx^ = J dxi . . . J dxn- The ' — ' stands for fermions and '+' for bosons. On the 
other hand, using the formal prescriptions (p.l9|) and ( p.21|) for ^J,, one can read 



dz{^{x)tlj{z)) { G"'^'^ 



6ip{z 

/ dy riyi ■■■ym) dz{^ij{x)^{z)) ^(±i)'-i(5(z - m) I G^ n Hy 



k'f^l 



dy] p-^iy,. . .yjJ2i±iy-\^ix)ij{m)) lG[,p]ll,p{y,,) ) . (A.3) 



Similar expression holds for J dz{ijj{x)^{z)) ( 's^jhr^bP] )• Putting latter two together we 
get precisely ( |A.2| ) . This confirms the validity of ( |2.18| ) . It is easy to persuade oneself that 



exactly the same sort of arguments leads to 
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mx)Fm = I dzmx)^{z)) ( ||^ ) (A.4) 

{r{^{x)F[^])) = I dz{r{^{x)^{z))) { T ( j^ I ) (A.5) 

(GMr(^(x)FM)) = / dz{T{^{x)^{z))) { cmr (^ I ) + 

+ ldzmz)^{x)) (^l^nFm) ) , (a.6) 

etc. 

with T being either the chronological or anti-chronological time ordering symbol. At 
this stage it is important to realize that from the definition of gj,-. directly follows that 
[wfe)' wt]=f — ( '^' holds for bosons and '+' for fermions). Indeed, 



s'^Fm ^^r f ,Y'\ .. i i 



S^{x)6^{y) 



^ ^ [ dx\ (a"(a;i . . . x . . . y . . .Xn) ± 



A 7 A ,. . ^'F\ 



± a"(xi...t/...i:...x„))(±ir+^' n i){x^) = T 

•m^i,j 



Similarly [^^; jf^]^ = 0. Analogously we might prove 

6^Fm 



6ilj{y)6ilj{x) ■ 

(A.7) 



6ip{x)6ip{y) 6'ip{x)6tp{y) ' 
and 



5tlj{x)5tlj{y) 6ij{x)6tP{y) ^^^ ^ ' 6tlj{x) #(y) 



(A.^ 
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The "p" is for bosons and 1 for fermions. With ( p.l8| ) and (A.4)-(A.6) one can easily 
construct more comphcated expectation values. For example, using ( p.l8| ) and (A. 4) we get 



(^(x)^(2/)F[^]) 

dzidz2 
2 



((^(x)^(;^i))(^(l/)^(^2)) + (-ir(^(x)^(^2))(^(l/)^(^i))) 



^^^^R^b^+^^^^^^^^^^^^f^^^ ^^-''^ 



Similarly, using (|2.18|) and (anti-)commutativity of the arrowed ^ , , . , we get 



&'il>{x) ' 



dzidz2 



X G[^] 



2 ((^(x)7A(;2i))(^(l/)^(^2)) + (-ir(^(x)^(;^2))(^(l/)^(^i))) 



(5^(zi)(5V^(2;2) 



+ / ^^((^(zi)^(x))(^(z2)^(2/)) + (-ir(^(^2)^(x))(^(^i)^(2/))) 



X ' 



2 



6i){zi)6i){z2) 



^^f^^) + <^^^)^^^)^<^[^™^- (^-^^^ 



We could proceed further having still higher powers of fields and variations. However, there 
is a quite interesting generalisation in case when we have (anti-)time ordered operators. 
Let us have F[^] = T(F [?/']), in this case 
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(Fm) = $^(/cixya'^(...)(r(n^(x.)) 

n=l ^ ' jj rnj^ij 

= /^V--.(rw<..)V<--.))> ( ,^,(yi- },.) ) + (f l°l> ■ (A.12) 

where F[ip] differs from F[ip] in the replacement a"(. . .) -^ '"' (n starts from 1 !). In 
comparison with (A.4)-(A.ll), the a°(. . .) (i.e. the pure T = contribution) does matter 
here. Note that a^{. . .) generally involves non-heat-bath fields with corresponding space- 
time integrations. Similar extension is true if F[iIj] = Tc{F[ip]), where Tq is the time path 
ordering symbol. In that case 

(^M) = 5^ (/ c^x) "«"(... )(Tc(nV^(x,))) 

n ^-^^ ' p=l 

I'^'FMA \ 



+ (i^[0]), (A.13) 

with| j(jdx = J^dt Jydx. and -£jpr = 5c{x — y) . Wick's theorem for the Tc'-oriented 
product of fields has an obvious form 



{Tciijix^) . . . ^(X2„))) = 5^£p(rc(^(x,)^(x,))) (Tc( n H^k))) . (A.14) 



^A contour J-function 6c{x — y) is defined as /^ dz6c{z — z )f{z) = f{z ), see [ p] , |ll9[ . 
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This can be directly derived from Wick's theorem ( |2.17| ), reahzing that 



Tc{ip{xi)...tlj{xm)) = ^spedtp^, . . . ,tpjilj{xp,) . . .tpixpj , (A.15) 

p 

where P refers to the permutation of the indices and 9c{ti, . . . ,tm) being a contour step 



function 731 defined as 



, . ,1 (ti,...,tm are 7^-oriented along C) /A1«^ 

"0 (otherwise) 



So for example ( |A.6|) may be written as 



{Tc{^P{x)F[m = l^dz{Tc{^{x)ij{z))) (Tc I ^^ | ) . (A.17) 



Particularly important is the Keldysh-Schwinger path [^ , [73|, |120[| (see FIG . 13.31) . In the 
latter case 

{Fm) = f dz,dz,{r{^{z,mz,))) ( ^.ffj^l . 

+ I dZidz2{T{ilj{zi)^{z2))) ' 



C2 \6^/j{z2)Sip{zi) 

5^Fm 



+ 2 dzi dz2{ip{zi)^lj{z2)) . 

Jc2 Jc^ \6tp{z2)6^{zi) 

+ im)- (A.18) 

Application to the product G['ilj]F[ip] with Flip] = TdFlip]) and G[^] = Tc,iG[ip]) is 
straightforward and reads 
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{G[^]F[ij]) = I dz,dz,{T{i^{z,mz,))) ( j;i;^^f[ 



+ I dzidz2{r{^{z,)^{z2))) { G[ij] ^'^ 



6i/j{z2)6^{zi) 



+ 2 I dz,dz,m^M'{z,)} ' ^V']SSFM 



6^Ij{zi) 6^1j{z2) 
+ (G[0]F[0]), (A.19) 

where the overhning indicates that we work with ^'"jf ^"'' instead of «""(. . .)/^™(- • ■)^ "we 



S4>{x) 



have also abbreviated f^, dz to J dz J^, dz. We should also emphasise that ^V-J^l used in 
( |A.19| ) is 6{x — y) rather than 5c{x — y). 



In Eg. ( p. 361) it has been used the inverted version of (|A.19|) , namely 






5'ilj{z2)5i){zi) 
+ / dz,dz2{ij{z,)i;{z2) { i^^i^) , (A.20) 



\6ij{z,) 6ij{z2) / ' 



Here {G[ip]F[ilj])' has the coefficients a"(. . .)/?"(• • .)^^ instead of a"(. . .)/5™(. • •)• Note, 
that the a°(. . .)/5°(- • •) does not contributes and thus we do not have any pure T = 
contributions. Eq.( |A.2(3|) has a natural interpretation. Whilst the LHS tells us, that from 
each thermal diagram (constructed out of {G['tp]F[tp])) with ^^^^^ internal heat-bath particle 
lines we must take n + m identical copies, the RHS says, that this is virtually because we 
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sum over all possible distributions of one heat-bath particle line inside of the given diagram. 
The pictorial expression of ( |A.2CI| ) is depicted in FIG. [A.1| . 



X (number of lines) = 2_^ 





Figure A.l: Diagrammatic equivalent of Eq.( |A.2d ). The cut separates areas constructed out of 

F[V'] and G[i)]. 



A. 2 Graphical formalism 



(This section is partially based on refs. [ff2| , |121|| .) 



In Chapter ^ we heavily used a graphical representation of the Dyson-Schwinger equations. 
In this section we provide a short derivation of such a representation and in the section to 
follow we perform a comparison with the functional Dyson-Schwinger equations of Section 
H and Appendix [A.l . 



The Dyson-Schwinger equations were originally constructed ||122| , |123|] with the motiva- 
tion that they could provide some information about the complete Green's functions outside 
the scope of perturbative theory. The basic philosophy is to directly use the equations of 
motion in order to construct the hierarchy of (integral) equations for full Green's functions. 
The usefulness of these equations is usually confined to the cases where some approximation 
(truncation) is available in order to bring them into manageable form. In this appendix we 
aim to derive the finite-temperature Dyson-Schwinger equations using the more intuitive 
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path-integral formalism. The alternative derivation in the non-equilibrium context (based 
purely on operatorial approach) is discussed in Section [4.1.1| . 

Let us start with the (T ^ 0) generating functional of Green's functions Z[J] (the 
partition functions with an external source J) 



Z[T]= / V(f)exp{i{S[(f);T]+ / rf^x J(x)0(x))} 

S[4>-T] = f d^xC{x). 
Jc 



(A.21) 



Here J^ d^x = j^ dxQ jy d^x with the subscript C indicating that the time runs along some 
contour in the complex plane. In the real-time formalism, which we adopt throughout, the 
most natural version is the so called Keldysh-Schwinger one [^,[3^, which is represented by 
the contour in FIG. |3.3| . Within the path-integral the c-number fields are further restricted 
by the periodic boundary condition (KMS condition) [Bl, RB],pB| which for bosonic fields 
reads: 

(j){U-t(3,^)=(j){ti,^). (A.22) 



The zero temperature generating functional may be recovered from ( |A.21| ) if we integrate 
over the close-time path (no vertical parts) and omit the KMS condition ( |A.22| ). Now, the 
LHS of ( [A. 211 ) is independent of 0, thus particularly it is invariant under infinitesimal point 
transformation 



{x) -^ (t>{x) +ef{x) 



(A.23) 



where / is an arbitrary (0 independent) function which fulfils the periodic boundary con- 
ditions 

/(t,-^Ax) = /(t„x). (A.24) 
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Prom ( |A.23| ) is obvious that the corresponding functional Jacobian is 1, i.e., V(f) = V(j)'. 
This immediately implies that 



Z[J] = I V(j)'exp{i{S[(j)' -ef]+ / d'^x J{x)(j)'{x) - e d^xJ{x)f{x)) 

Jc Jc 

-IS lV<p'\j d^x f ^(x) + J{x)] /(x)l exp(«^[0'] + I [ J0') + 0{e^) 



c 



\( 



c 



+ Z[J], 



.4. f5S, 



V<P{ I d^x{—ix) + Jix)\fix)^expitS[(j)]+tJ J<P) 



As ( [A.25| ) is true for any /(...) fulfilling the condition ([A.24|) one may write 



(A.25) 








/ ^s[M 



[x) + J{x) ) fix) 



d^x \ 



5 



i6J{x) 



+ J{x) ]Z[J]. 



(A.26) 



Employing the commutation relation: —ijjZ = Z{ 
more appropriate form, namely 



ijj), we may recast ( |A.26| ) into 



5S 
— J[x) = -— 



[x] — I 



5J{x) 



I^^ 



{x)+t f dhWix^z)-^ 
Jc H{z). 



I. 



So, for instance, for A $^ theory we have 



(A.27) 



d'y d^w d^z WJx, y) W(x, w) T^^) P^\y^ w, z) W,^{z, x) 



[U + ml)<j)^{x) - ^ { {cp^{x)f + i?,<P^{x) D^ Jx, x) - Wjil{x, x, x) ] 



:n + m^)0,(x) + -^©(fiJ^,a;,x) 



(A.28) 
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(no sum over q;!)0 or, using identity; {Ox + ^o) ^f i^^ v) — ~{^'iY^^{^ ''V)-! ^^ ^^-Y invert 
the differential operator and write equivalently 



d'yJp{yWF{y.^)=r{x) + ^jdSnf{x,y)n%p{y,y,y). 



(A.29) 



Eq. ( [A.29| ) has the following graphical representation 




Z^/3=l, 



2 a 



6 /^/3=1,2 a 




where the hatched blob refers to the (full) 1-point amputated Green's function, the dotted 
blob refers to the (full) 3-point amputated Green's function, the cross denotes the source 
J, while the heavy dot without coordinate indicates that the vertex must be integrated over 
all possible positions. The corresponding graphical representation in terms of the connected 
Green's functions reads (c.f. Eq. (|A.28|) ) 



In ( A.2q ) we have implicitly used the identities; 



1 



5"Z 



Z SJ{xi) . . . SJ{Xn) 

SJ{xi) . . .SJ{x„) 
6(j){x) 



= (-l)"+izD(")(:ri,...,x„) 
(-!)"+! B-(")(xi,...,a;„) 



SJ{y) 

S(f>y{z) 



ix,y) 



5"T 



r(")(a;i...a;„) 



S<j){xi) . . .(j}{Xn) 

dytdylB'^^,ix,yi)r^'^'^''{yi,y2,z)ni^,iy2,y). 
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a 



E, 



a 



J 



E 



2 Z^f3 a 



^ 
P 



- -T 



6 ^/3 a 



4E 



6 /-^f3,^,5,K 



Here the double-hatched blob describe the 3-point vertex function and the hatched blobs 
refer, as before, to the connected (truncated) Green's functions. 

To see how to obtain the corresponding Dyson-Schwinger equations for the 2-point 
(connected) Green's functions, let us perform in Eq. (|A.27|) (or, for simplicity's sake, in 
Eq.( |A.2^ )) variation w.r.t. Jpiy). This directly gives us 



'D, 



+ ml) w^p{x, y) - ^ Dil^(x, X, x,y) + ^ My) »L1(^, ^> ^) 



icr3)afiS{x-y) 



or, inverting the differential operator, we may equivalently write 



X 
3 

Xoi 



B''^f,{x,y)=BF^p{x,y) - '-^ d'z]D);^{x,z)B['^^^{z,z,z,y) 



+ -^ I d'zB;^{x,z)]D)^^l{z,z,z)My)- (A.30) 



Eq.( [A.30| ) is graphically represented as follows 
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x_ 
a 




y 
P 



X 

a 



y 

p 



-y 

6 ^ 



7=1,2 



6 2^7=1,2 



X 


/^ 


^7 V 


a 


Y^---^ 


3^_y p 


X 


1 


j^\^f-^.j 


a 


.V 


\ A-'^ P 



This may be reformulated purely in terms of the connected Green's functions and the 
vertex functions, indeed, performing variation w.r.t. Jpiy) in first two lines in Eq.( |A.28|) , 
and inverting the differential operator as before we obtain 



J 1 Y^ x_ 

R 2 Z-/7 a 



p 2 .^7 a 



y 



a Y. 



^ 6 Z^7i « \ 



2 ':-^7i a 




Note that if the symmetry is unbroken (i.e. when ttiq > 0) then after setting J = the 
second, fourth and sixth term on the RHS will vanish. 

It is important to mention that similarly as for the connected Green's functions we may 
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write the Dyson-Schwinger equations for the vertex functions. Those may be obtained 
directly from Eq. ( |A.27] ) provided that one performs the corresponding number of variations 



w.r.t. (f). So, for instance, for the A$^ theory the 1-point vertex function F, 



(1), 



X) 



sr 



-Ja{x) is graphically represented (c.f. Eq. (|A.28|) ) as 



5(j)" (x) 



X 

a 



1 \]/ 




-\- t: X 

2 a 



_L 1 V Xl 

^ 6 Z^/3t5 -•- 




Here the slash stands for the inverse free thermal propagator. Particularly important ex- 
amples are the 2- and 4-point vertex functions V^"^^ and T'^^\ respectively, which are indis- 
pensable in the renormalisation prescription (see Chapter ^. In the case of the 0{N) A$^ 
in the large- A^ limit the corresponding graphical representations were explicitly calculated 
in Chapter ^ (see Sections p.3.1|^73l^ . 

A. 3 Comparison 

In order to find a link between the previous two approaches, let us start with the full n- 
point Green's function at T = 0. It is well known that in this case, the full Green's function 
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may be expressed by means of the Gell-Mann and Low formula [p^, |70| , pT5| , |121 



{0\T{^H{xi)...ijHiXn)m 



T (ij{xi) . . .ip{xn)e'^ '^"'^^'^''^ 



(0|r (e^/'^*^^^W)|0) 



(A.31) 



which may be equivalently rewritten in the Schwinger form |124| as 



(0|r(^^(xi) . ..^H{xn)m = (O Tc U{x^) . ..^{xn)e'fcd'-^i(-) 



(A.32) 



where ipn is a field operator in the Heisenberg picture, while ip and Cj, are operators 
in the interaction picture, the sub-index C denotes the Schwinger-Keldysh time contour. 
The Schwinger (but not the Gell-Mann and Low!) formula can be directly extended to 
finite temperature as we have mentioned in Chapters § and ^, the only difference is that 
(0| . . . |0) -^ (•••)• The following results will be equally valid for both T = and T 7^ 
provided we take into account corresponding expectation values. In order to be concise, 
we shall restrict ourselves to bosonic fields. The necessary extension to fermions and gauge 
fields is straightforward. 

For one field ipn in the presence of an external resource J{x) we have 



x) = {^h{x)) = (Tc ('^(a;)e*/c'^'-(^/(-)+^(-)^(-)) 



(A.33) 



Setting in (|XT7|) Fl^p] = e*/c'^*^'(^H^)+-^(^)'/'(^)), we get 



{^PHix))+ I d'zB,ix,z)(Tc'^^'^'^'^^'^y^ 



c 



6ij{z) 



H 



d^zBF{x,z)J{z), (A.34) 



c 
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where i]D)F{x,z) = {Tc{^{x)^{z))). For instance, for the A^"^ theory (£/ = — ^<l>^) we 
have 

($h(x)) ~^f d'zBp{x,z){<l>],{z)) = -fd''z3F{x,z)J{z), (A.35) 

or in the thermal-index notation; a, /3, 7, 5 = {1,2}, 

0"(x) + ^jd''z Wp\x, z) nfl{z, ^^z) = - j d^z Bf{x, z).Js{z) . (A.36) 



In ( |A.36|) we can recognise the Dyson-Schwinger equation ( |A.29| ) for the A^'^ theory from 



the previous section. 

To get further equations in an infinite tower of couplet integral equations we may perform 
successive variations in ( |A.35| ) with respect to Jpiy), or we may alternatively set in ( |A.17| ) 



F[ip\ = ijj[y) e'^fc'^ x{Ci(x)+j(x)ii:(x)) ^ jj^ both cases we obtain for A$'^ theory the following 
identity 

B'^^{x,y) = Bf{x,y) -'-^ j d'zWj{x,z) {a,)s, {Tc{^U^)^H{y))V' 

+ 1 j d^zWj{x,z)Js{z)(l)^{y) 
W'^^ix,y)=Bf{x,y)-^ J dhB'fix,z)Bf^/iz,z,z,y) 



^ 



^ I d'zBfix, z) Dg) (z, z, z) cp^iy) . (A.37) 



Note that ( |A.37|) precisely coincides with the Dyson-Schwinger equation ([A.30|) derived in 
the previous section. 
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Let us mention one more point in the connection with Eq.( |A.37| ). If, instead of setting 
F[i/j] = ip{y)e^-^c^ x{Ci{x)+j{x)i}{x)) ^ ^^^ ^^ have done in the previous case, we would set 
G[ij] = 1, F[^] = e^/c°'*^(^/(^0+^(^')^(^')) in ( [Alol ) or ( lAllI ), we would get somewhat different 
identity, namely 



x,y) = ]D)p{x,y) - d zid Z2 ©^"^(x, zi) (0-3)^5 

V /t f n^Icd'xiCJix) + Jixmx))U \\'\ . 



which in the case of A^^ theory Eq. ( [A.38|) may be explicitly written as 



-^2J d'zid'z, B'P{x,z,){a,U {Tc{^Uz,)^'h{z2))Y' (^3).. Of (^2,1/) 
- j d'z,d''z2 B'Pix, z^) J^{z,) Bf{z2, y) Js{z2) . (A.39) 



It is, however, not complicated to persuade oneself that ( |A.39| ) is equivalent to (^ 



Indeed, using ( |A.36| ) we get the desired equality ( |A.37| ). In ( [A. 381 ) we can recognise the 



Dyson equation with (Tc r ^'^ '^ 5^(i[]]'!p'(z2) )) ^^^^§ *^^ self-energy grj i'S{zi,Z2). 

Introducing a proper self-energy S, which is a self-energy that cannot be separated into 
two pieces by cutting a single line (i.e. in the matrix form S = S + SD^rS + . . .) the series 
( |A.38| ) can be summed formally to yield 



B''''{x,y) = Bf{x,y)+fd^z,d'z2Bf{x,z,){-zi:s,{z,,Z2)W^{z2,y) 
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which is the usual form of the (thermal) Dyson equation (c.f. equation ( |3.43| )). So the second 
in the hierarchy of the Dyson-Schwinger equations is nothing but the Dyson equation for 
the full two-point Green's function. Analogous procedure can be repeated for the Dyson- 
Schwinger equations involving the higher-point Green functions. As a result we can see 
then that the usual Dyson-Schwinger equations emerge naturally as a special sub-class of 
identities derived in Appendix [A.l| . Another sub-class of identities which may be derived 
from the functional formalism of Appendix [A.l| are thermal Ward's identities. For example, 



for the 0{N)^^ theory we immediately get from ( [A.10| ) that ©°(^(a;,y) = D^^(x,j/) and 
D"jj (a;, y) = a a ^ b. We, however, do not intend to dwell on this point more. 



We thus find that our functional formalism derived in Appendix |A.1| provides a unify- 



ing framework embracing such diverse concepts as Wick's theorem, the Dyson-Schwinger 
equations and Ward's identities. 



Appendix B 

Surface term in Eq.(4.69) 



In this appendix we give some details of the derivation of Eq. (|3.69| ). We particularly show 
that the surface integrals arisen during the transition from ( |3.68| ) to ( p.69| ) mutually cancel 
among themselves. As usual, the integrals will be evaluated for integer values of D and 
corresponding results then analytically continued to a desired (generally complex) D. 
The key quantity in question is 

1 f d°q f 2q2 \ £(g„) „ , 2, 



+ 2y(^l(ij^j^^^'<«-'"''^» 



1 /■ rf^g / 2q 



2 J (27r)^-i \{D-1 



(5+(g2-m2(0)). (B.l) 



(B.l) = ArrK(T))-ArK(0)) 

+ tl 2(ihY, I (2^-r l^^_, "^ q "(^^ - mlm) OM (^ + 1 



-2 



As usual, ab = X]j=i ^i^i ^^^ ^r ^ is a (D — 2)-sphere with the radius R. The expressions 

153 



APPENDIX B. SURFACE TERM IN EQ.(4.69) 154 

for A/t and TV are done by ( p.70| ) . 

With the relation ( [B.2| ) we can show that the surface terms cancel in the large R limit. 
Let us first observe that 



- ii^^^pEi)-; (^^^P^^(^o-^ -^.m);^^^^3T 

- i- ^^^ifpEiy y^/,,,^^^0 ^.^^ - 0. (B.3) 

In 2-nd line we have exploited Gauss's theorem and in the last line we have used L'Hopital's 
rule as the expression is in the indeterminate form 0/0. The remaining surface terms in 
(IRI ) read 



lim / -%^ / ds q {e{q' - ml{T)) - e{q' - m,^(0))} e{qo) 
1^00 J (27r)^ JdsS-^ 



The last identity follows either by applying L'Hospital's rule or by a simple transformation 
of variables which renders both integrals inside of {. . .} equal. Expressions on the last lines 
in (|B.3|) and (|B.4|) can be clearly (single-valuedly) continued to the region ReD > 1 as they 



are analytic there. We thus end up with the statement that 

(^)=UT{m',{T))-U{mlm. 



Appendix C 

High— temperature expansion of the 
gap equation 



In this appendix we shall derive the high-temperature expansion of the mass shift 5rn?{T) 
in the case when fields $„ are massive (i.e. 771^(0) 7^ 0). 

Consider Eqs.( |3.40D and ( p.49[ ). If we combine them together, we get easily the following 
transcendental equation for 5rn?{T) 

5m\T) = Ao {M{ml{T)) - M{mlm + ^him^M + 5m\T))\ . (C.l) 

Here A^ and Iq are done by ( p.42| ) and (|3.78|) , respectively. 

Now, both Ao and M. are divergent in D = 4. If we reexpress Aq in terms of A^, 
divergences must cancel, as 5m? [T] is finite in D = A. The latter can be easily seen if we 
Taylor expand ^A, i.e. 

M{ml{T)) = M{ml{0)) + 6m\T) M\ml{0)) + >l(m2(0); 6m\T)). (C.2) 



Obviously, Ai is finite in D = 4 as A^ is quadratically divergent. Inserting ( |C.2|) to ( |C.1| ) 
and employing Eq.( |3.55|) we get 
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6m\T) = K \ MimliO); 5m\T)) + ^^(^^(O) + 5m\T)) 



(C.3) 



This is sometimes referred to as the renormahsed gap equation. In order to determine M. 
we must go back to (|C.2|). From the former we read off that 



M{ml{T)-5m\T)) 



So 



MimliT)) - 7W(m2(0)) - dm'iT) M'{ml{id)) 



rfi- 



D^ 



2(47r; 



D 



D^i 



327r2 



5m^{T) = A, 



m2(T))^"i-K(0))2 
m^^{T) 



5m\T){^ - I) {mlm^-^ 



m^.{T) In 



5m\T) 



m2(0) 
K(0) + 8m\T)) In (l + ^) - 8m\T) 



327r2 




(C.4) 



(C.5) 



Eq.( |(J.5| ) was firstly obtained and numerically solved in |^. It was shown that the solution 



is double valued. The former behaviour was also observed by Abbott et al. |T8[ at T = 



and by Bardeen and Moshe ||7^ at both T = and T 7^ 0. The relevant solution is only 
that which fulfils the condition 5m? {T) -^ when T — i> 0. For such a solution it can be 
shown (c.f. 1^, FIG. 3) that "^L-* <^ 1 for a sufficiently high T. So the high-temperature 
expansion of (|C.5|) reads 



6m^{T) = Xr 



2m2(0) 



(5mHT)f (5m^{T)f 
6m4{0) "■" 12m6(0) , "'- r 

32^^2 +2 ° 






KT^ 

24 



Xrmr{T) 

8^r 



T + C m2(T)ln 



rUriT) 
TArr 



(C.6) 



Appendix D 

Derivation of the Shannon 
(information) entropy 



(This section is based on refs. p6|, [28| ,p^, |T25[] .) 



In Chapter ^ we intensively used the concept of the Shannon (or information) entropy 
which we identified (up to a negative muhiphcative constant) with an "information content" 
associated with a (macro-) system in question. The objective of the present appendix is 
to formulate mathematically more clearly the rather vague notion of "informative content" 
and to find a link between statistical physics and information theory. As a next step we 
derive the basic properties of the Shannon entropy. 

Let us first introduce the concept of information in the context of probability calculus 
and of information theory. We shall consider a set of events (or ensemble of all possible 
messages) {xi, . . . , xn} with respective probabilities {pi, . . . ,Pn}- So, for example, if only 
single letter messages are transmitted 

{xi,...,xn} = {A,B,...,Z} , 

then the corresponding set of p's characterises a particular language^]. Because events {xj} 
^For instance, p{A) in English is 0.0703, 0.0645 in French is and 0.0693 in Czech, the least frequent 
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are all possible different messages which might be send from a sender to a recipient, the 
corresponding probabilities must sum up to one: J2i Pi = ^■ 

One says that Q^ is the amount of information conveyed by the message x^ if Qm is a 
non-negative and continuous function S^ of Pm defined on the range < pm < 1- Moreover, 
the more likely a message is, the less information is conveyed by the knowledge of its 
actual occurrence (the more stereotypical a message is (i.e. the larger is the probability for 
being received), the less information it imparts). This implies that Qm = ^{Pm) niust be 
decreasing function of Pm on the interval [0, 1], and particularly '^{pm = 1) = 0. 

In order to get more qualitative results let us consider first a simple system with only two 
possible messages xi and X2 and the corresponding language {^1,^2}- Since the probability 
of receiving the amount of information 53i is pi and that of '^2 is P2, the expected averaged 
amount of information received is 

I{.Pl^P2) =Pl^iPl) +P2^iP2), 

or, more generally, for A^ mutually different messages {xi, . . . , xn} with the language {pi} 
the average amount of information received is0 

I{pi,. ..,Pn) =Pl'^{Pl) + ■■■ +PN'^{Pn) ■ (D-1) 

Let us note that ( p.lj ) implies that /(...) is symmetric and continuous in all its arguments. 
As we shall see, these conditions will strongly restrict a possible class of feasible J's. A 

further, sever restriction on the possible form of I{pi, . . . ,Pn) is imposed by the, so called, 

letter Z has p{Z) = 0.0005 in English, 0.0006 in French and 0.0008 in Czech. In information theory is usual 
to call a set of p's a langue even if a message itself is not composed of actual letters 

^We may also say that the averaged amount of received information is the averaged gain of information 
associated with the transmitted message or equivalently, the ignorance of a receiver, that is removed by 
recipe of the message. 
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additivity law for information |Q. To understand the basic features of the latter, let us 
start with A^ = 3. 

In order to determine the degree of uncertainty about the system (or equivalently, the 
expected amount of information received), we may, e.g., first evaluate the degree of un- 
certainty connected with the fact that we do not know whether Xi or Xi (complement 
of Xi) occured and add the degree of uncertainty whether the message in xi is X2 or X3. 
The amount of uncertainty due to (or information conveyed by) the first determination is 
obviously 

HPuPi) = Pi^iPi) + Pi^iPi) ■ 
The amount of uncertainty due to the second determination is 

/(Xi,X2|^l) := I{p{x2\xi),p{x3\xi)) = p{x2\xi)Q{p{x2\xi)) + p{x3\xi)Q{p{x3\xi)) . 

Using the Bayes's formula for the conditional probabilities {p{XY) = p{X)p{Y\X)) and 
the fact that ii X CiY = X then p{XY) = p{X) {X and Y denote sets of events) we may 
write 

/(.„ .2|x0 = ^3 ('^) + ^9 (^) . (D.2) 

p{Xi) \p{xi)J p{xi) \p{xi)J 



However (|D.2|) is the amount of information conveyed only when the Xi event occurs, so 



the total amount of information conveyed on average must be 



HPi^P2,P3) = ^(Pl,Pi) +Pi/(x2,X3|Xi) 



Hpi,Pi)+PiI{?^,?^] ■ (D.3) 

\Pi Pi, 
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We shall now show that the additivity law ( p.3| ) together with symmetry and continuity of 
/(. . .) suffice to determine J(. . .) uniquely up to a multiplicative constant. The proof will 
be done in three steps. 

Step 1 

Let us assume that we have events {xi, . . . , xat-i, yi, . . . , yu} and the corresponding 
alphabet {pi, . . . ,pm, 9i, • • • , qu}- (The way we perform a splitting into x's and y's is 
actually irrelevant.) We shall show now that the following important identity holds 



/(pi,...,pjv-i,gi,---,gM) = /(Pi,---,Piv)+P7v/ (—,•••,— ) , (D-4) 

\Pn Pn/ 



where pn = ^j=i Qi- We may think of events {xi, . . . , xn-i} as a one composite event x. 
For M = 1 Eq.(|D]5) is then trivially fulfilled. If M = 2, Eq.(^) coincides with Eq.(^) 
and so is true as well. Let us now take an induction step and let assume that (p.4|) is true 
for a general M > 2 and let us prove that this must hold also for M + 1. Actually, due to 
the induction hypothesis we may directly write for M + 1 



/(pi,...,P7v-i,gi,-.-,gA/+i) = Iipi,...,PN-i,qi,p'N) +p'nI (-r,---,^T^ 1 , (D. 

\Pn Pn J 



5) 



with p'^ = X]j=2 5'i- Using now relation for M = 2 we may write the ffist term on the RHS 
of O as 



^(Pl,---,P7V-l,gi,P^) =/(Pl,---,P7v)+P7v/ ( — , ^^ ) , (D-6) 

Pn Pn 
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with pn- = qi +p'j^ = X]i=i 1i- Using the symmetry of /(. . .) and the induction hypothesis 
we may write 



Pjv''"' Pn J KPn'pnJ Pn \p'n'"' p'n 

^ Pn Pn J \Pn Pn J \Pn Pn ' 



Plugging both ( p.6| ) into ( p.7| ) into ( p.5| ) we receive the desired relation for M + 1, and so 
this proves that Eq.( p.4| ) holds. 

Step 2 

In this step we extend our reasoning to an arbitrary number, say n, of groups of the 
messages, i.e., we want to evaluate the amount of information 

/(Pl,l, • • ■,Pi,N^,P2,1, ■ ■ ■ ,P2,N2, ■ ■ ■,Pn,l, ■ ■ ■ ,Pn,N„) ■ 

Using the result (p.4|) we may directly write the former as 

r/ \ , T fPn,l Pn,Nn\ 

^Pl,l, ■ ■ ■ ,Pl,Nr, ■ ■ ■ ,Pn-l,l, ■ ■ ■ ,Pn-l,N„_r,Pn> + Pn 1 I ,•••, I , 

\ Pn Pn J 

where Pn = ^i=iPn,i- If we now shift pn to the very left in /(. . .) and iterate further we 



get 



^(Pn,Pl,l, • • • ,Pl,iVi, . . . ,Pn-l,l, • • • ,Pn-l,iV„_i) = I{Pn,Pl,l, ■ ■ ■ , Pn-2,7V„_2, Pn-l) 

Pn-1,1 Pn-l,Nn-i 



+ Pn-ll 



Pn-l ' ' Pn~l 



with Pn-l = X]i=i ^ Pn-i,i- Shifting Pn-i to the very left and repeating iteration we get 
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I{Pn,Pn-UPl,U ■ ■ ■ ,Pl,7Vi, • • • , Pn-2,iV„_2) = HPn, Pn-l, Pl,l, ■ ■ ■ , Pn-3,Nn--iPn~2) 

Pn~2,l Pn-2,N„-2 



+ Pn-2l 



Pn-2 ' " ' ' Pn-2 



with Pn-2 = J2i=i ^ Pn-2,i- Shifting Pn-2 to the very left and repeating iteration, etc. up to 
Pi we finally obtain 

^(Pi,i)---)Pn,Jv„) =/(pi,P2,...,Pn) + ^Pi/ I ^^,...,^^) • (D.8) 

i^l \ Pi Pi J 

Eq.( p.8|) is the desired result. 

Step 3 

In this last step we shall actually solve Eq.( p.8| ). Before we start let us define one useful 
function. If all the messages were equiprobable, i.e., all pi = 1/n with n being the number 
of all possible messages, then we define 

a{n):=l(-,. ..,-]; n > 2; /(I) = . 
\n n J 

The first noticeable fact about (y{x) is that it fulfils a very simple functional relation, namely 



\vnn van ) \n nj ^—^ n \m m ^ 

= cr{n) + a{m), (D.9) 



where in the first line we have used Eq. (p.8| ). The functional identity ( p.9|) has the well 



known solution]^; a{x) = k\n{x). The constant k is the only ambiguity which the solution 
possesses. We shall specify k latter on. 



A simple way how to solve (|D.9|) is to assume that (j{x) is a continuous function. If we were able to 
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Let us now assume that probabilities {pi, . . . ,pn} on the RHS of ( p.8| ) are rational 
numbers, i.e., we may write pi = Ni/N, p2 = N2/N, . . . , pn = Nn/N with A^ = J2i=i ^i- 
In this case 

I{puP2,...,Pn) =^l ^'•••'^ 



Inserting this back into Eq.( p.8|) we may rewrite (p.8|) as 



n 






N' ' N ^^ N 

i=l 






As /(...) is by assumption continuous function we may analytically continue result (p.lOj ) 
to irrational probabilities, and so generally 

n 

I{pi,...,Pn) = -k^Pi\n{p^). (D-11) 

1=1 

Note that the structure of Eq. (p.ll|) is precisely that as in ( p.ll) . Comparing both (p.l| ) 



and (p.ll| ) we may identify the amount of information ^rn of the message Xm as ^{pn 



find the solution of the relation u{xy) — <j{x) +<y{y) for continuous arguments we could at the end restrict 
our attention only to discrete ones. Assuming this continuity we may perform derivation w.r.t. x and get 

d(T{yx) , , 

— ^— = ya (yx) ^ a (x) . 

Setting yx = z we obtain the differential equation 

z <t' {z) — X u' {x) — k , 

where k is constant. So the solution is obvious: (t(x) — fcln(a;). 
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—kln{p„i)- Because 3(. . .) is a non-negative, decreasing function of its argument, the 
constant k must be positive^ 



The constant k clearly fixes the units of information. In information theory k is chosen 
to be 1 with the logarithm in the base 2, so /(. . .) = — X]jPi log2(Pi)- The reason for choice 
this is quite pragmatic. The messages are usually written in the binary characters ({0, 1}, 
yes/no). The most elementary message is composed of two equiprobable events (unbiased 
choice between two possible messages about which a receiver does not posses any further 
information). In this case the information entropy is 1(1/2,1/2) = cr(l/2) = 1. Thus 
the most elementary message carries the unit amount of information. This unit is called 
"bit" f\. If one transmits only ra-letter messages (i.e. Ill ... 11, 111 ... 10, . . . , 000 . . . 00) 
then the amount of information is clearly o"(2") = n. Consequently such messages convey 
information with n bits. 

As it was mentioned, in information theory Eq. ( p.ll|) refers to the situation after the 
reception of a message. Accordingly we interpret ( p.ll|) as the average gain in the infor- 
mation associated with the transmitted message (the greater the initial uncertainty, the 
greater the amount of information conveyed. If there is no initial uncertainty or doubt 
to be resolved, the alphabet {pi} shrinks to a single case {pj = l,Pkj^j = 0} and hence 
the Shannon entropy (gained information) is zero). Let us note that if events {xi} are 
equiprobable the information entropy equals to the expected number of binary (yes/no) 
questions whose answer take us from our current knowledge to the one of certainty^. 



"^Note: If fc = log2e the amount of information (D.ll) is called both the Shannon and information 
entropy. For a different choice of k only notion of the Shannon entropy is used. 

'^If the base of the logarithm is e {k ~ 1), then the information is measured in "nats". 

^For instance, for /(1/2, 1/2) the binary question may sound: is it 1 which is transmitted ? In the case 
of the n-letter messages we may ask whether the transmitted message has on the first position 1 (yes/no), 
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In physical systems, however, no message is send, so to speak. We may, nonetheless, 
in accordance with the previous derivation view ( p. Ill ) as the mean information about 
the system which is not transmitted yet (i.e. the recipient is waiting for it). From this 
standpoint /(. . .) appears as entropy of the language, that is as uncertainty (or ignorance) 
about the ensemble of all possible messages which may be received. The only difference is 
that the uncertainty (ignorance) cannot never be completely removed, as in the case when 
the message is received, but it may be removed partially, namely when one performs a 
measurement on the system. We shall return to this point in a while. 

Now, question stands what k we should choose in real (macroscopic) physical systems. 
It is clear that in statistical physics, the systems are too complex and the amount of all 
possible transmittable information is so vast that k must be chosen very small in order for 
one to get judiciously large numerical values of the Shannon entropy in typical processes. 
Because k should be a constant valid for all systems, its numerical value depends only on 
the choice of physical units and hence may be determined via arbitrary, but suitably chosen 
system. For example, we may define as a unit entropy the entropy corresponding to one 
mole of spins of free valence electrons in a piece of iron. Assuming that all the 2^^ spin 
configurations are equiprobable (Avogadro's number A^^ = 0.6024 x 10^^ mol~^), then this 
yields /(. . .) = cr(2^'*) = kNA = lmol~ . From this reasoning we would get k ~ 10~^^. 
In order to obtain a connection with the usual von Neumann-Gibbs entropy of statistical 
physics we may note that the Shannon entropy coincides with the von Neumann-Gibbs one 
provided we fix A;/ln(2) = ke = 1.3804 x lO^^ajK-i. 

In connection with statistical physics we may define a notion of information content 

^ inherent to a system of interest. Let us assume that the system had originally the 
on the second position 1 (yes/no), ... , at the n-th position 1 (yes/no). 
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information entropy /o(. . .) (i.e., entropy of the ensemble of all possible transmittable results 
or all possible results of a measurement). Entropy Iq is often called a-priory entropy. 
Then measurement is made but because of experimental errors or impossibility to measure 
all phenomena there is a whole new ensemble of values, each of which could give rise to 
the one observed. The information entropy, say /i(. . .) may be defined also for this a- 
posteriori ensemble. The latter expresses how much uncertainty still left unresolved after 
measurement. Let us now define the informative content ■^ as an amount by which the 
uncertainty about a system has been reduced, i.e. 

or equivalently Ji = Jg — '^. After discarding the (constant) additive entropy Iq the lat- 
ter leads to the statement that the informative content is equivalent to negative entropy 
( "negatropy" ) . That is, as our information about a physical system increases, its entropy 
must decrease^ . This result is due to Szilard ||106|| and Brillouin ||105 |. 



The passage to quantum mechanics is simple. If the macro-state of a system is repre- 
sented by the density matrix 

PH = ^Pn\'4'n){i^n\, (D.12) 

n 

then the information entropy turns out to be (for simplicity we omit from now on the 
sub-index H) 

J(p) = -Tr(plog2p), (D.13) 

and the information content 

<5 = Tr(plog2p). (D.14) 



^Note that in information theory /i(. . .) = 0. 
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In the language of information theory the ensemble of all possible messages is set of all 
possible results of a measurement of a given system of observables. The corresponding 
alphabet is {pi,P2, ■ ■ ■}■ 

Some properties of the Shannon entropy 

• Concavity: The Shannon entropy is concave on the set of p's on a given Hilbert 
space. 

Concavity of I{p) means that for any pair pi and p2 and < A < 1 we have 

/(Api + (1 - A)p2) > A/(pi) + (1 - A)/(p2) . 
This may be proven very simply with a help of the inequality 

Tr(Xlog2 Y) - Tr(Xlog2 X) < (TrX - Trr)log2e . 

(use the spectral decomposition of X and Y and the inequality Inx < (1 — x)) The gener- 
alised concavity identity reads 

i{J2kp)>J2xj{p), 

i i 

where Aj > and J2i ^i = ^■ 



Maximum: If the possible kets in the spectral decomposition ( |D.12| ) span a finite 
W-dimensional subspace of the Hilbert space then 

/(p)<l0g2W, 
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with the equahty only in the case when all probabilities in ( p.l2| ) are equal, i.e. Pi = p = 
1/W. 

To prove this we may look at the difference 

I{p)-\og,W = Trip {log, p-'W)). 

Taking the spectral decomposition of p together with the inequality Inx < (1 — x) (note the 
equality is fulfilled only when a; = 1) we obtain I{p) — logaW < 0. The latter is equal to 
if and only if W/pi = 1 for all pi. 

• Minimum: The Shannon entropy has a minimum equal to zero. This happen only 
when p describes a pure state. 

Because — log2(. . .) is a convex function, one may use Jensens's inequality of statistical 
mathematics [^,0]: if /(. . .) is a convex function then (/(X)) > /((X)). Thus /(p) = 

-(log2P) > -log2(p) = 0, so 

lip) > . 

I{p) = only if there is no uncertainty about a message, i.e. when alphabet {pi} shrinks 
to a single case {pj = l,Pij^j = 0}, i.e. when p describes the pure state. 



Appendix E 

Some mathematical formulae 



E.l Integrals in D dimensions 

(This section is based on refs. [B2|,B^,E0, 



Throughout our dissertation we frequently apply dimensional regularisation; i.e. we replace 
the dimension 4 by a lower dimension D where the corresponding (loop) integrals are 
convergent. Bellow we provide a short list of integrals which we found useful during our 
calculations (cf. Sections |3^ , |0| , |3^ and ^^ , Appendices |B| and |C|). 



The paradigmatic integral of the dimensional regularisation is 



(ql + XY r(n) 

(use D-dimensional polar coordinates and the fact that ^ dVt = S^^^ = 2tt~ /T{^) ) Eu- 
clidean regime is defined via Wick's rotation as; g° = iq^, q = q^;, d^q = id^qE- Although 
the LHS as a D-dimensional integral is senseful only for integer values of D, the RHS has an 
analytic continuation for all D G C with D ^ 2n (so namely ioi D = 4: — 2e {e > 0,6 ^ 0)). 

Performing the change of variables qE ^ Qe + Ie we get 

^Note: All the quantities entering formulae bellow are dimensionless!! 
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D^ 



d^QE D/2r(^-T 



D 



(g| + 2/^g^ + X)" '^ r(n) ^^ '^^ 



Successive derivatives with respect to 1% then yield 



(g2 + 2lEqE + X)" ~ fH 






X <jr(n-f)/^/^ + -5^^r(n-i-f)(x-/|: 



1 

2' 
The analytical continuation to Minkowski regime (i.e. Wick's rotation of both q% and 1%) 

together with Eq.( ^.75| ) gives 

d°<i ^ _ r(l-f) , 2ND/2-1 



[m 



{27c)D q^ - m^ ± le (47r)^/2 



2 

m 



IGvr 
with 

A = 7 + ln47r . 

e 

(for convenience we have introduced the usual factor l/(27r)^) Previous results together 

with the Feynman parametrisation: 1/ab = J^ dt l/[at + 6(1 — t)]^, yield 

d^q i 

{27i)^ (g2 — 77^2 _|_ ie)(^(^q + p)"^ — m? + ie) 

^ 'a- [ dt\n{m^ + p^t^ - t)) + 0{e] 



167r2 

d^q i q^ 

(2,71)^ [kP' — IV? + ie){{q + p^ — tti^ + ie) 



327r2 V ^0 
d^q i q^q^ 1 



^"^ ' A- / citln(m2+p'(t'-t)) + 0(£; 



(27r)-^ (fc2 — rn? + ie){{q + p)^ — -m^ + ze) IGvr^ 



{g'"'A{p'^,m) +p''p''B{p'^,m)) 
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with A(p^,m) and B{p'^,m): 

A{p^,m) = - i -m^ |A-|lnm2 + |-/ dtln{m'^ + p^{t^ - t)) 



+ y 



A + f-/ dt\n{m^ + p\r - t)) 



+ 0{e] 



B{p^, m) 



3p2 



m 



A - ilnm^ + 3 - / c?t ln(m^ + p^{t^ - t)) 
Jo 



A + 1 - f dt\n{m^ + p\t^ - t)) +0{e)\ 



Note: the integral f^ dt {. . .) might be evaluated explicitly, the result reads 



dt\n{m^ + p\r - t)) = ln(m^)-2 + 2 



4- a 



arctan 



V4-a/ 



with a = p /m . 



E.2 Special functions and important relations 



(This section is based on refs. [§6|,|T],||,||,|7|,^,^,g].) 



The gamma function T{x) and the Riemann zeta function ^(x) are defined as follows: 



Tix) 



dt e'H 



tsX-l 



Rex > 



C(x) = ^'^ ^ Rex > 1 



n=l 



The above definitions converge only in the specified regions of the complex plane, but they 
can be analytically (single-valuedly) continued. The following important relations (used 
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in Section |3.4| ) are valid in the entire complex plane (save for the points x = —n, {n 
0, 1,2,.. .) where the simple pole residue is ~^| ) 

r(x + 1) = xT{x) 



For n being integer 



rfxird - a;) 



TV 



sin(7ra;) 



r|i + .)r(i-x 



TT 



22x-l 

r(2x) = — ^r(x)r(x+2 



cos(7ra;) 



r ( 2 ) ^'H{x) = r ( ^- ) TT— c(i - x) 



Ti--n 



-ly 



2"^/? 



(2n-l)!! 
Important numerical values of ({x) used in the text are : 



X 





3 

2 


2 


5 
2 


3 


4 


5 


C(a:) 


1 

2 


2.612 


7r2 
6 


1.341 


1.202 


90 


1.037 



(note: only numerical values of C(2?t- + 1) are available) 
Important numerical values of T{x) used in the text are : 



TT 



X 


1 

2 


1 


5 
4 


3 

2 


7 
4 


r(x) 


V^ 


1 


0.906 


iv^ 


0.919 



APPENDIX E. SOME MATHEMATICAL FORMULAE 173 

Both ({x) and T{x) appear in expansions of the following definite integrals used in the text 
{n> 0): 

/ dt = (1 - (1 ^ l)2-")r(n)C(n) (Einstein's integrals) 

Jo e* ± 1 

/ rft =2(i-2-")rHCH 

Jo smht 

"^^^ fc=0 

/ rft— -^ = 2^-'^(i-22--)rHC(n-i). 

Jo cosh t 



Euler's ip function (or digamma) was used both in Section p.4| and Section |4.4| . ip{. . .) is 
defined as the logarithmic derivation of F function: 

1 dV{x) 



T{x) dx 

p-Zt \ 

dt. 



oo /g_i ^-zt 



iIj{x) 



^ t 1-e- 
The former directly implies that 

ip{x + 1) = ip{x) H — 

X 

7A(i) = V^(l)-21n2, 
or recursively {n is integer) 

^(x + n) = ip{x + 1) + ^ + + . . . + 7 -T . 

^ ' ^ ^ 1 + x 2 + x {n-l)+x 

Defining the Euler-Mascheroni constant 7 = —^(1) (the only numerical value is known: 

7 = 0.5772156649 . . .), we get directly 

n—l ^ 

^(n) = -7 + ^ 

k=l 

^(i + n) = -7 - 21n2 + 2(1 + ^ + ... + ^^) ■ 
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In Section ^^ we use functions Kn, which are the Bessel functions of imaginary argument 
of order n (we deal only with n being integer). Kn is defined as: 

v^(ix)" '•°° 



K4x) = ^f^ /°°dte— ^*sinh^"t 
^(r^+2) Jo 

9) Jo 



T(n + 



The former implies the important relation: Kq{x) = Ki{x). In Section ^^ we use the 
following relations: 

dt — , = e 2 A I -^ 



v^t(t + a) V2 



/ dt , = ae^KAap) . 

Jo Vt^ + 2at ^ ' 

The limiting form for small arguments x (n fixed) reads: 

Kn{x) r. \T{n){\x)-'^ . 

Some miscellaneous functions used in the text: 

• Beta function B{z;y) (see Section 13.4]): 



B{x;y)= f dtr-\l-t)y-^ = f 
Jo Jo 



,0 JO (l + t)^+^ Tix + y) 



Gauss' hypergeometric functions 2-F'i[- • •] (see Section |3^) : 

pr . 1 r(c) ^ r{a + k)r(b + k) x'' 



1 /"^ 

-/ dtt''^i(l-t)^-'-^(l-tx)-" 

5 c 0) Jq 



B{h-c- 
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(note: 2F1 converges for |x| < 1 with a branch point at x = 1, for c = —n{n = 0,1,...) 
2-Fi is undetermined) 



The following relations for 2-^1 are used in Section p.4| : 



2Fi[a,6;a-6+l;-l] = 2-"v^ 



T{l + a-b) 



2Fi[a,b; |a + |6 + |; |] 



r(l + ia-6)r(i + ia) 
r(l + ia + |6) 



r(i + ia)r(l + |6) 



Bernoulli numbers -B^ (see Section p.4| ) are defined through the series : 



X 



T-^«' 



a 



a=0 



X 
^ 



Ixl < 2tt . 



Important numerical values of B^ used in the text are 



a 





1 


2 


4 


Ba 


1 


1 

2 


1 

6 


1 
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